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ABSTRACT 
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DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS 

SOME ASPECTS OF WIND TUNNEL MAGNETIC SUSPENSION SYSTEMS WITH SPECIAL 
APPLICATION AT LARGE PHYSICAL SCALES 
by Colin Paul Britcher 

Despite over a quarter of a century of development, wind tunnel 
magnetic suspension and balance systems (MSBSs) have so far failed to 
find application at the large physical scales necessary for the majority 
of whole model aerodynamic testing. Recent developments, such as the 
cryogenic operation of wind tunnels and advances in superconducting 
electromagnet technology have greatly reduced the apparent cost of a 
large, or perhaps more specifically, high Reynolds Number MSBS (LMSBS). 
Many difficulties remain, however, and three are addressed in this 
thesis. 

A powerful method of magnetic roll torque generation is essential for 
any LMSBS. Two variants of the new Spanwise Magnet scheme are studied 
herein. Spanwise Permanent Magnets are shown to be a practical method 
and are experimentally demonstrated using the Southampton University 
MSBS, though precise evaluation of maximum torque capabilities has not 
been possible. 

Extensive computations of the performance of the Spanwise Iron 
Magnet scheme indicate potentially powerful capability, limited 
principally by current electromagnet technology. Some experimental 
verification of the computed performance at low applied field levels 
is presented. 

Aerodynamic testing at extreme attitudes is shown to be practical 
in relatively conventionally configured MSBSs. Preliminary operation 
of the Southampton University MSBS over a wide range of angles of 
attack is demonstrated. 

The impact of a requirement for highly reliable operation on the 
overall architecture of LMSBSs is studied. It is shown that the 
system's cost and complexity need not be unduly increased, provided 
certain of its unique characteristics are exploited. 
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SYMBOLS 


A 


b 


C 



a 

E. 

1 

F. 

1 

F 

g(b) 

H. 

1 



(SI units used throughout, except where noted) 
x-wise coordinate of element 
transformation matrix (field components) 
wing root chordwise cross sectional area 
y-wise coordinate of element 

transformation matrix (field gradient components) 

magnetic flux density (magnitude, in tesla, direction 
and/or gradient denoted by 
subscript) 

magnetic flux density (vector) 

peak energy product for permanent magnet material 
z-wise coordinate of element 
lift coefficient 

roll moment due to sideslip coefficient (US notation) 

parameter locating E/M axis 

even function (of subscript) 

magnetic force (in subscript direction) 

magnetic force (vector) 

unknown function of b 

magnetic field strength (magnitude, in A/m, direction 
and/or gradient denoted by 
subscript) 

coercive force (A/m) 

intrinsic coercive force (A/m) 

effective coercive force (A/m, see Section 4.3) 

magnetic field strength (vector) 

current (in E/M i) 

2 

current density (in A/cm ) 

magnetic polarization (magnitude, in tesla, direction 
denoted by subscript) 

current density in E/Ms generating through-wing field 
current density in E/Ms generating wing magnetizing field 
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• J and J 
r m 

remanent induction (Section 4 only) 
magnetic polarization (vector) 
roll torque (balance axes) 

roll moment due to sideslip coefficient (UK notation) 
mach number 

magnetization (magnitude, in tesla, direction denoted 

by subscript) THIS IS NOT TO SI STANDARD 

magnetization (vector) 

magnetic dipole moment (M x <$V) 

magnetization (traditional, in A/m, see Appendix 4) 

number of E/Ms in an array 

yaw torque (balance axes) 

odd function (of subscript) 

position vector of element 

wing semi-span 

wing area 

magnetic torque (vector) 
volume of incremental element 

volume of whole core or non-incremental element thereof 
arbitrary vector 

balance axis or model axis (fixed model or model only) 

model axis (model free in balance axes) 

balance axis or model axis (as above) 

model axis (as above) 

balance axis or model axis (as above) 

model axis (as above) 

angle between magnetic force vector (in xz plane) and 
x axis 

incremental element of i 
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pitch attitude (Euler angles) 
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SUBSCRIPT NOTATION 

x,y,z component in x,y,z direction 

x' ,y' ,z' component in x^y'/z' direction 

o at origin of balance axes (field properties only) 

t value at wing tip 

Magnetic subscript notation dealt with fully in Appendix 4 

SUPERSCRIPT NOTATION 

' in model axes 

p peak value 


wing sweepback angle 

pitch transformation matrix (model axes to balance axes) 

magnetic permeability (used in Section 5 to represent y^) 

-7 -1 

magnetic permeability of free space (41T x 10 Hm ) 
relative magnetic permeability 
parameter locating E/M axis 
roll attitude (Euler angles) 

roll transformation matrix (model axes to balance axes) 
magnetic susceptibility 
yaw attitude (Euler angles) 

yaw transformation matrix (model axes to balance axes) 
pitch torque (balance axes) 
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ABBREVIATIONS 


A.C. alternating or unsteady current 

A/D analogue to digital (converter) 

APO Automatic Power Off 

AR wing aspect ratio 

CCD charge coupled device 

CID charge injection device 

D/A digital to analogue (converter) 

D.C. direct or quasi-steady current 

E/M electromagnet 

EPS Electromagnetic Position Sensor 

LHe liquid helium 

LN^ liquid nitrogen 

LMSBS Large Magnetic Suspension and Balance System 

MSBS Magnetic Suspension and Balance System 

N/A not applicable 

PR Passive Redundancy 

PWM pulse width modulation 

ReCo rare-earth cobalt (permanent magnet) 

RSS Rapid Shutdown Sequence 

SIM Spanwise Iron Magnet(s) 

SmCOg samarium cobalt (permanent magnet) 

SPM Spanwise Permanent Magnet(s) 

SUMSBS Southampton University Magnetic Suspension and Balance 

System 

t/c wing thickness to chord ratio 

+, X, L, V etc.classical representations of MSBS geometries 
(see Section 1.2.1) 
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INTRODUCTION 


1.1 Magnetic Suspension and Balance Systems (MSBS) with emphasis 
on Large Systems (LMSBS) 

The first MSBS for wind tunnel application was constructed in 
Prance by ONERA and first reported by Tournier and Laurenceau in 
1957 (1) . This system successfully demonstrated 5 component 
control of a range of simple models in work such as wake and drag 
studies up to intermediate supersonic Mach numbers. The potential 
advantages of wind tunnel testing without mechanical model supports, 
including the elimination of support interference and more straight¬ 
forward dynamic testing, fired the imagination of researchers in 
numerous institutions and quite soon many systems conceptually 
similar to the early ONERA system were under construction. 

A total of eleven MSBSs of various designs could be considered as 
constituting this first wave of effort and are tabulated briefly 
below. All bar one of the institutions involved presented papers 
at the First International Symposium on MSBSs in 1966 (2) . 


Table 1.1 

MSBSs pre-1970 

(listing in 

alphabetical order) 

Institution 

Degrees of 
freedom 
controlled 

Approx, test 
section size 

/ T 

a/ cm 

Application 

AEDC/NASA Langley 

5 

30 

Wake studies /R&D 

Univ. Michigan 

1 

4.5 

Low Re. No. sphere 
drag 

MIT (A) 

5 

10 

Static/dynamic 

MIT (B) 

5/6 

9-15 

Static/dynamic/ 

R&D/Magnus 

NASA Langley 

1 

11 

R&D 

ONERA (A) 

5 

5.5 - 8.5 

Drag/base pressure 

ONERA (B) 

5 

26 

Base pressure/heat 
transfer 

Princeton Univ. 

3 

13 

Wake studies 

R.A.E. Farnborough 

5 

18 

Sting effect/Magnus 

Univ. Southampton 

6 

13 - 18 

Static/dynamic/ 

R&D/Magnus 

Univ. Virginia (A) 

3 

9 

Cone & sphere drag 


(data from Ref. 3 and various) 


- 1 - 



Two further MSBSs were commissioned in the early 1970's, 
including the first system utilising superconducting electro¬ 
magnets. A Second International Symposium was held in 1971 (4) . 

_ Table 1.2 MSBSs post-1970 ______ 

Oxford Univ. 3 8 Low density sphere drag 

Univ. Virginia (B) 3 13.5 Dynamic stability/R&D 

Contemporary with these developments there was an increasing 
recognition in the aerodynamic testing community that the majority 
of existing wind tunnels were of too low a Reynolds number capability 
for the test purposes to which they were put. Further, it was clear 
that the bulk of production testing over the coming years would be 
at Mach numbers no higher than low supersonic ( <2.5), thus 
including the regimes where Reynolds number and support interference 
tend to have their greatest influence. Interest in MSBSs waned, 
due principally, it is thought, to the realisation that scaling of 
existing technology MSBSs to the large wind tunnels then known to be 
necessary, would be impractical for one reason or another, for 
instance E/M power requirements, with associated high capital and 
running costs, or lack of adequate roll torque capability. Indeed 
the modern day resurrection of the cryogenic wind tunnel concept, 
now universally accepted as the best approach to high Reynolds number 
testing in the crucial transonic regime, came about as a part of 
efforts to reduce the size of wind tunnels necessary for a high 
Reynolds number MSBS. 

Relatively recently however, various factors have acted to increase 
the apparent viability of a LMSBS. Concern over the accuracy and 
validity of support interference corrections is now widespread, 
particularly and rather paradoxically in the case of the new breed 
of high Reynolds number cryogenic wind tunnels, since these tunnels 
are generally operated at high pressures. 

MSBSs of the type considered to date are inherently capable of 
dynamic testing, which otherwise presents severe difficulties with 
mechanical model supports, and several small systems have been 
successfully used for this purpose. 

A related capability is that fixed model attitudes may be 
selected or changed rapidly and with little restriction, potentially 
manifesting a highly productive system. 
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Further potential advantages of a LMSBS are thought to lie in 
the areas of high angle of attack (extreme attitude) testing, where 
support interference can become rather intractable, and as an 
alternative solution to the difficulties involved with operation of 
complex conventional strain gauge force balances in the low 
temperature environment of the cryogenic wind tunnel. 

The above factors, coupled with the recent rapid developments 
of large scale superconducting E/Ms for power generation and nuclear 
physics research indicated that revived effort aimed at developing 
a LMSBS would be worthwhile. This effort was and is principally 
directed by NASA Langley Research Center with the intention of 
developing a LMSBS at that institution. 

1.2 Status of MSBS development circa 1978 

Many of the aforementioned small MSBSs were in no way intended to 
be representative of large systems, rather being specialised designs 
tackling specific and specialised problems. Only two systems are 
known to have been constructed with scaling to larger systems 
specifically in mind, namely MIT(B) and UVa(B). Nevertheless much 
existing technology was thought to be applicable to a future LMSBS 
and some of the relevant design features will be reviewed here. 

The reader must carefully note that the comments below refer only 
to the state of developments at the time that this project commenced 
(late 1978). Considerable developments are reported in later 
chapters herein and it must be assumed that corresponding or 
complementary developments have taken place elsewhere. 

1 .2.1 E/M configurations 

Most current MSBS E/M configurations can be considered as falling 
into one of three broad categories as follows (See Fig. 1.1): 

I) Geometrically orthogonal simple E/Ms 

e.g ONERA (A) & (B), Univ. Southampton 

II) Convoluted E/M e.g MIT(B) 

III) Orthogonal force, isotropic model e.g UVa (A) & (B) 

Class I above is characterized by separated E/Ms often of simple 
solenoidal form arranged in a relatively straightforward configura¬ 
tion around the test section, generally with the total number of 
E'/Ms equal to or only marginally greater than the number of degrees 
of freedom controlled. 
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Class IX features E/Ms, other than the axial field E/Ms, usually 
of non-solenoidal form, arranged in relatively complex ways around 
the test section with the intention of achieving relatively high 
uniformity and symmetry of all fields as well as high magnetic 
efficiency. The number of E/Ms is not strongly influenced by 
the number of degrees of freedom controlled. 

Class III is inherently limited to three degree of freedom 
control (zero magnetic torques), being originally developed to 
perform dynamic testing. This class cannot be considered viable 
for a LMSBS where principally static testing is to be performed 
since there is no possibility of selecting arbitrary model attitudes. 

Numerous subclasses of Class I exist, mainly classified by the 
E/M disposition in a cross section through the test section, such 
as "L", "V", etc. It seems certain that any future LMSBS will 
exhibit considerable E/M symmetry for various reasons, including 
E/M efficiency and system reliability, so only subclasses such as 
the classical "+" and "X" need be considered. 

Considerable efforts have been expended over many years in 
attempts to derive "optimum" E/M configurations but early design 
studies of superconducting E/Ms for a LMSBS have indicated that 
E/M configuration and geometry is likely to be forced predominantly 
by technological limits of superconductors and by requirements for 
roll control (5) . All configurations studied herein fall into 
Class I above, but further discussion of E/M configurations may be 
found in Ref. 6. 

1.2.2 Power supplies 

The class of supply almost universally used to date has been the 
phase controlled multiphase rectifier, with thyristor switching 
elements in all but the most venerable units. This type is easily 
scaled to high powers. One institution is known to have experimented 
with rotating machinery for primary power supplies (MIT) although it 
is difficult to achieve high output slewing rates with conventional 
generator designs. 

The high reactive powers that would be experienced in dynamic 
operation of LMSBSs make some form of energy storage in the supply 
most attractive, in order to raise the system's overall energy 
efficiency by regenerating power from the load E/Ms. The only MSBS 
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operated with superconducting E/Ms to date (UVa(B)), utilised power 
supplies that were specially constructed but can now be regarded as 
of relatively conventional design, a notable feature being 
capacitive load energy storage and regeneration (7) - 

Bipolar operation will be essential in at least some of the 
suspension E/Ms in any realistic LMSBS but relatively few small 
MSBSs have utilised bipolar supplies, mostly using gravitational 
or D.C» current biasing to achieve bidirectional control field 
capability. 

No power supply of the required characteristics for LMSBS 
application is known to have been demonstrated at design power 
levels remotely approaching the LMSBS requirement. 

1.2.3 Position sensors 

Optical methods of model position sensing have been employed in 
the overwhelming majority of cases. In turn, the majority of these 
methods are detail variants of the simple analogue shadow movement 
detection systems described in Appendix 1. These systems are 
inapplicable to LMSBSs because of the lack of immunity to light path 
degradation.. ONERA developed a target scanning/tracking system 
using vidicon tubes (8) , which is thought to represent a broad class 
of sensor suitable for LMSBS use, albeit with the magnetic field 
sensitive vidicons replaced with more modern solid state CCD or CID 
imaging arrays. Oxford University also used a scanning method (9) 
more closely approaching conventional T.V. techniques although the 
system was only used to detect the centroids of small spheres. 

The MIT Electromagnetic Position Sensor (EPS) (10) has often been 
viewed as a first choice for a LMSBS position sensor since the non- 
optical approach accrues, in principal, considerable flexibility and 
versatility. At the time of writing that view must be considered 
still valid since the system has amassed hundreds of hours of operation 
at small scale, although the problem of electromagnetic noise pickup 
is severe enough to demand very careful design for a large system. 
Numerous advanced optical methods have been proposed and some sub¬ 
jected to preliminary evaluation. None has so far been demonstrated 
on a comprehensive and fully representative basis at any scale. 

1.2.4 Control systems 

These have almost universally consisted of simple analogue feed- 
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back networks. Digital controllers have recently been demonstrated 
with an idealized single degree of freedom system (MIT) but no 
comprehensive application has been attempted. There does not 
appear to be serious difficulty with the control hardware for a 
LMSBS, the technology of reliable high speed digital flight control 
computers developed for civil, military and space applications 
appearing adequate. 

Existing control algorithms generally do not encompass the aero¬ 
dynamic characteristics of the model, nor the magnetic or inertial 
cross-couplings between.model degrees of freedom which exist in all 
MSBSs, rather controlling each degree of freedom separately with 
aerodynamic loads being small with respect to model mass. This 
approach is unlikely to yield the best performance from a LMSBS 
since, firstly, the aerodynamic loads may be large with respect to 
model mass and, secondly, the aerodynamic and magnetic cross-couplings 
may be severe, particularly at extreme model attitudes. More 
sophisticated control algorithms, perhaps including self-adaptation 
or self-optimisation with changing test conditions, thus seem 
necessary for effective operation of a LMSBS. No such algorithms yet 
exist. 

1.2.5 Miscellaneous 

No MSBS yet exhibits any viable hardware redundancy or utilises 
any failure tolerant features and, principally for this reason, all 
systems tend to be relatively unreliable in operation. 

Several methods of roll control exist but all are unsatisfactorily 
weak in torque capability. Thus, whereas the generation of adequate 
forces and torques in degrees of freedom other than roll appears 
practical at large scales, the generation of adequate roll torque 
does not. 

1.3 Contents of this thesis 

The underlying emphasis is towards the concept of a LMSBS though 
the majority of contents are unspecific in their application. At 
the commencement of the current phase of investigations, there 
appeared to be several specific and fundamental obstacles to be 
surmounted, or shown to be surmountable, before detail design of 
a LMSBS could commence. These include (not in any particular order 
of importance): 



1) Generation of adequate roll torque. 

2) Provision of adequate overall system 

reliability and integrity. 

3) Design of very large capacity power amplifiers 

(loosely referred to in the context of 

MSBSs as power supplies). 

4) Design of large multipurpose A.C. capable super¬ 

conducting E/M arrays. 

5) Development of large scale, precise and versatile 

position sensors. 

6) Development of advanced control algorithms including 

full magnetic decoupling and some self-adaptive 

features. 

Of course only a few of these topics could be addressed within 
the scope of a study of this kind. 3) and 4) above are quite 
beyond the capability and experience of this University. It had 
been felt initially that 6) above would represent a worthwhile and 
significant early advance but this view was quickly modified. 

Other problems, such as l)-5) quoted, were recognised as representing 
potential blocks to the development of a LMSBS, whereas advanced 
control algorithms need not be made available until quite late in 
the construction of such a system. Since it seems certain that, 
for reasons of reliability and versatility, the major part of the 
control system of a LMSBS would be implemented with some form of 
digital computer, merely realistic estimates of required computer 
power would be required in the early stages of design. 

Some attention was paid to 5) above but it was eventually 
concluded that the technological advances of certain sensing devices, 
such as solid state photo-detector arrays, currently occurring for 
application in other areas, may alleviate this problem in the 
relatively near term. 

The bulk of this thesis is therefore directed towards problems 
1) and 2) above, also the matter of extreme attitude testing in MSBSs 
which, although potentially very advantageous for a LMSBS does not 
represent a crucial block to its development. 



2. THE SOUTHAMPTON UNIVERSITY MAGNETIC SUSPENSION AND BALANCE 
SYSTEM (SUMSBS) 

2.1 Historical aspect 

Studies of magnetic suspension systems began at Southampton 
in 1959, construction of a system beginning in 1962, with the 
principal objective of performing dynamic testing in conventional 
wind tunnels without mechanical model supports. Work was under 
the direction of Dr. M. Judd, latterly and to date Dr. M.J. Goodyer. 
It was quickly shown that dynamic testing was indeed feasible 
although the high frequencies of oscillation required at small scale 
led to the development of the two-mass "tuned" model (11) . 

Extensive Magnus force testing was performed on ballistic-type 
models with reverse Magnus force, absent in corresponding sting 
supported tests, being discovered (12) . Several winged models 
(AGARD-G etc.) were tested in 6-component suspension and methods 
of rapid data acquisition were developed (13) . Studies of scaling 
of MSBSs (14) revealed difficulties that were also being identified 
by other groups of researchers, particularly that the power con¬ 
sumption of conventional copper conductor E/Ms in a MSBS scale as 
a high power of characteristic dimension. It became clear that 
systems designed for wind tunnels of realistic scale and Reynolds 
number capability, ambient temperature operation being mandatory at 
that time, would be quite out of the question, consuming megawatts 
of steady power. 

SUMSBS fell into temporary disuse but was re-commissioned in 1978 
and adapted to demonstrate operation with the University's 0.1m 
low speed cryogenic wind tunnel, the first known operation of a 
MSBS with such a tunnel (15, 16 ) . 

2.2 Modifications to SUMSBS to achieve current status 
2.2.1 E/M configuration 

As previously mentioned it seems certain that any future LMSBS 
will exhibit considerable E/M symmetry. It thus appeared desirable 
to reconstruct SUMSBS, previously of the unusual „ J_„ 

configura¬ 
tion (11) to a more representative arrangement. Subsequently, it 
was realised that extensive E/M symmetry would be advantageous to 
suspension of models at extreme attitudes (Section 7). Further, 
spanwise magnet rolling moment generation systems (Sections 4, 5, 

6) require in general at least a quadrupole field disposed in the 
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cross-sectional plane of the test section. Thus, in order to 
demonstrate such roll control systems, suspension at high angles 
of attack and representative LMSBS configurations, the previous 
E/M layout of SUMSBS was augmented and modified by the addition of 
two new "vertical" E/Ms, identical to the six "vertical" and 
"lateral" E/Ms existing, and provision of twin symmetrical (fore-and- 
aft) "drag" E/Ms, replacing the single E/M previously used (Fig. 2.1; 
Ref. 11; Appendix 1) . 

The system appears to fall into the "+" subclass of Class I as 
defined earlier but in order to accommodate spanwise magnet roll 
control, 6-component suspension takes place with the model's axes 
inclined at 45° to the horizontal, in the roll sense (Fig 2.2), 
wherein the system falls into the "X" subclass. 

2.2.2 Power supplies 

In common with the majority of contemporary MSBSs, SUMSBS utilised 
monopolar three-phase thyristor E/M power supplies. These were 
satisfactory insofar as they were relatively simple and had operated 
almost faultlessly for over a decade and a half. The low frequency 
switching inherent to the mode of operation does introduce con¬ 
siderable ripple in the output current and although this was not 
unduly troublesome in normal suspension there is an undesirable 
frequency response cutoff associated with the firing frequency 
(150Hz with U.K. mains supply). The monopolar operation presents 
obvious restrictions and the single quadrant capability of the 
supplies necessitated the installation of ballast resistors in series 
with each E/M to shorten the time constant of current decay. This 
also eases control system design, reducing the lag between E/M 
current and applied voltage, but would be quite unacceptable in a 
LMSBS due to the enormous power wastage that would result. 

Since some bipolar supplies were necessary for SUMSBS effectively 
to demonstrate suspension at high angles of attack (Section 7) it 
was decided to take the opportunity to solve many problems simul¬ 
taneously and incorporate a more sophisticated type of power supply, 
that is Pulse Width Modulated (PWM) transistor switching supplies, 
discarding the thyristor types completely. The supplies chosen 
(Appendix 1) were designed for use with industrial D.C. servomotors 
and feature the ability to store some regenerated energy from the 
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load. The switching frequency is 5kHz which permits more 

effective filtering of high frequency components of load current 

than was previously possible and renders the frequency response 

of the supplies near infinite as far as the control system is concerned. 

2.2.3 Position sensors 

The position sensing systems used in SUMSBS since its inception 
all fall into the category of simple analogue optical shadow move™ 
ment detectors and have never been capable of tracking a particularly 
wide range of model movement, perhaps * 15 degrees pitch/yaw, 
without mechanical realignment (11) . Other types of sensor exist, 
capable of tracking a wider range, notably the MIT electromagnetic 
position sensor (EPS) (lO) , and other types are proposed with 
promise of still wider capabilities such as photodiode or CID array 
target trackers (5, 17 ) . SUMSBS requires 6-component sensors 
capable of tracking ideally 90° of pitch rotation with simple models. 
Despite long deliberations and much experimentation it was concluded 
that in the context of the current research none of the more advanced 
sensors could be incorporated into SUMSBS due to various time and/or 
physical space and/or financial constraints. Thus a new set of 
traditional sensors has been constructed with widerange pitch 
tracking capability effected by arranging that the appropriate sensors 
may be rotated about the model's pitch axis, conventional fixed 
sensors being used where extreme model attitudes are not required 
(Appendix 1). The obvious should be stressed, that is that these 
sensors are quite inappropriate for use in a LMSBS and do not 
represent current thinking in this area. 

2.2.4 Control systems 

The simple all analogue control systems previously used with 
SUMSBS (11) were not capable of straightforward or realistic 
adaption to the requirements for suspension over wide ranges of 
attitude or for incorporation of more advanced control algorithms. 

The major part of the control system was thus abandoned and 
replaced by a mini-computer with local A/D and D/A conversion and 
limited analogue pre- and post-processing (Appendix 1). 
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ROL LING MOMENT GENERATION SY STEMS FOR MSBSs 


3 .1 Introduction, the need for roll control 

Historically, one of the key obstacles to the development of 
MSBSs for application to large scale wind tunnels was the 
development of a rolling moment generation system of adequate power. 
Nevertheless numerous small MSBSs have operated over many years 
without any active roll control system, indeed to date only two 
MSBSs have been demonstrated with such control (11, 18, 19) . It 
is pertinent to examine how this can be and why operation at large 
scales would be impractical without active roll control. 

Firstly, much useful testing could be, and has been, performed 
in small scale systems on bodies of revolution, notably sphere or 
missile drag measurements and Magnus force investigations, where 
active roll control was not required. Secondly, a considerable . 
percentage of the research effort expended with MSBSs has been in 
developing the systems themselves, particularly features such as 
position sensors, where suspension of models other than bodies of 
revolution is unnecessary, where fixed roll attitude was required 
it frequently proved relatively easy to induce a preferred roll 
attitude, with slight positive stiffness about that attitude, by 
arranging the model's magnetic centre to lie above its centre of 
gravity. 

The overwhelming majority of testing in a LMSBS would certainly 
involve winged aircraft models. Small "tare" roll moments are 
likely to exist due to slight magnetic, aerodynamic or mass 
asymmetries but presumably these could be countered by designing 
in a preferred roll attitude in the classical fashion. In general, 
however, whenever a winged model is tested with non-zero yaw angle 
(w.r.t. tunnel axes) significant roll moments will be generated, 
typically via the L aerodynamic derivative. In order to permit 
testing under these conditions a magnetic torque must be generated 
to oppose the net aerodynamic torque. Further, at high angles of 
attack the onset of stall is frequently asymmetric, producing steady 
roll torque of either sign or perhaps oscillatory torque, such as 
with asymmetric vortex shedding. Damping of roll oscillations is 
essential under all the above-mentioned conditions and is best 
introduced by appropriate manipulation of a method of positive roll 
control. 
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It is difficult to estimate the magnitudes of torque that 
might be met in a future LMSBS because the aerodynamic character¬ 
istics of the models are of course unknown. However, Ref. 20 
includes the following data for an F-94A : 

C (U.S. equivalent notation for L ) in the range 

Lg V 

-0.07 (high speed) to -0.05 (low speed) 

Span = 11.37m Wing area = 22.2m 2 

Scaling to a 4-foot span model in an atmospheric 

pressure tunnel at M = 0.85 : 

L - 11 Nm/degree Yaw 

Data from Ref. 21 relating to an arbitrary design resembling a 
Hunter gives a somewhat larger result (43 Nm/deg). It is clear 
that considerable roll torques are to be expected from this source. 

In the case of asymmetric stall we may take an arbitrary case of 
an aircraft with one wing completely stalled and the other com¬ 
pletely unstalled. A stalled of approximately unity might be 
expected with the unstalled perhaps twice this value. Without 
reference to any specific geometry or test conditions it is 
immediately seen that in the (worst) case of an untapered uniformly 
loaded wing the resultant lift vector will act some one-sixth of 
the semi-span outboard of the aircraft centreline along the unstalled 
panel. 

There thus seems little possibility of realistic operation of 
a LMSBS without a powerful method of magnetic roll torque generation. 

3.2 Methods of magnetic roll torque generation 
3.2.1. Introduction 

Various methods of generating roll torque, falling into three 
main categories, have been proposed and evaluated over a number of 


years. These include: 

D.C. field systems 

1) "Bent" fuselage core (22) 

2) Shaped fuselage core (11, 23) 

3) Through wing magnetized wing 

cores (11, 18) 

4) Active model mounted coils (11) 

5) Passive model mounted coils (11) 

6) Spanwise magnets (24) 
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A.C. systems 

7) Planar conducting loop (19, 22) 

Aerodynamic systems 

8) Active aileron control (11) 

9) L (C T ) control (11) 

v s 

Types 4, 5, 8, 9 are currently considered unsuitable for LMSBS 
application and types 1, 2, 3 have shown inadequate torque capa¬ 
bility in previous analyses. Type 7, the A.C. system, promises 
high torques but significant complication of overall MSBS design 
seems necessary to accommodate the strong A.C. fields required. 

An investigation of the two subclasses of type 6 is presented 
later. No studies of the current status of other methods have been 
undertaken so direct comparison between methods cannot be included. 

3.2.2. D.C. methods 

3.2.2.1 Theory of D.C. torque production 

At any point in a magnetized body, where magnetizations and 
applied fields are quasi-steady, the forces and torques acting on an 
incremental element of material can be written as follows : 

6f = (M . V) h Sv 6t = (M x H) Sv 

If M represents permanent magnetization (polarization) there is 
no difficulty and H may be taken directly as the external or applied 
H field., Where M is composed wholly or partly of induced magnet¬ 
ization there is some ambiguity since M will be a function of applied 
Hi If the total force and torque on an isolated body are required 
then it is valid to consider the induced magnetization as permanent 
and the effective H field as that field that would exist if the body 
in question were removed. A body of magnetically soft material will 
require a low demagnetizing factor in at least one direction for the 
product (M x H) to be large. 

Resolving using the co-ordinate system of Pig 3.1 : 

r = a£ + bj_ + ck 

Fig. 3.1 

Axis system for magnetic 
calculations 
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we have : 
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Integrating roll torque over the volume of the core 


(3.1) 


E {6t + (r x fip ) ■> 


/ (M H - M H ) + (M H +MH 

vyz z y x xz y y: 

M H ) c dV 


3.2.2.2 Review of some existing D.C. methods 

Shaped fuselage cores exploit the terms M^ and H^. In 

normal suspension / M H dV and / M H dV represent the 
v x xz v x xy 

"lift" and "lateral" forces respectively and rolling moment can be 

generated essentially by redistributions of these forces over the 

volume of the core such that /(-M H c + M H b) dV is 

v x xy x xz 

non-zero. The core cannot be axisymmetric to satisfy this condition. 

Drawbacks of the method are that either (H ) or (H ) the 

xy z xz y 

vertical and lateral gradients of the conventional "sideforce" or 
"lift" fields respectively, must be of large magnitude with slender 
fuselages and that valuable fuselage core volume must normally be 
omitted to satisfy the condition of non-axisymmetry. 

Bent fuselage cores introduce transverse components of magnet¬ 
ization, M^. Principal variations of the method proposed to 

date seek to exploit the terms M H or M H by arranging suitable 
y z z y 

distributions of transverse magnetization and applied fields to make 
the relevant integrals non-zero, for example the E coil method (22). 
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Large transverse magnetization components will be required to 
generate high torques but this can only be achieved in general by 
severe camber of the core magnetization or by use of more complex 
multiple reverse camber core geometries which in turn necessitate 
rather awkward field distributions. 

The previous wing mounted magnet scheme used at Southampton 

University introduces asymmetrically disposed vertical components 

of magnetization, which, with a suitable applied field distribution, 

produces roll torque via the M H term above. Historically, the 
z y 

system was chosen to suit an existing E/M arrangement. The self 
demagnetizing factors of vertically magnetized wing cores are very 
high in all conventional aircraft geometries, making iron cores 
difficult to magnetize and permanent magnet cores susceptible to 
demagnetization. 

3.2.2.3 General comments 

Moment generation via integration of fuselage forces (shaped cores) 
seems intuitively unattractive for conventional geometries due to the 
rather short moment arms available within a slender fuselage. 

However, there are at least three cases where this method may at 
least contribute usefully to an overall torque figure. These are : 

(1) Wide body transport. 

A wide body transport will exhibit simultaneously a high ratio 
of fuselage volume to first moment of wing area and a conservative 
test attitude/speed envelope. It may thus be argued that if 
sufficient forces and moments (apart from roll torques) can be 
generated with slender aircraft geometries then some of the 
fuselage core volume may be sacrificed, in the wide body case , 
for roll torque production. 

(2) F-15 style fighter. 

Several modern fighters are now exhibiting partial bifurcation 
of the rear fuselage to accommodate well separated twin engines. 

Such peculiarities of aerodynamic design create a non-axisymmetric 
fuselage of rather high moment of volume to volume ratio . 

Suitably exploited, therefore, such designs may yield useful roll 
torque without sacrifice of any core volume. 

(3) Underwing stores, engine pods, etc. 

These create a similar effect to that described under (2) above. 
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However, production of torque by exploitation of these character¬ 
istics could not be regarded as generally applicable for obvious 
reasons. 

Torque generation via transverse components of magnetization 
seems appropriate but attempts to introduce transverse magnetization 
in soft iron fuselage cores are hampered by the usual requirement 
for powerful axial magnetization. Permanent magnet fuselage cores 
would need to be carefully configured to avoid strong mutual de¬ 
magnetizing effects between adjacent orthogonal magnetization 
components. 

It is logical to concentrate the spanwise magnetizations in the 
wings of the model, though parts of the fuselage may be used, since 
the self-demagnetizing factors of cores of the general proportions 
of conventional wings, magnetized along their long axis, are very 
low. 

Examination of equations 3.1 indicates that cross couplings 
between degrees of freedom are likely if transverse components of 
magnetization are introduced. Cross couplings exist in all 
practical MSBS designs but it has been standard practice to arrange 
for each degree of freedom to be controllable by relatively simple 
spatial distributions of applied field with cross couplings well 
suppressed at the datum model position. It is not thought that 
elimination of all significant cross couplings is a realistic goal 
for LMSBS design since a wide variation in model geometries, 
including magnetization distribution, is expected and a modestly 
sophisticated control system should be capable of handling quite 
severe cross couplings. There is considerable advantage in E/M 
array design, however, if primary forces and moments can be 
generated by simple field distributions. 

Simple unidirectional transverse magnetization cannot be ruled 
out simply on the basis of its severe cross coupling into pitch, 
nor because it destroys the mirror symmetry of the model, but it 
will exhibit a null torque line unless vertical fields can be 
generated relatively independently in the areas of axial (fuselage) 
and transverse (wing) magnetizations. 
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The null torque line exists in uniform 
applied field gradients. 


Symmetrically disposed transverse magnetization components 
will generate roll torque with application of a vertical field 
with transverse gradient as Fig. 3.3. 
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There are two subclasses of this system which will be dealt 
with separately and are referred to here as the Spanwise Permanent 
Magnet (SPM) and the Spanwise Iron Magnet (SIM) schemes. 

3.2.3. A.C. Methods 

The planar conducting loop method proposed by Stephens and 
since successfully demonstrated at M.I.T. is an ingenious approach 
to the problem and does have the advantage of minimal coupling 
between roll and other degrees of freedom. However, the strong 
A.C. fields required would be difficult to generate in a large 
scale system and may not be possible at all with current technology 
superconducting E/Ms. The predictions of roll torque capability 
appear promising but it is felt that the practical difficulties of 
operating an A.C. system at large scales would lead LMSBS designers 
to choose a D.C. system, if one were available with sufficient 
torque capability, in preference to the existing A.C. method. 
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SPANWISE PERMANENT MAGNETS (SPMs) 


4.1 Introduction 

If transverse magnetization components are generated by 
installation of permanent magnet cores, predominantly in the 
model's wings, then roll torque is directly generated by appli¬ 
cation of the appropriate through-wing fields as Fig. 3.3. 

In all good permanent magnet materials at low applied fields 
the polarization (J) is practically independent of applied field 
strengths and may also be considered uniform over the volume of 
the core in the case of rare-earth cobalt (ReCo) materials, this 
class of permanent magnet exhibiting by far the highest known 
values of energy product and coercive force (25). A serious 
disadvantage of all known materials suitable for SPM application 
is poor mechanical properties, particularly high brittleness, 
that would preclude the use of permanent magnet cores as major 
load bearing members in the model's structure. Only some fraction 
of the available wing volume could thus be utilised for the SPM 
cores. 


4.2 Theory 

4.2.1 Ideal applied fields and simple spanwise magnetization 

Where the transverse magnetizations are entirely in the y 
direction (spanwise) then equations 3.1 reduce to: 

F = / J H dV 


If it is assumed that the applied field gradients are constant 
over the volume of the core then since: 


/ J dV = 0 


(by symmetry about the xz plane) 
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equations 4.1 collapse to: 

F = F = F =M = 0 
x y z 

L = / J (H + H b) dV 

v y z yz 

N = /-J (H + H b) dV 

v y X xy 

further,if it is assumed that applied fields exhibiting large 
components, (axial, magnetizing, lift) will be symmetrical about 
the xz plane (Appendix 5), then: 

/ J H dV - 0 
y x 
v 1 

Thus we have: 

F = F = F =M=0 
x y z 

L = / J y (H z + H yz ,b) dV - (4.2) 

N * /-J H b dV 

v y xy 

It is immediately seen that, since H is the principal applied 
field gradient for sideforce generation via conventional axial 
(fuselage) magnetization, there exists one fundamental cross coupling 
due to the presence of the spanwise magnetizations, that is coupling 
into yaw from a sideforce demand. It is not thought that this 
presents any serious difficulty since the magnitudes of aerodynamic 
forces and moments in the lateral plane are generally considerably 
smaller than in the vertical, whereas the magnetic force and moment 
capability due to the axial magnetization will often be approximately 
equal in' these two planes. 

Examining a SPM wing core with ideal applied fields we have (see 
Fig. 4.1) : 

From Eqns. 4.2: 

s 

L = / J Ag(b)(H + H b) db - (4.3) 

b=-s Y y 

Now H may be written as H b where H is the value of H at the 
z z t s 2 t 

tip. Assuming symmetry about the xz plane: 

L = 2 / J A g(b) H — db 
y o z s 

o t 
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If may be considered constant: 



The integral is the first moment of volume of the core and depends 
only on its geometry. is material dependent but does not vary 

dramatically amongst the best permanent magnet materials. The 

maximum usable value of H depends on material properties and the 
z t 

internal demagnetizing influences of the core and can vary 
considerably between different materials and core geometries. 

If the cross sectional area of the core is constant with b then 
g(b) =1 and we have: 

L = 2H z J yV - < 4 * 5 > 
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equations 4.1 become: 

F = / 0 / 0 / 0 dz dy dx = O 

x x y z 

x y 1 z 

F = / E /E / 0 dz dy dx = O 

y x y z 

x y • z 

F = / E /o / E dz dy dx = O 

z x y z 

x y J z 

L = r E r E r E dz dy dx (J H term) O 

X X y y z 2 yZ 

+ / E / E / E dz dy dx (J H b term) f O 

x X y y z Z y yz 

- / E / E / E dz dy dx (J H c term) f 0 

x x y y z z y yy 

M = / °X ^ °y / E z dz dy dx (J y H xy c term) = 0 

" / ° x / O y / E z dz dy dx (J y H yz a term) = O 
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N = - / O / E / O dz dy dx (J H term) = O 

x y J z 1 y x 

x y z 

+ f O / E / O dz dy dx (J H a term) - O 

x J y z * y yy 

x y z 

- / O / E / O dz dy dx (J H b term) = O 

J x J y J z yxy 

x y z 

Thus equations 4.1 collapse to: 


Of course c is small (Fig 4.1), thus equations 4.10 reduce 
further to: 

F = F = F = M = N = O 


It is immediately seen that there exist no primary couplings due to 
the transverse magnetizations with typical applied roll fields. 

There may, however, exist couplings due to the applied roll field 
with the conventional axial (fuselage) magnetization or due to 
other applied fields and the transverse (wing) magnetizations. 
Couplings between the applied roll field and conventional axial 
magnetization will be small with all realistic geometries, since the 
roll field exhibits a null along the x axis (Fig. 4.2), and will be 
zero with axisymmetric fuselage cores. 

The full calculations as shown above will not be carried through 
here for all cases but the important results are as follows: 

Table 4.1 Characteristic symmetries of applied fields 
(Appendix 5) 


Field classification 


Primary field or 
field gradient 
component 


Field component charac¬ 
teristic symmetry matrix 
(as Eqn. 4.8) 


0/ 
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Lift 


° / 


O O / 

From equations 4.1 etc. the following conditions can be identified: 

F ¥ 0 where H = 0 E 

x xy y z 

F ¥ 0 where H = OE 

y yy y z 

F ¥ O where H = O E 

z yz y z 

L ¥ 0 where H 


or 


H 


¥ 0 where H 


O where H 


yy 

or H 

xy 

These criteria can be related to the symmetry matrices in 
Table 4.1 either directly or by using: 



H = E (O )=> (implies) H 


t (O v ) and H. = E k (O k ) 

VV 


The conditions above become: 


& 0 with the applied 

field symmetry 
matrix as: 


'X E E 
X O E 
XXX 

- Satisfied by magnetizing and drag fields only 


X = don't care 


f O with / X X X 1 

X E E 
\X X X/ 

- Satisfied by yaw and sideforce fields only 

f O with / X X X 1 

X O O 
\X E Ej 

- Satisfied by pitch and lift fields only 


L / 0 with 

' | 

'x X X^ 
XXX, 

l ^ l 

lx X x| 
x E O ■ 

l ° r | 

/x X X 1 

X E O 


,X ° E ) 


[ X 0 El 


[x X X 1 


- Satisfied by no fields 
M ?0 with 


X O O I or 
XXX 


/X X X \ 
X O O 
X E E / 


- Satisfied by pitch and lift fields only 
N / O with 


fx O E ! 


lx X xl 


|x 0 E \ 

XXX 

1 or 1 

X E E 


X E E 

X X X I 


[x X xj 


( X x xj 


- Satisfied by yaw and sideforce fields only 


The origin of many of these cross couplings is easily identified 
as core asymmetry in the yz plane. If this condition is abandoned 
and a fully symmetric core considered then further couplings dis¬ 
appear, leaving: 
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with magnetizing field 


F f O 

x 

F^ f O with yaw field 

F^ O with pitch field 

L =0 with all fields 

M ^ O with lift field 

N f O with sideforce field 

4.2.3 Theory. General spanwise magnetization 



F ig. 4.3 Sweptback SPM core 

In the majority of cases of interest the model's wings will be 
sweptback as shown in Fig. 4.3 above. To preserve low demagne¬ 
tizing factors for the SPM cores the direction of magnetization 
must be sweptback also, perhaps to lie along the axis of centroids 
of core cross-sectional area. 

Full analysis of this case becomes complex since both and 
components exist. However, if the SPM cores are assumed slender, 
whence applied fields may be considered uniform over any core 
cross-section, the fundamental effects of the x-wise magnetization 
components can be identified: 
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J A Q g(b) Sin ? 
J A g(b) Cos Z> 


From Eqns. 3.1: 

L = / J H +(JH + JH) 

v y z x xz y yz' 

L = / J (H + b H - cH)+J(bH - c H ) dV 

v y z yz yy x xz xy 

The terms may be broken down as before, using Eqns. 4.9: 

( dz dy dx / 0 


But since c is typically small, we may write: 
L — f J (H + b H ) + J v b(H_) dV 


yz 


xy 


The J term is identified as the sum of the moments of the "lift" 
x 

force (F^) components over the wing volume, the terms being as in 
Eqn. 4.11. 


4,3 Theoretical performance 

Direct evaluation of the integrals in Eqns. 4.1 etc. is not 
generally possible since field and field gradient components can 
seldom be represented as simple functions of x, y and z. Where 
the MSBS E/Ms are air cored, the applied field at any point within 
the model core may be calculated analytically or numerically 
(Appendix 2), permitting evaluation of the total forces and torques 
acting on the model, where the magnetization distribution is known, 
effectively by numerical integration of Eqns. 3.1 over the core 
volume. The magnetization distribution may be assumed known in 
the case of ReCo materials below the onset of irreversible de¬ 
magnetization, then typically being constant and uniform, otherwise 
tending to be difficult to compute. 

A key requirement for LMSBS application is the identification of 
the absolute maximum torque capability of a particular core, which 
in the SPM case will be determined by the onset of irreversible 
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demagnetization of the core. 

The total field within the core is the sum of the applied (E/M) 
field and the internal self-demagnetizing field of the core. The 
self-demagnetizing field of cores of arbitrary shape and known 
magnetization distribution may be calculated numerically and certain 
idealised forms yield to direct analysis. Where the total field 
acts parallel to the local direction of magnetization the onset of 
irreversible demagnetization may be predicted directly from the 
materials BH (or JH) curve, ReCo materials generally exhibiting 
near straight-line curves in the 2nd quadrant (Fig. 4.4). 



Unfortunately, the total field will seldom act as described 
above, and whereas some data exists concerning the demagnetization 
of single crystals of certain materials with total fields not 
parallel to the direction of magnetization (for instance see Refs. 
26, 27, 28), it is not thought that this can be applied to bulk 
materials, demagnetization then being powerfully affected by 
processes principally dependent on the material's microstructure, 
such as domain wall movement (28, 29, 30). No relevant data 


- 33 - 



concerning bulk ReCo materials has been discovered (see Ref. 

31 for treatment of particular ferrites), rendering predictions 
of maximum torque capability of general SPM systems impossible at 
present. However, a crude estimate of the order of magnitude of 
available torque may be made by the following argument. 

In slender cores magnetized along their long axis, the internal 
self-demagnetizing field is small over the majority of the volume 
and may be neglected here. The effects of the magnetization 
reversal at the wing root (Fig. 4.1) and of the detail geometry of 
the cores is beyond the scope of this treatment. The magnetic 
field strengths internal to the cores may thus be regarded as pre¬ 
dominantly due to the applied field. The conventional parameters 
characterising the performance of permanent magnet materials, such 

as H , H ., J , BH etc. (25), are insufficient for identification 
c ci r max 

of incipient irreversible demagnetization. Several authors resort 
to use of the parameter H fc , representing the demagnetizing field 
strength that reduces the intrinsic polarization to 90% of the 
remanent value (25, 32, 33), illustrated in Fig. 4.5. 



Fig. 4.5 Definition of 
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Values of for relevant materials are sparsely quoted, 
reflecting the arbitrariness of the parameter's definition, some 
available data being reproduced in Table 4.2 below: 


Table 4,2 Some properties of bulk ReCo materials 


Material 

J 

r 

(T) 

H 

c 

H k 

-(kA/m)— 

• 

H . 
cx 

Reference 

Sintered 

Sm Co,. 

0.85 

650 

900 

1200 

25 

(1977) 

"Commercial" 

Sm Cc> 5 

0.84 

560 

488 

800 

32 

(1971) 

Sm Co,. based 

"RECOMA 20" 

0.90 

700 

1000 

1200 

33 

(1976) 


Current material developments are expected to result in 
improvements in H c and with tending to comfortably exceed 
H . At this stage it would not appear unreasonable to assume a 
usable value of of some 800 kA/m (VM^-IT) and of 0.85 T. 

If this value of can be taken to represent the peak total field 
strength that can be applied to the material before the onset of 
irreversible demagnetization, with the working polarization taken 
as 90% of the value of quoted immediately above, then some 
progress may be made. 

The applied field is a sum of many contributory fields ("lift", 
"drag", etc.), but it should be noted that the two fields (neg¬ 
lecting roll) likely to involve the most powerful fields within 
the test section need not cause the exposure of the wing core 
to those high field strengths (Fig. 4.6). 

If it can thus be assumed that the principal component of applied 
field within the wing cores is the applied "roll" field then H 

z t 

in Eqns. 4.4, 4.5, etc. becomes H^. Taking the values assumed 
above, we have from Eqn. 4.4: 

2 4 x lO 6 s 

L " “— - / A g(b) bdb - (4.13) 

S o ° 

and from Eqn. 4.5 for rectangular slab cores: 

L ^ 1.2 x 10 6 A q s - (4.14) 
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Fig. 4.6 Components of lift and drag fields 
at the wing core 

Eqn. 4.14 corresponds to 600,000 Nm/m 3 of core. For the represen¬ 
tative SPM core studied in Section 4.3.1 below, this would indicate 
a maximum capability of approximately 500 Nm. It may be, however, 
that applied fields of sufficient strength to realise this capability 
cannot be generated with existing technology E/Ms. In fact, taking 
the roll torque calibration calculated below (Section 4.3.1), maximum 
torque would occur at around 16,000 A/cm 2 in the specified E/Ms. 

This is well beyond current E/M technology (see Section 5.8.2). 
Present apparent LMSBS roll torque requirements (5) may be 
approachable with more carefully configured E/M arrays. 

4.3.1 Performance of representative SPM configurations 

Since accurate analytic predictions of performance (below the 
onset of demagnetization) are not practical, numerical examination 
of representative SPM configurations is appropriate, the computer 
program FORCE (Appendix 2) being suitable for such studies. Some 
preliminary results are presented here. 
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to be 


The geometry of E/M array and model is chosen so as 
directly comparable to the baseline geometry used in the studies 
of the SIM system (Section 5) and readers should refer to Section 
5.4 for explanation of the evolution of these geometries. The 
wing core is chosen to occupy one half of the available cross 
section (Fig. 4.7), the aerodynamic and other loads being assumed 
as principally carried by the surrounding structure. A 
quadrupole field elegantly provides the required through-wing 
field, in turn being supplied by the appropriate quartet of E/Ms 
from the baseline SIM array (Figs. 4.8, 5.4), though the entire 
SIM array is used for examination of 360° roll capability. 

Confirmation of some of the predicted couplings between axial 
or spanwise magnetization and the applied roll or other fields is 
straightforward. Conventionally configured but arbitrarily 
dimensioned E/Ms are specified to provide lift, pitch, etc. fields 
(Figs. 4.9, 4.10). It should be noted that the overall E/M array 
is realistically proportioned although the axial and main E/Ms do 
(unintentionally) spatially conflict. 

The magnetic forces and torques generated on the model and their 
variation with simple displacements or rotations of the model are 
shown in Figs. 4.11 - 4.19. Confirmation of the predicted magnetic 
couplings due to the spanwise magnetizations can now be made. 

These are (from Section 4.2.2): 


Force/torque 

Applied field 

Figure (zero displacements) 

F 

X 

Magnetizing 

4.19 

F 

Yaw 

4.18 

y 



F 

Pitch 

4.17 

z 



L 

None 

— 

M 

Lift 

4.14 

N 

: 

Sideforce 

4.16 


The wing core's symmetry about the yz plane is conveniently 
destroyed by a displacement in the x direction. Consulting the 
relevant Figures therefore, the extra predicted couplings for a 
core lacking symmetry about the yz plane may be confirmed also: 
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Dimensions in 
metres. All core 
polarizations = 

1 Tesla 


x—t-j— 



\ v 

SIM cores (Section , 
5) this size 



Non-magnetic 2 
envelope (7.2cm ) 


Detail view on A showing direct comparison 
to baseline SIM geometry (Pig.5.5) 


Magnetic core (7.2cm ) 

Equivalent airfoil 


Fig. 4.7 Model magnetic configuration for SPM calculations 



Fig. 4.8 Roll E/M configuration for SPM calculations. 


All dimensions in metres 













Fig. 4.11 Couplings for wing in roll field 


E/Ms as Fig. 4.8, 4.10 
Model as Fig. 4.7 
Roll field defined by : 

+X 11 = “ I 12 =+I 13 = " I i4 = 1 » 459 » 240 A ( 100 ° V 01 " 2 ) 



All force and torque components not shown are exactly or 
effectively zero. 

Fig. 4.11a Fig. 4.11b 



x (m) y (m) 
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Fig. 4.13 Couplings for fuselage in roll field 


E/Ms as Fig. 4.8,4.10 
Model as Fig. 4.7 
Roll field defined by : 



Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 


Fig. 4.13a Fig. 4.13b 



x (m) y (m) 
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Fig. 4.14 Couplings for wing in lift field 


E/Ms as Fig. 4.9, 4.10 
Model as Fig. 4.7 
Lift field defined by : 

+I i =-i 3 = -i 5 =+i 7 = 1,600,000 A (1000 A/cm 2 ) 



Model rotations and torques as Fig 4.11 

All force and torque components not shown are exactly or 
effectively zero. 

Fig. 4.14a Fig. 4.14b 
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Fig, 4.15 Couplings for wing in drag field 


E/Ms as Fig. 4.9,4.10 
Model as Fig. 4.7 
Drag field defined by : 


+I g =-X 10 = 2,500,000 A (1000 A/cm ) 

1 i- I 8- Bdl n- I i4-° 


Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 


Fig. 4.15a 



Fig. 4.15b 
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Fig* 4.16 Couplings for wing in sideforce field 

E/Ms as Fig. 4.9,4.10 
Model as Fig. 4.7 
Sideforce field defined by : 

-I 2 = +I 4 = +I g = -I Q = 1,600,000 A (1000 A/cm 2 ) 



Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 
effectively zero. 



50 








Fig. 4.17 Couplings for wing in pitch field 

E/Ms as Fig. 4.9,4.10 
Model as Fig. 4.7 
Pitch field defined by : 

-I 1 =+I 3 =-I 5 = +I ? = 1,600,000 A (1000 A/cm 2 ) 



Model rotations and torques as Fig. 4.11 

All force and torque components not specified are exactly or 
effectively zero. 

M is non-zero but less than 0.003 Nm in Figs. 4.17a-£ 


Fig. 4.17a Fig. 4.17b 
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Fig. 4.18 Couplings for wing in yaw field 

E/Ms as Fig. 4.9,4.10 

Model as Fig. 4.7 

Yaw field specified by ; 

-I =+x =-i =+I = 1,600,000 A (1000 A/cm 2 ) 

2 4 6 8 



Model rotations and torques as Fig. 4.11 


All force and torque components not shown are exactly or 
effectively zero. 

Fig. 4.18a Fig. 4.18b 



0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 

x (m) y (m) 
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Fig. 4.19 Couplings for wing in magnetizing field 


E/Ms as Fig. 4.9,4.10 

Model as Fig. 4.7 

Magnetizing field defined by s 


+I g =+I 10 = 2,500,000 A (1000 A/cm ) 
*1 " X 8 and *11 - X 14 = ° 


Model rotations and torques as Fig. 4.11 

All force and torque components not shown are exactly or 


effectively zero. 

Fig. 4.19a 



Fig. 4.19b 
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Force/torque 

Applied field 

Figure (positive x 
displacement) 

F 

X 

Drag 

4.15a 

F 

y 

Sideforce 

4.16a 

F 

z 

Lift 

4.14a 

L 

None 

— 

M 

Pitch 

4.17a 

N 

Yaw 

4.18a 


It should be noted that the two couplings involving generation 
of pitching moments (from applied lift or pitch fields) are very 
weak, due to the short moment arms of wing core elements about the 
y' axis in the chosen configuration. 

The full information presented in Figs. 4.11 - 4.19 indicates 
that further couplings will occur when the model is displaced or 
rotated from its datum position, for instance a lift field is seen 
to generate pitch and yaw torques when the wing core is rolled. 

The principal result of these will be the manifestation of magnetic 
Stiffnesses (positive or negative) in various model degrees of 
freedom when the classical fields considered are applied. These 
effects will somehow disturb the magnetic stiffnesses and dampings 
generated by conventional control algorithms and may require action 
to counter adverse effects. 

However, study of the Figures reveals that these extra couplings 
occur in a relatively straightforward manner and few are of large 
magnitude. Whilst further treatment of this topic cannot be pre¬ 
sented here, it is felt that modest adjustments to and expansions 
of conventional control algorithms will restore full required 
performance. 

Since no attempt has been made to ensure representative relation¬ 
ships between the strengths of the classical applied fields, 
further interpretation of the data discussed above should be made 
with extreme caution. 

4,4 Discussion 

The performance of SPM systems is critically dependent on the 
precise magnetic properties of the core material, particularly the 
intrinsic polarization and useful coercivity (see Section 4.3). 
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There seems good reason to expect considerable improvements in 
the values of these parameters exhibited by the best available 
materials in the near future. Fig. 4.20 shows how the highest 
reported values of BH max » the classical parameter of magnetic 
performance, have risen consistently over the past decades, the 
performance of ReCo materials still being well below the theoretical 
limits at least in the respect of intrinsic coercivity (28, 34). 

The torque capabilities calculated heretofore are thus expected to 
rise as time passes within the limits of the assumptions made. 

Further, it is widely accepted that as the temperature of ReCo 
materials is lowered (below 300K, say) their magnetic properties 
tend to improve steadily. The parameters of interest here, that 
is intrinsic polarization and coercive force, are typically found 
to increase by a few per cent and many tens of per cent respec¬ 
tively, during cooling from room temperature to around 80K (28, 30,35, 
36, 37 and Fig. 4.21). The force and torque capability of a ReCo 
cored MSBS model would thus be significantly greater at the typical 
lowest operating temperatures of a cryogenic wind tunnel than at 
"normal" temperatures. 

4.5 E'/M configurations for SPM roll control 

The E/M configuration chosen so far is a fully symmetric quadrupole 
(Fig. 4.8). There is no particular reason to regard this as being 
in any way optimum. 

From Fig. 4.11 it is seen that roll torque capability falls 
rapidly with increasing roll displacement from the selected datum 
orientation. This would be a serious limitation in practice, being 
due to the inability of the chosen array to develop H , along a 
y' axis displaced by 45° in the roll sense. This problem may be 
completely overcome by including extra E/Ms in the roll array. 

These may perhaps be shared with other force/torque generation duties 
but will be regarded here as entirely specialised for simplicity. 

The complete 8 E/M array used for the SIM calculations (the 
SPM 4 E/M array being, of course, derived from this) provides roll 
torque at any roll angle, the variation of torque capability with 
angle being shown in Fig. 4.22. 

If 360 roll angle capability is deemed unnecessary, then the 4 
E/M array may be adapted to increase the angular capability, certain 
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All currents are 1,459/240 A 
(1000 A/cm 2 ) 

Current directions chosen 
to create positive torque 

at all angles „ ; _ 

Contribution of s 



Pig. 4.22 Variation of roll torque with roll displacement for 
complete SIM E/M array (SPM cores) of Fig. 5.4 
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possible adaptions being illustrated in Fig. 4.23. Comprehen- 
sive treatment cannot be attempted here but Figs. 4.24a-d illustrate 
that some improvements are possible, using the techniques of Fig. 
4.23. 


d may be positive or negative 
E, may be greater or less than 45° 



Fig. 4.23 Possible modifications of 4 E/M SPM roll 















modified quadruplet geometries. 





5. SPANWISE IRON MAGNETS (SIMs) 


5.1 Introduction 

Symmetrically disposed transverse magnetization components can 
relatively easily be induced in a magnetically soft wing core by 
application of a symmetrical field as shown below: 



Rolling moment will be generated with application of through 
wing fields as in Fig 5.2., Section 3.2.2., and Section 4. 
Rotation of the induced spanwise magnetization vectors by the 
through wing fields will be inhibited by the high demagnetizing 
factors in the through wing direction. Calculation of the 
performance of SIM systems is not straightforward, the induced 
magnetizations not being directly analytic. 
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H (positive) 



Fig. 5.2 SIM through wing field 

Thus, whilst the SPM equations may be approximately applied if 
the induced magnetizations remain predominantly spanwise and the 
value of polarization is replaced by some mean value of induced 
magnetization, that value of magnetization and the detail perfor¬ 
mance of the SIM system cannot be adequately predicted without 
further analysis. 

However, many magnetically soft materials exhibit saturation 
inductions of approximately 2 Tesla and high permeabilities to around 
that value. Operation with higher levels of spanwise magnetization 
than those achievable with SPMs (approximately 1 Tesla) is thus 
possible, leading at least to higher torques per unit through wing 
field. Further, it is assumed (Section 5.4) that magnetically soft 
materials with good mechanical properties are available, allowing the 
entire wing volume to be utilised if required, and since no irrever¬ 
sible demagnetization effects exist in the SIM case it is to be 
expected that the absolute maximum available torque for any config¬ 
uration may be greater than with SPMs. 

Since soft magnetic materials generally behave entirely isotropically, 
available material data is sufficient to permit, in principle. 
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calculation of the performance of the SIM system at all applied 
field levels. 

5.2 Solution of iron-air-conductor magnetostatic systems 

This class of problem has attracted considerable interest over 
recent years with the application of superconductors to various 
fields such as nuclear research and with the drive to improve per¬ 
formance and efficiency of electrical machines. 

General magnetostatic problems involve the solution of Poisson's 
equation or, outside current carrying conductors, Laplace's equation. 
These equations may be solved in principle by a number of methods 
including (after 38) 

1) Analogue 2) Graphical 3) Analytic 4) Numerical 

Predictions of roll torque from SIMs require treatment of 
saturation effects, leading to non-linear solutions, and are inherently 
three dimensional. In most practical cases methods 1) and 2) above 
are not able to handle the non- linear problem, indeed graphical 
methods are generally restricted to two dimensions. Analytic methods 
are available for non-linear 3D problems but only for highly restricted 
geometries of conductor and iron. Application to general problems is 
currently quite impractical. Again following (38) existing numerical 
methods in this field may be divided into four principal categories: 

1) Finite difference 2) Images 

3) Integral equations 4) Variational formulations 

although other classes of solution do exist, such as the Monte-Carlo 
method (39). 

It would appear that image methods are inapplicable to non-linear 
problems within the bounds of the present formulations. Method 4) 
above may be considered an energy method, somewhat analogous to virtual 
work methods in structural problems, whereas 1) above tackles Poisson's 
or Laplace's equation directly but both generally require that the 
complete volume of the problem is meshed with a computation grid. 

The characteristic geometry of the SIM roll system is very "open" 

(Fig. 5.3), that is a small iron region separated from the conductors 
by large air gaps. Methods 1) and 4) above would thus require meshing 
of considerable volumes of air, at least enclosing all the conductors, 
leading to large computation times. Integral equation methods 
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(3 above) need only require meshing of iron regions and thus appear 
appropriate here. 

5■3 Program GFUN 

This program, developed at the Rutherford Laboratory, Didcot, 

Oxon, U..K. since 1970 by Trowbridge et al, is an example of the use 
of integral equation methods. GFON has been applied to a wide 
range of problems over many years (38, 40, 41, 42, 43) and has 
exhibited consistently good accuracy, consequently gaining a high 
international reputation. A disadvantage inherent to the solution 
procedures employed is a somewhat awkward representation of field 
within iron regions, as the vector sum of the fields from external 
currents and induced magnetizations, expressed as constant within 
suitably shaped elements (such as tetrahedra). In high permeability 
regions the two contributions to the field may be nearly equal and 
opposite, leading to rather poor resolution of the internal field. 

GFUN in fact handles the simultaneous integral equations as a single 
matrix equation, solution of the latter generating eigen solutions for 
the internal field. Clearly, prediction of the detail geometry of 

the internal field is affected by the choice of element distributions. 
The lack of precise information concerning the internal field of the 
SIM cores is not thought to be especially critical. It must be noted 
that the nature of the solution procedure implies that the foregoing 
adverse comments do not apply outside the iron region, however para¬ 
doxical that may seem, and need not apply to the resolution of forces 
and torques provided methods of field integration over control volumes 
external to the iron are chosen. 

Acce>ss to GFUN was granted by the British Science Research Council* 
under Grant No. GR/B/3691.5. Modifications to the program were carried 
out by Simkin to permit the full symmetry existing in many of the 
required SIM cases to be exploited (reducing computation time) and to 
provide a torque integration option. The torque integration scheme 
is conceptually similar to the well established methods for force cal¬ 
culations but the fact that previous users have apparently not required 
torque information must be seen to represent a major possible source 
of systematic error in data included hereafter, indeed some diffi¬ 
culties were experienced before a consistent integration scheme could 

* Now the Science and Engineering Research Council. 
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be identified. To counter this uncertainty, various low field 
torque measurements have been made and limited cross checking under¬ 
taken (Section 5.7). Saturation of iron regions,which could not be 
achieved in corroborative experiments, should not directly influence 
the reliability of torque predictions due to the nature of the 
methods used. 

5.4 Baseline GFUN geometry and configuration 

A baseline geometry and configuration is required from which the 
effects of variations of various parameters, such as wing aspect ratio, 
may be examined. At the time of commencement of this study there 
existed no clear specific choice of geometry of either model or E/Ms, 
necessitating evolution of a baseline geometry on the following, 
somewhat arbitrary, basis. 

The favoured scale for studies of LMSBSs has been a test section of 
approximately 8 ft x 8 ft cross section (such as NTF). Calculations 
are therefore made at this scale, but the scaling of results to 
different tunnel sizes is quite straightforward (Appendix 3). The 
test section is assumed square with no corner fillets. Clearance is 
allowed around the aerodynamic cross section for structure, plenum 
chamber etc., and is chosen to be one foot (5). A similar allowance 
is made for the thermal insulation and structure surrounding each E/M 
(5). The most uncertain characteristics of the E/Ms are the winding 
shape and maximum usable current density. It appears (5) that manu¬ 
facturers prefer circular windings where possible due to reduction of 
difficulties associated with conductor stressing. An idealized 8 E/M 
configuration has been chosen (Fig. 5.4). The problem of optimizing 
the E/M array must be dealt with separately, being heavily influenced 
by particular requirements for forces and moments in other degrees of 

freedom. Maximum usable overall current density (J) for super- 

2 

conducting E/Ms varies from order 1500 A/cm for cryostable conductors 
to order 15000 A/cm^ for adiabatically stable conductors, within the 
limits of present technology. Doubts exist (5) as to whether adia¬ 
batically stable conductors could be applied to the 8 ft scale case so 
generous winding cross sections are allowed in the E/Ms, permitting 
partial saturation of the wing cores at suitably low values of J. 

High J data is included for the purpose of identifying trends at high 
field levels. The current levels in the E/Ms are represented by 

and J as defined by Figure 5.4. 
m 
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The geometry of the baseline SIM wing core is simply chosen as a 
10% thick, 10:1 aspect ratio rectangular hexahedron (rectangular "slab" 
core) spanning slightly under one half (49.2%) of the test section 
width (Fig. 5.5). 

The choice of core material (or its magnetic characteristics) 
presents some difficulty at this stage due to the necessity in a final 
design to achieve an optimum compromise between many magnetic, electrical 
and mechanical properties. The SIM wing cores are visualised as 
occupying nearly all the aerodynamic volume of the wings with only those 
detail fittings and features potentially subject to changes during test 
programs being added from non-magnetic material. It is recognised that 
this view may require qualification in the light of detailed results and 
experience in model design, but it is nonetheless clear that the core 
material requires good mechanical properties such as high yield point, 
low brittleness and high Young’s Modulus, perhaps also at low tempera¬ 
tures. High electrical resistivity may be preferred in order to 
suppress eddy current flow in the cores but is not essential to funda¬ 
mental operation of the system. Low coercive force is necessary to 
avoid torque calibration being dependent on the recent past magnetic 
history of the cores, although if cores are operated far into saturation 
this feature becomes of lesser importance. Ultra high permeability is 
not necessary since in any airgap dominated system, iron regions of 
medium and higher permeabilities tend to behave as if infinitely 
permeable. The key magnetic parameter is undoubtedly saturation 
induction. Studies of LMSBSs (5) have indicated that valuable economies 
can be made in E/M size, hence cost, by utilising fuselage core materials 
with the highest available saturation induction. It is logical that 
similar criteria should apply to the SIM case since operation of the 
core beyond the saturation point is anticipated. 

One class of materials that appears promising is the cobalt iron 
alloys, classically represented by the 50:50 Iron-Cobalt alloy 
"Permendur". These materials have not found widespread industrial 
application due to their relatively high cost but this is unlikely to 
be a problem in LMSBS applications. Indeed the cost of a fuselage 
core for an 8ft tunnel has been estimated at U.S. $ 5000 (1981 $), 

(5). This is small compared to a typical manufacturing cost of a 
large wind tunnel model. 

Magnetic and mechanical properties of material of this type are 
dependent on the precise alloying constituents, heat treatment and 
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preparation of samples but typical properties for some commercially 
available materials are shown below. 


Material 
(Trade Name) 

Permendur Vanadium Permendur 

Vacoflux 50 

Source 

Western Electric Co. Ltd. 

Vacuumschmelze 
GmbH 

Reference 

44 

44 

45 

Density 
(kg/in ) 

8300 

8300 

8150 

Initial y 

800 

800 

1000 

Saturation Induction 
(T) 

2.45 

2.4 

2.35 

Static coercivity 
(A/m) 

159 

159 

110 

Resistivity 

0.7 x 10" 6 

0.26 x 10~ 6 

0.35 x 10” 6 

Youngs Modulus 

GPa 

- 

- 

230 

Yield strength 

MPa 



400 


It should be noted that the saturation inductions, in the range 
2.35 - 2.45T (exceeding 2.5T in laboratory specimens) are substantially 
higher than the corresponding value for high purity iron of some 
2.158T (room temperature. Ref. 44). 

The peaks in the permeability versus magnetizing field strength 
curves for the above materials (Fig. 5.6) present some problems to 
GFUN. The permeability of each iron element in GFUN is assumed 
constant throughout the element and is updated at each iteration. 

Sharp rises in the permeability of elements with relatively small 
increases in the magnetizing field acting on those elements (and vice 
versa) occurring from one iteration to the next may cause local 
oscillations of the iron's interior field and consequently slow 
convergence. These potential problems may be avoided in early work 
by choosing a permeability versus magnetizing field strength charac¬ 
teristic that falls monotonically. This has been done by arbitrarily 
fitting intermediate points between the initial constant permeability 
line and the terminal saturation boundary, thus establishing the base¬ 
line characteristic shown in Fig. 5.7. A slightly conservative value 
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of 2.3T is chosen for M , . The effects of variations of some 
sat 

material properties on the torque capability of the SIM system is 
studied later. 

The baseline wing element distribution is chosen along well 
established principles. Tetrahedral elements have proven to be the 
most reliable choice for GFUN, the total number of elements is the 
largest that can be handled by the convenient batch version of the 
program and the spatial distribution of the elements is chosen such 
that the elements are relatively numerous in the region of strongest 
anticipated magnetization. This non-uniform element distribution 
also yields superior convergence. The control volume for use with 
the torque integration schemes is dimensioned such that its surfaces 
lie close to the model core, thus yielding best accuracy. 



Fig. 5.7 GFUN baseline BH curve 
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5.5 Simple theoretical approximations 
5.5.1 Torque on ellipsoids 


Some insight into the behaviour of the SIM system may perhaps 
be gained by study of some magnetostatic torque producing system 
that is analytic. Unfortunately, as previously mentioned, few such 
systems exist and few of these yield to straightforward analysis. 

One suitable example is torque on magnetically soft uniform ellipsoids 
in uniform applied fields. The general equations applying are as 
follows 



H = H 


N M 

-Ji~y 


= (1—N ) M + U H 


x o x 


Im| 


(5.1a) 

(5.1b) 


(5.2a) 

(5.2b) 

(5.3) 
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Combining (5.1) and (5.3) 


M 

X 

P H 
o x e 

1 

4 

(5.4a) 




M 

y 

y « 

° Y e 

( ^ ) - 
1 + W 

(5.4b) 

torque 

on the 

ellipsoid is as follows: 


T 

(M H 

X Y e 

- M H ) V 

y x e 

(5.5) 


Following Fig. 5.8 we have from (5.1): 
H. Cos0 = H Cos<|> - N M Cos0 


H. Sin0 = H Sin<|) 


(5.6a) 

(5.6b) 


Combining (5.6a) and (5.6b): 


1 

H g CoscJ) 

2 

H g Sirej) \ 

H. + N M 


H. + N M 

1 — 1 


1 1 -X_ 

\ / 


\ % / 


(5.7) 


This equation gives from any given H g , (j> and M. Below satura¬ 
tion M, H and X m ate functions of each other (Eqn. 5.3). Solution 
of (5.7) thus involves an iterative process using a material's 
particular BH curve. When is found, X m is known from the BH 
curve and torque can be found from (5.4) and (5.5). This procedure 
is somewhat inconvenient. Where x m is large (permeable material 
below saturation) we have HL-0 and Eqns. 5.1 collapse to: 


M „ y H Cos<J> 


(5.8a) 

(5.8b) 


and from (5.5): 


Directly, for maximum torque at any given applied field strength 
(j> = 45° and: 

T a applied field 2 (at any given $) -(5.10) 
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Where the material is far into saturation M = M , = constant, 

sat 

ana “ ^at 

"o| H i| 

Prom (5.6): 

H Cos$ _ N M 


Cos0 

Re-arranging: 


H Sind) N M 

e ^ - y Sc 

Sin8 y 


(() = Sin 


(N - N ) M Sin : 


And from (5.5): 

CT _ M sat H e Sin and using (5.11): 


(5.13) shows that(Tl 


is independent of H and has the value: 


I V/Ma: 


X Max 


sat 
2 ^o 

Equation (5.12) becomes: 


Sin 


y x' sat | + 45 

2 U H 


It should be noted that (5.16) and (5.12) break down where H is 


again: 


CosQ = 

[ H i + N y M sat j Sin8 

Cos<|> 

\ ^o / Sin(|) 
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c/a 

Fig. 5.11a Torque per unit volume (maximum) for ellipsoids. 



c/a 

Fig. 5.11b Maximum torque for ellipsoids 




The breakdown occurs where HL = O =°°). At this point, from 
(5.17): 

Tan <J) = Tan 0 

Fig. 5.9 shows the variation of the direction of induced 
magnetization, after Eqn. 5.12, for a variety of cores. It is seen 
that with a powerful applied field the induced magnetization lies 
nearly parallel to that field. 

Fig. 5.10 shows torque per unit volume for a particular core 
(Eqn. 5.13) and clearly illustrates the existence of a finite and 
approachable absolute maximum torque. For any given geometry and 
applied field direction it is seen that the torque per unit applied 
field is far from constant. In fact, at high values of <j) with 
powerful applied fields, increases in applied field strength may 
apparently result in reductions in torque. 

The variations of the maximum torque and torque per unit volume 
with ellipsoid geometry are of interest and are shown in Fig. 5.11, 
using data from Ref. 47. The significant feature here is that 
relatively slender ellipsoids are capable of generating high torques, 
perhaps implying that slender wing cores (SIMs) will remain relatively 
capable torque producers. 


5.5.2 Linear SIM theory 

The field at any point P, in the principal plane of a symmetric 


quadrupole may be written a 
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Fig. 5.12 SIM quadrupole magnetizing field 
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H , = I m (k 1 Cos2<|> + k 2 Cos4<J) +-) 

y p “ r r 

H , = -I (1 Sin24> + 1 Sin4(J> +-—) 


where k and 1 are functions of r. For an identical quadrupole 
r r 

displaced by 45° about the x axis (giving the through-wing field 


H , = l r (k ' Sin2(J) + k 2 Sin4<f> +-—) 

y p r r r 

H , = I (1 Cos2<|> + 1 Cos4(J> +-•) 


If a linear magnetization characteristic is assumed for the wing 
core (permeability constant or »1) and a suitable radius chosen, we 


where D^, and D^, are effective demagnetizing factors in the y' and 
z 1 directions respectively. If H , is assumed to be fairly constant 
along any particular radial we have: 


Equation 4.4 becomes 


A g(b)b db + M , term 


For constant geometry and neglecting the M , term: 


Applying a magnetizing field and a through-wing field simultaneously: 
L <f fl (k Cos2<J> + —) + I (k Sin2<f> +-) 1 x 


Cos2<j> +-) - I (1 Sin2<f) +-) 


Neglecting terms above 2$: 

L <i (I - I ) Sin4(j) + I I Cos4<J) 
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Fig.5.13 Variation of roll 
torque with 4>, 1 ^, I m « 
Simplified linear theory. 
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This equation is illustrated in Fig. 5.13. Where <j) = O, it is 
seen from Eqn. 5.23 that: 

L « 1 1 - (5.24) 

r m 

Thus, in a fully symmetric configuration, the maximum torque for 
given total available ampere-turns will occur at: 

I m = I r = °) 

5.6 GFUN results 

5.6.1 Baseline geometry - Effect of variation in material properties, 
low and intermediate applied fields 

Where the permeability of the core is constant (linear solution, 

typical of low applied field levels) the effects of the magnetizing 

(^J ) and through-wing (J ) fields are independent. Since, under these 
m r 

conditions, the spanwise magnetization is everywhere proportional to 
J m and the through-wing field proportional to J^, the roll torque is 
expected to vary as the product J J . This is confirmed in Fig. 

5.14 - 5.16, being as predicted by Eqns. 5.9 and 5.24. 

The variation of torque with core permeability, permeability held 
constant within each solution, is of interest. Fig. 5.17 showing com¬ 
paratively low sensitivity to permeability variations at high values 
of permeability, the core then behaving as if infinitely permeable. 

Using the baseline BH curve it is clear that there exists some 
maximum level of applied field commensurate with the whole core lying 
in the initial constant permeability region of the BH curve. At 
higher applied field levels the permeability of certain strongly 
magnetized volumes of the core will progressively fall, eventually 
the bulk of the core settling onto the terminal (saturated) region of 
the BH curve. During this process the magnetization is no longer 

proportional to J , rather reaching some limiting value. It might 
m 

therefore be expected that torque becomes proportional to alone. 
However at high values of through-wing field the induced magnetizations 
may no longer be predominantly spanwise, rather turning to lie more 
nearly parallel with the direction of local (applied) field (see 
Fig. 5.9). Torque may therefore reach some limiting value, or 
continue to rise (or fall, see Fig. 5.10) with rising applied fields 
as some function of and (Fig. 5.18). 
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Fig. 5.17 Torque versus permeability. Baseline core and E/Ms. 


Linear solutions. 



t/c and taper 

GPUN solves linear (constant permeability) cases directly, 
without recourse to iterative procedures, hence relatively economi¬ 
cally. The variations of the initial (constant permeability) 
torque capability of the system with various geometrical parameters 
as defined in Fig. 5.21 are shown in Figs. 5.22 - 5.24, the solid 
symbols representing the baseline core. 

These figures require some explanation. It is clear that the 
torque for a given applied field is far from being a constant per 
unit core volume and in fact does not obey any simple relation to 
geometry (such as first moment of volume about the x axis). It 


is believed that this effect is due to the fact that the effective 
spanwise demagnetizing factors are predominantly determined by the 













Fig. 5.19 Torque versus m » Baseline core geometry and E/Ms. 



Fig. 5.20 Square-root torque versus As Fig.5.19 
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Fig. 5.22 Torque versus chord. Baseline E/Ms. Core as Fig.5.21. 
Linear solutions. 



Linear solutions, 







slenderness of the core, being lower for more slender cases. Thus 
the removal of wing volume by reduction of chord or thickness lowers 
the effective spanwise demagnetizing factor, hence increasing the 
level of spanwise magnetization of the core for any given applied 
field level. This may be partly justified by inspection of Figs. 

5.25 - 5.26, showing the increase in the peak element spanwise 
magnetization with increasing slenderness. 

Tapering the core under the criteria chosen does not affect the 
core volume but displaces volume from the most effective regions 
near the tips to less effective regions near the roots. The slender¬ 
ness of the tip portions of the core is increased and since these lie 
in the regions of strongest applied fields the peak value of element 
magnetization tends to be increased as shown in Fig. 5.27. 

5,6.3 Effect of the presence of fuselage and axial magnetizing field 

It is to be expected that the presence of an unsaturated iron 
fuselage should act to increase the mean level of spanwise magneti¬ 
zation in the wing cores since it provides an easy flux path at the 
wing root (Fig. 5.28). A saturated fuselage may or may not act 
similarly. With the fuselage geometry as defined in Fig. 5.29, the 
wing core span remains unchanged, the root now being enveloped by the 
fuselage core. The table below shows some results for this geometry, 
confirming the expectation of augmentation of torque with fuselage 
present. The effect appears weak. 


Fuselage Dimensions (cms) 

Torque (Nm) 

Absent 

109.78 

50 x 5.2 x 5.2 

110.22 

50 x 7.4 x 7.4 

113.53 

100 x 7.4 x 7.4 

114.40 


If a soft iron fuselage core is used it will generally require an 
axial magnetizing field. The effect of an axial field on the wing 
cores will tend to be to rotate the spanwise magnetization vectors in 
the plane of the wing, in the sense of sweepback or sweepforward. 
Where the core permeability is constant, the magnitudes of the span- 
wise components will be unaffected. A reduction in the spanwise 
components, hence torque, will be expected where the wing core 
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permeability is falling with rising magnetization. With the fuse¬ 
lage dimensions fixed at 50 x 7.4 x 7.4 cms (Fig. 5.29) and a near 
uniform axial field of representative strength applied by a large 
Helmholtz E/M pair of otherwise arbitrary dimensions this expectation 
is confirmed by the results shown in Fig. 5.30. 

The magnitude of the axial magnetization components induced in the 
wing core by the axial field will depend to a large extent on the 
axial slenderness of the core, which will in turn be most strongly 
affected by the core's thickness to chord ratio, relatively thick, 
narrow chord cores being least powerfully affected. Very thin cores, 
such as the F-16 case studied later, may thus be seriously affected 
although no relevant data currently exists. 

5.6.4 Effect of sweepback 

Sweepback of the wing cores as defined in Fig. 5.31 rotates the 
easy axis of wing core magnetization away from the spanwise direction, 
but does not affect the core volume nor the position of the centroid 
of volume. The induced magnetization vectors with a purely spanwise 
magnetizing field are thus expected to themselves be sweptback, at 
some angle probably less than the geometrical sweep angle. This 
would, in general, lead to lower torque per unit applied field. 

However, at low sweep angles it would appear that the increasing 
slenderness of the core, caused by the chosen geometry, dominates, 
leading to slightly augmented torques at modest sweep angles, as shown 
in Fig. 5.32. 

When an axial magnetizing field is applied, components of that field 
act along the easy axis of magnetization of the sweptback core, in¬ 
creasing or decreasing the spanwise magnetization components (hence 
torque) depending on the field polarity. For a sweep angle of 30 
degrees Figs. 5.33, 5.34 show the effect to be significant. At 
relatively low spanwise fields it is seen that powerful axial fields 
of either polarity reduce the torque for particular spanwise and 
through-wing field levels. 

5.6.5 Behaviour at high levels of roll torque generating field 

At high applied roll field levels the core becomes saturated over 
most of its volume. The induced magnetization components then behave 
as vectors of constant strength but variable direction, as predicted 
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Fig. 5.33 Effect of axial field on sweptback core. Baseline E/Ms« 
Core as Fig. 5.31. Sweep = 30°. 5 point BH curve. 



As Fig. 5.33. 



for ellipsoids in Section 5.5.1. It is not immediately obvious 
whether, under these conditions, increases in applied field should 
lead to increases in torque generated and, indeed, whether or not 
some absolute maximum torque capability will exist in this case for 
each particular core geometry. 

The studies in Section 5.5.1 suggest that an absolute maximum 
should exist and that this maximum can be realised with finite 
applied fields of appropriate sense and direction. 

Various high applied field cases have been computed with GFUN 
(Figs. 5.35, 5.36) with somewhat inconclusive results, a clear limit 
or limiting trend not being identified. It should be noted here 
that the torque integration schemes in GFUN are expected to become 
progressively less reliable at increasing applied field levels. 

This is due to the difficulty of accurately resolving on the surface 
of the external control volumes (Section 5.3, Fig. 5.5 and (43) ) 
the "model" field (due to the core's induced magnetization and now of 
essentially fixed magnitude) from the total field, which becomes 
mostly due to the applied field from the E/Ms. 

The peak torque levels achieved in Figs.5.35, 5.36 are, however, 
at least one order of magnitude in excess of apparent existing LMSBS 
requirements (5). 

5.6.6 Pseudo F-16 wing core performance 

It is understood that a representative aircraft type for use in 
LMSBS design studies has previously been chosen to be the F-16 fighter. 
This type presents a considerable challenge to the SIM roll scheme 
since its wing thickness (hence volume) is very low, the blockage 
effects of the fuselage would necessitate choice of model wingspan 
considerably below the 50% of test section width used heretofore and 
the extreme taper both in chord and thickness leaves relatively little 
volume in the magnetically most effective regions of the wing (the 
tips). 

Each wing panel has been crudely represented with GFUN as a hexa- 
hedral slab, uniformly tapered in both thickness and chord, with 
approximately the same span, total core volume and moment of volume 
about the chordwise centroids of volume as a typical F-16 model (Fig. 
5.37). The element distribution, particularly the element aspect 
ratio, within the core is at the limits of what is generally 
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acceptable with GPUN, hence the results must be regarded as subject 
to increased uncertainty (see Section 5.7) but again comfortably 
exceed apparent LMSBS requirements (5) at moderate applied field 
levels (Figs. 5.38 - 5.41). Ref. 5 was, in fact, based upon a 
requirement for 140 Nm roll torque at the scale computed herein, this 
in turn being based upon a requirement for testing at atmospheric 
stagnation pressure, Mach. 0.85. 

5.6.7 Effect of variations of E/M geometry 

Non-circular E/Ms, for instance those in Fig. 5.42, exhibit 
improved packing around the wind tunnel test section compared to the 
circular baseline E/Ms. The field distribution in the region of the 
model will be altered also, although the effect tends to be slight 
due to the relative remoteness of the E/Ms. The performance of 
different E/M configurations may be approximately normalized by an 
appropriate measure of the field in the region of the model, but 
calculated with the model absent. This has been done by computing 
the mean field level along the y axis, taken over the span of the 
model. It is seen in the tables below that the effect of E/M 
geometry on the generated torque under these conditions is small. 


E/M geometry 

Normalized J 

r,m 

Torque (y = 1000) 

Fig. 5.4 

1000.0 

120.51 

Fig. 5.42 

894.97 

119.78 


For reasons other than production of rolling moments, the classical 
array of 8 E/Ms distributed in the yz plane may not be preferred. A 
16 E/M system, shown in Fig. 5.43 has therefore been computed as an 
example with E/M performance normalized as above: 


E/M geometry 

Normalized J 

r,m 

Torque (y = 1000) 

Fig. 5.43 

869.65 

112.82 


It is clear that the performance of the SIM system is not strongly 
affected by the detail geometry of the E/Ms, hence permitting con¬ 
sideration of alternative E/M geometries and configurations with the 
existing GFUN results being approximately applicable provided E/M 
performance is normalized by the model's near field using Fig. 5.44. 
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5.7 Verification of GFUN data 
5.7.1 Alternative torque computations 

No direct alternative computations of the performance of any 
representative SIM system have been attempted since this would 
necessarily require access to an alternative computer program of 
comparable power and sophistication to GFUN, preferably solving 
the magnetostatic system by an entirely different method. Such 
programs probably do exist. 

However, as mentioned earlier, the reliability of GFUN's pre¬ 
diction of magnetization of iron regions has been considered good 
for many years and there seems no particular reason to regard the 
computed wing core magnetizations as being subject to any more 
than the usual levels of uncertainty (43). Were this the only 
source of uncertainty, the computed torques could be regarded as 
likely to be accurate, within the limits of geometrical and other 
assumptions, to plus or minus very few percentage points at all 
computed applied field levels for the baseline and similar geometries, 
with somewhat greater uncertainty in certain cases, such as the F-16 
core, where the iron element distribution was sparser than desirable 
(see Section 5.7.3). 

The main potential source of uncertainty must be regarded as the 
torque integration schemes themselves, since these have been specially 
developed for the SIM computations and hence not subjected to such 
extensive testing and verification as the rest of GFUN. 

It is possible with GFUN to arbitrarily fix the induced magneti¬ 
zation in iron elements, thus effectively converting them to permanent 
magnet material. If this is done in such a way as to approximately 
preserve the typical spanwise magnetization in the SIM cores (Fig. 
5.45) then direct and representative verification of the torque 
integration schemes is possible. Alternative calculations have been 
made using the computer program FORCE (Appendix 2) which calculates 
forces and torques by elementary numerical integration of the relevant 
vector products of applied field and core magnetization over the 
volume of the core. This method differs fundamentally from the 
methods used in GFUN. Results for the geometry of Fig. 5.45 are as 
follows: 
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Torque integration scheme 

Predicted Roll Torque 
(Nm) 

GFUN by Maxwell field stress 
integration over surfaces of 

113.94 

external control volume. 

140 x 40 x 20 cms. overall 

FORCE by vector product 


integration over volume of 
core 

116.82 


The discrepancy of approximately 2.5% is considered acceptable 
though leaving scope for future improvement. 

Certain comparisons have been made of predicted torques over 
differing GFUN control volumes with typical results as follows: 


GFUN control volume overall dimen- 

Relative torque to 

sions (cm) 

standard volume 

.140 x 40 x 20 (assuming symmetry 

1.0 

in the yz plane) 


140 x 40 x 20 (no symmetry) 

0.999993 

.180 x 70 x 50 (yz symmetry) 

0.9856 


These results are similarly considered acceptable. 

As mentioned in Section 5.6.5 it is thought that the accuracy of 
the GFUN torque integration schemes will fall with rising applied 
field level when the SIM cores are well into saturation. No direct 
high field computations were made with permanent magnet cores 
specified and this was a serious omission but has been partially 
rectified by more detailed analysis and computation using existing 
data. Specifically, since the publication of Ref.48, the induced 
magnetization distribution of two high field GFUN cases has been 
used as input data to FORCE, each GFUN element now representing a 
uniformly magnetized tetrahedra of permanent magnet material. 

The vector product integration procedures carried out by FORCE 
are not likely to yield accurate estimates of the "true" torque pro¬ 
duced in these particular cases, indeed such procedures would other¬ 
wise be used in GFUN, being relatively straightforward. The reasons 
for this cannot be fully explained herein, but are principally 
founded on the fact that GFUN’s prediction of local external field 
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is markedly superior to that of local.internal field, hence magneti¬ 
zations . In these cases the predicted internal magnetization dis¬ 
tribution will be somewhat erroneous, but the field external to the 
core should be well resolved (see (41, 42) ). Fig. 5.46 shows that 
discrepancies do exist between GFUN's torque predictions and the 
FORCE estimates from GFUN magnetization data, and that the magnitude 
of these discrepancies rises with rising applied field level. In 
the principal region of interest (J <10,000 A/cm 2 ) it is considered 
that the agreement is sufficiently close, bearing in mind the comments 
above, to manifest significant extra confidence in the GFUN torque 
integration scheme, although further verification would clearly be 
advantageous. 

5.7.2 Experimental measurement of torque with low applied fields 
and correlation with GFUN predictions 

Experimental verifications of the bulk of GFUN's predictions are 
not possible without an array of powerful (high field) E/Ms. Such 
an array was not available. However, the 8 main E/Ms from SUMSBS 
became available during the reconstruction of that system, enabling 
some representative low field (approximately constant permeability) 
torque measurements to be made. These E/Ms are not, unfortunately, 
axisymmetric but it is believed that the SIM system is not particu¬ 
larly sensitive to detail E/M geometry, rather to the mean applied 
field levels in the region of the SIM cores (see Section 5.6.7). 

The geometry and characteristics of the experimental SIM cores is 
shown in Fig. 5.47, the E/M layout in Fig. 5.48 and the experimental 
layout in Fig. 5.49. Torque could be applied to the wing cores via 
a fine thread, pulley, scale pan and weights (Fig. 5.49), the 
experimental procedure being to measure the equilibrium angular 
positions of the wing with known E/M currents and applied torque. 

There are, as predicted by Eqn. 5.23, generally two such positions 
in each octant, one stable, the other not. Fig. 5.50 shows a 
typical calibration curve for one octant, the general form of which 
is repeated for all combinations of currents. Taking the (measured) 
reference angle for zero torque at zero roll current as 91.6°, Eqn. 
5.23 may be used to generate a set of theoretical calibration curves 
with insertion of an appropriate scaling factor. This factor is 
chosen so as to zeroise the cumulative torque error over all measured 
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Fig. 5.47 Geometry and characteristics of SIM cores for experimental 
verificatio n of GFUN predictions. 
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data points, taking the measured angular positions as correct. 

Fig. 5.51 shows the comparison between these theoretical curves 
and all measured data. 

Certain irregularities in the experimental data are inevitable. 

The non-axisymmetry of the E/Ms is expected to result in slight 
loss of left-right symmetry of the calibration curves, taken about 
their nominal peaks, in Fig. 5.51. Some evidence of this is 
apparent but the effect appears weak. Friction and irregularities 
in the wing assembly support bearings, together with the usual 
uncertainties in measurements, introduce significant random and 
systematic errors, which are considered to be of comparable magnitude 
to the observed residual discrepancies in Fig. 5.51. It is thus 
concluded that this Figure effectively verifies the simplified 
theory of Section 5.5.2. 

The GFUN representation of the E/Ms was, of necessity, somewhat 
idealised (Fig. 5.52), and a small correction to E/M current density 
proved necessary to achieve correct predicted field strengths in the 
region of the model. Complete measurement of the field of the 
experimental E/M array and comparison with the GFUN representation 
proved to be too laborious, therefore the current density correction 
was calculated by measurement of the through-wing field components 
along a datum spanwise axis of the experimental SIM cores, with 
the cores absent. This axis corresponds to the reference axis of 
the E/Ms (^ = 90 , see Fig. 5.51) rather than the spanwise axis of 
the SIM cores with the cores in their zero roll torque, zero roll 
current orientation (p = 91.6 ). The comparison between experimental 
and predicted (corrected current density) fields is shown in Fig. 

5.53. Since the magnetizing field E/Ms are nominally identical 
to those creating the through-wing field, the data shown effectively 
verifies the complete E/M array, although detail discrepancies in 
the field distribution are still possible. GFUN's predicted field 
distribution (corrected) is shown more fully in Fig. 5.54. 

Comparison between selected performance curves from Fig. 5.51, 
two GFUN predictions and fitted linear theory (Section 5.5.2) is 
shown in Figs. 5.55 and 5.56. Only two GFUN points were computed 
due to the lengthening of the calculations caused by the complex E/M 
geometry. Taking into consideration the extensive computational 
idealisations and various possible sources of experimental error. 
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the agreement is considered to be excellent. The apparent dis¬ 
crepancies may, in fact, be unrepresentatively small (see Section 
5.7.3), being less than the estimated experimental error alone. 

It is concluded that major systematic errors in the predictions 
of GFUN at low or moderate field levels are unlikely. 

5,7.3 An assessment of likely accuracies of GFUN results 

Considering previous results with GFUN (40, 41, 42, 43) regarding 
prediction of magnetization levels, the agreement between the GFUN 
torque integration scheme and an alternative method applied to a 
representative low field case (Section 5.7.1.) and the performance 
of GFUN in predicting torque for the experimental system (Section 
5.7.2.), it is thought that the error magnitude in any GFUN prediction 
for the baseline or similar geometries at low applied field levels 
should not exceed 10%. Typical error may be significantly less than 
this figure, perhaps 5%. The error will tend to be mostly systematic 
in nature, available data showing low random content. This implies 
that trends in performance should be reliably identified provided 
fairly consistent formulations of the problem (for instance element 
distributions) are used. This was in fact done. The effects of 
geometrical and other idealisations are not included in the above 
figures. It being anticipated that increased inaccuracy may occur 
with rising applied field levels and a predominantly saturated core, 
it would seem appropriate to increase the figure for peak anticipated 
error to perhaps 20% for intermediate applied field levels (arbi¬ 
trarily 2000 < m < 10000A/cm 2 ). At still higher applied field 
levels the predictions become progressively less relevant to immediate 
requirements (5) since the torques predicted are high and the peak 
E/M fields required to achieve those torques are outside the limits 
of existing technology (Section 5.8.2.). The high applied field 
results (J^ ^ > 10000 A/cm 2 ) should therefore, perhaps, be regarded 
for the time being as somewhat speculative and requiring further 
verification, such as more detailed testing of the GFUN torque 
integration scheme at high field levels. 

Significant departures from the baseline core geometry (sweep, 
etc.) are achieved only by utilising iron element distributions that 
are undesirably sparse. A version of GFUN exists (using the sub¬ 
program GETM 400) that can deal with up to 400 independent iron 
elements, rather than the lOO in the standard batch program, which 
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could tackle most cases herein with nearly ideal element distri¬ 
butions. This facility did not become available during the period 
of this study and anyhow would require much greater CPU time for a 
solution, CPU time for solving the magnetizations being limited and 
varying approximately as the square of the number of iron elements. 
Thus, cases such as those involving core sweepback, the effect of 
axial fields and the F-16 cores, are subject to increased uncertainty. 

Estimates of the likely peak error magnitudes in the above cases 
are principally based on judgement but realistic estimates are 
thought to be: 


Case 

Peak antici¬ 
pated error. 

Comments 

Baseline 

low fields 

10% 


high fields 

20% 

Rising with applied field 

very high fields 
(J > 10000) 

50% 

Requires further veri¬ 
fication. Ditto. 

Sweepback < 20° 
low fields 

20% 

Standard integration volume. 
No usable symmetry. 

Sweepback > 20° 
low fields 

25% 

Increased integration 
volume necessary 

high fields (positive 
Axial) 

30% 


high fields (negative 
Axial) 

40% 

Poor convergence of 
solutions. 

Axial field with fuselage 
low fields 

20% 

Large integration volume. 

high fields 

30% 

No usable symmetry. 

F-16 

low fields 

J high fields 

40% 

60% 

Undesirably sparse element 
distribution. 


The principal idealisation inherent in all cases herein is the 
representation of wing chordwise cross sections(airfoil sections) as 
rectangles, since it has not (Section 5.4) been thought that slab 
cores buried inside non-magnetic aerodynamic envelopes would be used 
in practice. In order that the existing GFUN predictions be appli¬ 
cable to MSBS cases with true airfoil sections some appraisal of the 
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effects of this idealisation is necessary. However, it is not 
obvious on what basis, apart from cross sectional area, GPUN data 
can be matched to true sections. In the computed i’-16 cases the 
moment of volume about the centroid of area was chosen arbitrarily 
to generate a GFUN representation. The thickness so chord ratio, 
defined conventionally, could also be used, leading to generally 
similar results with classical streamline wing sect ions. Typical 
comparison between a GFUN section and a true section is shown in 
Fig. 5.5 7 : 


GFUN slab core 



Fig. 5.57. Comparison of GFUN and true core cross 
sections 

Since the spanwise slenderness is scarcely affected, it is 
expected that the mean spanwise magnetizations should be near 
identical in both cases, though the detail distributions must be 
different. The chordwise and through wing slenderness will be 
affected by the change of section, though the precise value of 
through wing slenderness is not thought to be critical at low or 
moderate fields. The (increased) chordwise slendcrness with the 
true section will lead to somewhat greater sensitivity to axial 
fields, the amount by which the effects of the axial field are 
amplified being perhaps in the range O - 50%. 

5.8 Disc ussion and conclusions 

5.8.1 _F urther use of GFUN in SIM computations 

GFUN 1 s creators could not have envisaged that it would see use 
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in this type of problem. Because of this, several detail features 
are inconvenient'for SIM computations. Some of these, such as the 
lack of a proven torque integration scheme, have been partially 
rectified (by Simkin at Rutherford Laboratory) during the course of 
this study, but some remain. Examples are the inability to exploit 
all the symmetry existing in all problems, the lack of a torque 
integration scheme with the high element number version of GFUN, and 
the doubts concerning torque predictions at high applied field levels. 
It is believed that these difficulties can be overcome relatively 
easily by further improvements or expansions of the program code. 

If this were done it seems likely that GFUN could provide torque (and 
force) predictions to a much higher level of accuracy than those 
estimated for the results herein. Definition of E/M and model core 
requirements for specific performance demands (all degrees of freedom) 
for specific LMSBSs should then be possible, to adequate accuracy for 
preliminary and intermediate LMSBS design and cost studies. It is 
difficult to imagine GFUN or similar programs becoming sufficiently 
accurate for precise predictions of performance (say better than 1% 
accuracy) in all cases of interest, due principally to the geometrical 
idealisations required in the formulations of the problem, though this 
level of accuracy should be attainable for certain simple cases. 
Empirical calibration of practical systems would therefore appear 
mandatory. Universally high accuracy may, however, be considered 
unnecessary, for instance where significant over-capacity is incor¬ 
porated in LMSBS E/Ms. 

5.8.2 Application of the SIM system to LMSBSs 

Even taking account of the relatively large errors considered 
likely to exist in the computed results, it is clear that the apparent 
torque capability of the SIM system considerably exceeds those pre¬ 
dicted heretofore for other rolling moment generating systems. For 
instance, with the baseline geometry (5 point BH curve), the value of 
600,000 Nm/m 3 of (permanent magnet) core predicted in Section 4.3, 
perhaps 300,000 Nm/m 3 of wing volume, is exceeded at approximately 
4000 A/cm 2 in all E/Ms. 

Subject to the provision of adequately powerful E/Ms the SIM scheme 
must be considered a viable contender for LMSBS application. The 
matter of E/M design requires further analysis, but data is included 
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(Fig. 5.58), showing the baseline magnetizing field at 1000 A/cm . 

The baseline through-wing field is similar but indexed by 45°. 
Examining this data it is difficult to imagine the specified E/Ms 
being operated above, perhaps, 3000 A/cm^, with existing technology 
superconductors, whence the peak field within their bores becomes 
approximately 6.6 T. That value of field may not be increased by 
geometrical adjustments, only by improvements in superconducting 
E/M technology. Whilst certain optimisations of E/M geometry 
remains possible, it is seen that a fundamental limit to the available 
torque from a particular SIM system exists. 

Where peak performance is required it is easily seen that the E/Ms 
should generally be located as close to the model as possible 
(Fig. 5.59). 

The matter of cross couplings between roll and other degrees of 
freedom due to the presence of the spanwise magnetizations, or the 
(applied) wing core magnetizing or through wing fields, has not been 
directly addressed here. It is felt that the analysis presented in 
Section 4 will remain approximately applicable though the presence of 
axial magnetizations and/or the lack of fore-and-aft symmetry in 
swept wing (such as the F-16) cases will result in relatively complex 
couplings, mostly second order in magnitude. 

This study has indicated that at least with respect to available 
roll torque, the SIM system is viable for application to LMSBSs. 
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Fig. 5.58 Field intensity around baseline magnetizing field E/Ms. 

"* 2 ’ 

J = 0, J = 1000 A/cm . Fields in tesla. Dims, in metres 



Local field fairly 
independent of E/M 
location 


Field at model 
increased with E/M 
location (B) compared 
to initial location (A) 


Fig. 5.59 Improvement in E/M effectiveness with closer packing 
to wind tunnel wall. 
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SPANWISE PERMANENT MAGNET ROLL CONTROL WITH SUMSBS 


6._1 I ntroduction 

Since the Spanwise Permanent Magnet roll control technique was 
viewed as potentially the most suitable method for use with SUMSBS 
in the foreseeable future, preliminary suspension and minimum cali¬ 
bration of a pilot model has been attempted in order that this view 
might be verified. 

6.2 Spe cial features of SUMSBS for roll control 

SPM roll control requires that SUMSBS be used in its X con¬ 
figuration, achieved by rotating the model's datum axes by 45 in 
the roll sense (Fig. 6.1). The E/M geometry thus generated does 
not directly correspond to previously discussed SPM E/M arrays 
(Section 4) due to the fore-and-aft separation of the two quadrupole 
groups of main E/Ms, but the required through-wing fields are still 
generated relatively effectively (Fig. 6.2). 

The pilot SPM model merely consisted of a standard axisymmetric 
model (Appendix 1), crudely adapted to carry vaguely representative 
permanent magnet wing cores (Fig. 6.3 and Appendix 1). All magnetic 
cores were Alnico V. The awkward appearance of this model must be 
excused but it should be noted that the wing core cross sections and 
volumes are in approximately representative proportion to the fuselag 
core for an aircraft of relatively high wing volume, such as an A-10. 
The measured torque capability may, of course, be satisfactorily 
scaled to differing wing sizes (Section 6.4). 

The roll attitude sensing system, described more fully in Appendix 
1, consisted of a low power laser directed onto a model-mounted 
mirror, the reflected beam falling onto a light spot position 
detecting device. Some difficulties were experienced with stray 
reflections of the laser beam, but the sensing system is not neces¬ 
sarily seen as representing anything more than a short term solution. 

All-digital control was employed (Appendix 1) with the roll 
stabilisation loop incorporating, purely for convenience, only one 
phase advance network (see 11, 59), instead of the two (in series) 
used in all other degrees of freedom. This corresponds to previous 
practice with SUMSBS. 

The magnetic couplings occurring due to the presence of the span- 
wise magnetizations were expected to require adjustments to the 
control algorithms to alleviate their effects, but accurate estimates 
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Fig. 6.3 Geometry of pilot SPM model. Dimensions in mm.s 
All magnetic cores Alnico V 

Polarization (J) assumed as 1.25T for FORCE calculations 
Non-magnetic envelope not shown 
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>f coupling effects were not possible for two principal reasons: 


1) The iron cores of the main E/Ms preclude accurate 
representation of these E/Ms with FORCE (Appendix 2) 
and measurement of the actual field distributions in 
SUMSBS was deemed to be too lengthy. Therefore, the 
line conductor approximation to SUMSBS' configuration 
and geometry used in Section 7 (Figs. 7.5, 7.18) was 
employed for this study, with the current levels in each 
E/M adjusted so as to match the calculated field with 

the measured field at the origin of balance axes (Fig. 6.5). 
The field distributions around this origin are not expected 
to be particularly well represented. 

2) The polarization distribution of Alnico V cores will 
not be constant and uniform, as is assumed in all FORCE 
computations. The relatively low coercivity of this 
material typically results in weakening and splaying of 

the polarization vectors, particularly around the extremities 
of the core, partly due to self demagnetizing effects and 
partly due to local demagnetization by applied fields, as 
illustrated in Figl 6.4. 



Idealized Practical 


Fig. 6.4 Idealised and practical polarization 
distributions for Alnico V cores . 

For FORCE calculations the core polarization was assumed constant 
and uniform at 1 Tesla. This figure is below the saturation polar¬ 
ization of Alnico V (-1.25T, Ref. 25) but may represent an over¬ 
estimate of the effective polarization level. The forces and 
torques generated due to various classical demands (calculated) 
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are shown in Tables 6.1 and 6.2, with the E/M geometry and current 
directions clarified in Fig. 6.5. Pure roll torque or drag force 
components are generated by classical roll or drag demands respec¬ 
tively and there exist no couplings into roll torque or drag 
force from other demands (lift etc.). It is concluded that de¬ 
coupling of the controller is unnecessary for preliminary work 
since no coupling terms are of the same order of magnitude as the 
relevant primary terms (Table 6.2), with only one coupling 
appearing as significant, that being the expected coupling into 
Yaw from a Sideforce demand. Further, generous stability margins 
are available. 

Table 6.1 Required current directions for generation of classical 
demands with X configuration (SUMSBS) 

E/Ms as Fig. 6.5. 

Current directions shown as positive (+) or negative (-) as Fig. 6.5. 


Demand 
(model axes) 

T 1 

J 2 

*3 

X 4 

*5 

*6 

J 7 

*8 

Lift 

+ 

+ 

- 

- 

- 

- 

+ 

+ 

Sideforce 

+ 

- 

- ■ 

+ 

- 

+ 

+ 

- 

Pitch 

- 

- 

+ 

+ 

: - 

- 

+ 

+ 

Yaw 

+ 

" 


+ 

+ 

- 

- 

+ 


Table 6.2 Coupling terms for SUMSBS classical applied fields 
with SPM pilot model 


Fields as Table 6.1. 
Model as Fig. 6.3. 
FORCE calculations 


Demand 
(model axes) 

P , 
z* 

F , 

y 

M' 

N' 




-5 


Lift (F z .) 

0,394 

0 

10 

0 

Sideforce (F yl ) 

0 

0.394 

0 

-0.01 

Pitch (M') 

-0.03 

o 

0.03 

o 

Yaw (N') 

0 

-0.08 

0 

0.03 

Primary terms (e.c 

j. lift force 

from lift 

: demand) 

shown underlined, 

Forces given as N 

per ampere 

(relevant 

E/Ms t in 

SUMSBS) 


Torques given as Nm per ampere (similarly). 
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L 




Currents shown 
in E/Ms this side 
only 



Fig. 6.5 Geometry of 
X E/M configuration 
for FORCE SPM 
calculations. 

I 1 " I 8 =19734 * 8 A 
“ X 10 =11809.3 A 

equivalent to 20 A cable 
current in SUMSBS 
This geometry derived 
by rotation of + 
geometry of Fig. 7.5 


Dimensions in cms, 




6.3 Suspension and calibration 

Stable suspension was achieved very easily, following trimming 
of D.C. offsets within the controller (to provide standing currents 
to support the model's weight) and adjustments to various loop 
gains. This latter requirement was due to the large increase in 
the deadweight and moments of inertia of the model compared to the 
unwinged case. Suspension quality was good. 

Time only permitted the absolute minimum calibration to be 
attempted, that being measurements of roll torque versus E/M 
currents for a model at one position and attitude (the usual datum). 
For convenience, roll torque was applied via a lightweight torque 
rod attached to the extreme nose of the model (Fig. 6.6). This 
technique has not been widely used in the past and has certain dis¬ 
advantages, perhaps principally the likelihood of applying axial 
forces to the model as well as roll torques. 



Only one E/M current could be accurately monitored with available 
equipment, but it is assumed that the current increment in all E/Ms 
due to the applied roll torque will be nearly equal. The restriction 
of the model to one fixed position and attitude (position/orientation 
error integrators operating in all control loops) further validates 
this approach. The line conductor approximation to SUMSBS (Fig.6.5) 
was again used to generate an estimate of the roll torque per ampere 
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(8 E/Ms) with the specified model. 

Results are presented in Fig. 6.7 with Figs. 6.8-6.10 
illustrating certain aspects of performance. 

6.4 Discussion 

The torque capability demonstrated is at least ten times higher 
than has been possible previously, under comparable conditions, 
with other SUMSBS roll control systems (11, 23). 

The calibration curve of Fig. 6.7 exhibits slight nonlinearities 
and a slope approximately 50% less than the computed estimate shown. 
However, the shortcomings in the experimental technique, the gross 
simplification of E/M (hence field) geometry in the computed E/M 
configuration and the partly unknown model magnetization distri¬ 
butions are thought to be sufficient to account for these discre¬ 
pancies and neither point is regarded as a matter for any concern 
at this stage. 

Scaling the measured calibration curve, using Eqn. 4.4, to a 
realistic aircraft geometry, taken to be the Cessna Citation I (a 
convenient choice with zero sweepback) indicates a torque capability 
of the order of : 

8.2 gm.cm./Ampere ; perhaps + 80 gm.cms. normally usable: 

- the model span being taken as 4*5 inches with representative core 
volumes. This is lower than the pilot SPM model due to the lower 
thickness and span, also the taper of the Citation wing. 

SPM roll control is concluded to be entirely practical for use 
with SUMSBS. 
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Fig.. 6.7 Experimental measurements of roll torque versus E/M current 


increment for pilot SPM model. 



Fig. 6.10 SPM model in suspension, model rolled from 
datum orientation in sense of starboard 
wing up. 
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7. EXTREME ATTITUDE TESTING WITH MAGNETIC SUSPENSION AND 
BALANCE SYSTEMS 
7.1 Introduction 

Contemporary small MSBSs have not previously been capable of 
suspending models over a particularly wide range of attitudes, the 
axis of the model's fuselage typically being limited to approxi¬ 
mately 30 degrees excursion from the major axis (wind axis) of the 
MSBS (MIT(B)) or less. Whilst restricted ranges of sensitivity 
of position sensing systems represents the primary limit to such 
excursions in most cases, it is undoubtedly true that none of the 
existing MSBSs (apart from the modified SUMSBS) have been configured 
with the intention of suspending models significantly outside this 
range, indeed in the majority of cases the MSBSs would be magnetically 
incapable of realistic operation over a substantially wider range. 

It would be highly desirable for a LMSBS to be capable of 
supporting and restraining typical test models over a wide range of 
test attitudes under representative test conditions. This would 
enable rapid exploration of the aerodynamic characteristics of aircraft 
or missile models over their complete attitude envelope, including of 
course the regimes where the model may be thought of as being wholly 
or partially stalled. Such regimes are currently of great interest 
in military work and tend to be rather intractable analytically. 

Current test methods in this class generally require a range Of model 
supports, each functional over a restricted range of attitudes. 

Support interference corrections can be virtually impossible since 
the supports often lie in regions of grossly turbulent wake. 

In the context of this report "extreme" attitudes are classified 
as any attitude outside the model axis excursion limit (30 deg.) 
mentioned above. It is thought that extreme attitude capability for 
a MSBS should include model axis angular excursions up to and 
including 90 degrees in at least one plane, permitting, for instance, 
model testing over the range of angle of attack of -90 to +90 degrees. 
Yaw and roll excursions need not be over such wide ranges but ideally 
would be so. Since the model's magnetic characteristics are com¬ 
pletely insensitive to the wind direction, combined pitch, yaw and 
roll excursions over a -90 to +90 degree range would encompass, with 
"reversal" of the model in the MSBS axes, the complete range of 
possible attitudes. 
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Several fundamental difficulties arise when considering model 
suspension at extreme attitudes. These include : 

1) Identification of E/M array geometries and configura¬ 

tions capable of generating, via field and field 
gradient components, forces and torques on the 
model in the required senses and magnitudes over the 
full range of model attitudes. 

2) Synthesis of control algorithms capable of accommodating 

large changes in model aerodynamic characteristics and 
magnetic couplings to the E/Ms. 

3) Design of position, attitude and other sensors to 

monitor wide ranges of model motion. 

Part of 1) above is addressed here, that is, the inclusion of 
adequate versatility into the E/M array configuration. Sizing the 
E/Ms thus specified to satisfy particular absolute force and torque 
requirements must be performed separately. 

7.2 Theoretical background 

7.2.1 Required field and field gradient components 

Forces and torques with conventional slender, axially magnetized 
models are assumed to be predominantly created by the following 
field and field gradient components (Appendix 5) : 


Table 7.1 Required field ant 
Force/torque 

F^, (axial force) 

F^, (sideforce) 

F^, (normal force) 

L' (rolling torque) 

M* (pitching torque) 

N' (yawing torque) 

(magnetizing field) 

(magnetizing field) 

Axes as Fig 


field gradient components 
Component 

H xx' 



H , (spanwise magnets) 

yz 



H , (soft iron fuselage 

o core) 

H , (soft iron spanwise 
y o magnets) 
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In "normal" suspension model and balance axes coincide and 
these components correspond, neglecting the cases of roll and 
magnetizing fields, to: 


Fig. 7.1 Model and tunnel axis 
systems 



Pitching or yawing the model through 90 degrees translates 
these components into: 


respectively. 


It may be immediately noticed that all nine primary field 
components : 
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are required independently at the origin for the full range of 
model attitudes to be usable 

The generalized couplings between field and field gradient 
components in tunnel and model axes may be calculated as follows: 


The transformation matrices for Yaw (V), Pitch (0) and Roll 
($), in tunnel axes are as follows: 


I Cosip 

-SintJj 

0 

' 

CosO 

0 

Sin0 


1 

0 

0 

Sinip 

Cosip 

0 


0 

1 

0 


0 

Cos<j> 

~Sin<j) 

0 

0 

i : l 


-Sin0 

0 

Cos0 j 


0 

Sincj) 

Cos<j) 


The couplings for a vector in model axes are given by the product 
of the reverse coupling matrices : 

Iv J IcosQcosip CosOSintJ; -SinG 1 | v I 


V , :: Sin0Sin<|>Cosip 
y -Cos<j>Sinip 


Cos<{>Cos<jj Cos0Sin<j) 

+Sin<J>Sin0Sinip 


V 


y 


(7.1) 


V . 

z' 


I 


Sin<J>Sinlp 

|+Cos<J>Sin0Cosi|> 


Cos<J>Sin0Sinip 

-Sin<|>Cosip 


Which may be written as: 


V' = 


AV 


The couplings for field components are given straightforwardly by: 


H' = 


AH 


Whereas for field gradient components : 


V' = AV and using (7.2) 

V'.H* = AV.AH = BVH 
where the reduced form of B is given by : 


H , 1 

XX ' 


ai 2 

2aja2 

2aia 3 

a 2 2 

2a 2 a 3 

a 3 2 

H , 

xy' 


aibj 

aib2 + 
a 2bi 

aib 3 + 
a 3 bi 

a2b2 

a 2 b 3 + 

a 3 b 2 

a 3 b 3 

H , 

xz’ 


a i c l 

ajC2 + 
a2Ci 

aic 3 + 
a 3 c j 

a 2 c 2 

a 2 c 3 + 

a 3 c 2 

a 3 c 3 

H . 

yy' 


bi 2 

2bjb 2 

2b,b 3 

b 2 2 

2b 2 b 3 

b 3 2 

H , 1 

yz' 


bjci 

biC2+ 

b 2 ci 

b,c 3 + 

b 3 c. 

b 2 c 2 

b 2 c 3 + 
b 3 c 2 

b 3 c 3 



ci 2 

2cic 2 

2c, c 3 

C2 2 

2c 2 c 3 

c 3 2 
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The matrix A being written as : 


ai a 2 a 3 
bi b 2 b3 


jci c 2 c 3 J 

It is seen that field components behave as vectors during axis 
rotations whereas field gradient components do not and that the 
couplings from applied field gradients to field gradients in model 
axes are extremely complex. Certain simplifications are possible 
where rotation is restricted to one plane, say pitch, where : 


^ » <)> 

= 0 


[», 


Cos0 0 

-Sin0 


[»» 





H , 

y 

= 

0 1 

0 


H 

y 

- (7. 

2 a) 





Sin0 0 

. 

COS0 


»>l 



and: 












Cos 2 0 


0 

-2Cos0Sin0 

0 

0 

Sin 2 0 | 

H 

XX 

H , 
xy' 


! 0 

Cos0 

0 

0 -Sin0 


0 

H 

xy 

H xz’ 


Cos 0 Sin 0 

0 

Cos 2 0-Sin 2 0 

0 

0 

-Sin0Cos0 

H 

xz 

H , 

yy> 


0 


0 

0 

1 

0 


0 

H 

yy 

H , 

yz' 


0 

Sin0 

0 

0 Cos0 


0 

H 

yz 

H . 

zz' 


Sin 2 0 


0 

2Sin0Cos0 

0 

0 

Cos 2 0 j 

H 

zz 






- 

(7. 

3b) 




Model forces and torques can be given by any of the above using s 

F = v /M.VH dV and T = x H + r. x(M.Vh) dv 

- following Eqns. 3.1 

The force and torque couplings will differ from 7.2, 7.3 etc., due 
to the effects of the spatial variations of the applied fields over 
the volume of the model's core. Where the applied fields are 
relatively uniform over the core these effects will be of second 
order and the force and torque couplings will be approximated by 
7.2, 7.3 etc. 

It is considerably more convenient in many cases to represent the 
couplings by transforming model magnetizations into tunnel axes 
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rather than field and field gradient components into model axes : 



model magnetizations being given by : 
M = 4> 0 ¥ M' 


1.2.2 E/M configurations for multiple independent field and 
field gradient component generation . 

The requirement to generate 9 field or field gradient components 
independently necessitates at least 9 independent E/Ms. The desire 
for symmetry in the E/M array (Section 1.2.1) acts to increase this 
figure. 

A symmetric quadruplet of E/Ms as shown in Fig. 7.2 can generate 


four field gradient components 


: the originjH (Ij, I 4 =-12* I 3 )/ 


tti, I 2 , I 3 = U), H (I lf I 2 , 13 = I 4 ) r and 


H (Ii, I 2 , 13 = I 4 ),but it is immediately seen that H , H 
zz o XX o yY 0 

and H are not independent. Two field components at the origin 


H (I a , I 2 - -I 3 , I 4 ) and H (l lf I 3 = -I 2 , l 4 ) f may also be 


generated. Modified geometries (Fig 7.3), in fact corresponding 
to the "vertical" or "lateral" E/Ms in SUMSBS, are found to be 


relatively weak : 


shown. If H (or H ) were regarded 


as a prime component of field for this sub-configuration and the 

'stray'components H and H (H ) were countered by some other 
xx yy zz 

means, the modified quadruplet would be a useful generator of four 
independent field or field gradient components f H , H , H , 


A symmetric pair of E/Ms can generate one field and one field 
gradient component independently at the origin H , and H as 


shown in Fig. 7.4, with stray components H and H 
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Fig. 7 0 2 Symmetric quadruplet 
of E/Ms 



Fig. 7.3 Modified 
quadruplet of E/Ms 



Ct 


Fig. 7.4 Symmetric 
pair of E/Ms 


n 



Conventional E/M configurations can be considered as an assembly 
of quadruplets and pairs as defined above and it is straightforward 
to arrange for all required field and field gradient components to 
be somehow provided. 

In fact, current MSBS configurations (eg. MIT(B), SUMSBS) are 

frequently capable of generating 8 field or field gradient components 

(H , H „ H , H , H , H , H , H ), effectively though not 
x y z xx xy xz yy zz 
o o o o o o 1 o o 

completely independently and typically lack capability only in H yz . 
It may thus be expected that conventional configurations are 
capable of sensible operation over a wider range of attitudes than 
has been exploited previously. 

7.3 Identification of maximum force and torque capabilities of 
a MSBS 

Any extreme attitude capable MSBS will incorporate a number of 
independent E/Ms, perhaps between 10 and 20. Each E/M tends to 
create all nine independent field and field gradient components at 
general points. Thus, with a suspended model at a general attitude 
and position, each E/M creates force and torque components along or 
about each axis (x, y, z). It is clear that there will seldom exist 
a unique solution for E/M currents with any particular force and 
torque requirement. The absolute maximum force and torque capa¬ 
bility of a particular system may not, therefore, be directly 
analytic. 

With 6 independent E/Ms however, there is no serious difficulty 
since, neglecting demagnetization or saturation effects (see Section 
7.7), the required senses for force and torque specify variables 
ki-ks such that : 

klF y = k 2 F z = k 3 L = k 4 M = k 5 N = F x _ (7<5) 

The effects of the 6 E/Ms may be written : 


1 F 

X 


l a i 

a 2 

a 3 

a 4 

a 5 

a 6 


l h 

F 

y 


b 1 

b 2 

b 3 

b 4 

b 5 

b 6 



F 


Cl 

c 2 

c 3 

c 4 

C 5 

C 6 


I 3 

‘‘I 

| 

dl 
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d 3 

d 4 

*5 

* 6 . 


14 

M 


e i 

e 2 

®3 

e 4 

e s 

e e I 


I 5 

N / 


i f > 

f 2 

f 3 

f4 

fs 

‘•I 




where aj-f 6 are dependent on the systems geometry and magnetic 
characteristics and are functions of model position and attitude; 

Il~ Ie representing the current levels in each E/M. It is immediately 
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seen that where ai~f 6 are mostly non-zero, I 2 -l 6 (say) are dependent 
variables with only Ii unspecified. Since all E/M currents are now 
fixed relative to each other, the absolute maximum capability is 
found by increasing all currents until one reaches its limiting 
value, that value being determined by the E/M design. 

Useful results may be obtained for the case of 7 E/Ms, starting 
with the equations for a 6 E/M system operating with one current 
(say Ij) at its limiting value (Ii P ), as above : 

I F 1 jai a 2 a 3 a 4 a 5 a 6 Ii P 

I F y bi b 2 b 3 b 4 b 5 b 6 I 2 

F C1 C 2 C 3 C 4 C5 C 6 I 3 

L di d 2 d 3 d 4 ds dg I 4 

M ei e 2 e 3 e 4 es e$ Is 

N | fl f2 f 3 f * f *l I 1 * 

1 2 —1 6 being determined using Eqns. 7.7 and 7.5. Introducing 
another E/M (I7), the equations may be reformed : 

F ai a 2 a 3 a 4 as a 6 a 7 1 I 2 

F y bi b 2 b 3 b 4 bs b 6 b 7 13 

F z -I^p ci _ c 2 c 3 c 4 C 5 c 6 C 7 I 4 

L di d 2 d 3 d 4 d 5 d 6 d 2 I 5 

M ei e 2 e 3 e 4 e 5 e 6 e 7 I 6 

n/ jf,| |f 2 f 3 £, £5 f 6 £ 71 |l7 

Since the coefficients ai-fi and Ii P are constants for any 
particular model attitude and position, resolution of 7.8 follows 
a similar procedure to 7.6. This leads to a second current 
reaching its limiting value, say I 2 P . 

If the number of E/Ms is represented by n, then the cases above 
yield the criteria that where n=6 or n=7, n-5 E/Ms should be 
operated at their peak design currents to achieve a particular 
maximum force and torque capability. It is now argued that this 
criteria holds for all n>6. It is not, however, obvious which 
n-5 E/Ms should be chosen to be at their respective peak currents, 
but this can be resolved by computing all permutations. 

Cases where the matrix in 7.6 etc. is sparse may be dealt with 
by reduced forms of the analysis above. 




154 



7.4 "+" E/M configuration study 

7.4.1 Geometry 

The revised configuration for SUMSBS can fall into this category 
(see Sections 1.2.1 and 2.2.1). The performance of this system 
cannot be explicitly computed by FORCE (Appendix 2) since the main 
E/Ms are iron cored, however a system having similar proportions 
but arbitrciry exact dimensions has been computed to yield an 
indication of the usable attitude range of SUMSBS (Fig. 7.5). 

This representative system is chosen to be at ten times the scale 
of SUMSBS for convenience. Single loop E/Ms are used to simplify 
calculations but are expected to provide adequate approximations of 
the performance of SUMSBS' E/Ms, the current levels being adjusted 
where appropriate to approximate the relative field capabilities of 
the E/Ms in SUMSBS. 

7.4.2 Pe r f ormance 

As might be expected fromEqns. 7.3 etc., the generalized couplings 
between E/Ms and model core are likely to be extremely complex. 

Since the inherent limitations of the optical position sensing system 
used in SUMSBS would, at least initially, restrict large angular 
displacements to one plane, specifically the xz (pitch) plane, all 
computations presented here are made with the model restricted to 
this plane, with its centroid fixed at the origin of balance axes. 

Thus x=y=z=ij>=cf>=0. 

With or without the above restriction, the magnetic couplings 
between representative E/Ms (Fig. 7.5) and a typical model (Fig. 7.6) 
are easily found for given model attitudes and positions using FORCE. 

It is not, however, clear in what form the couplings should be 
presented; as couplings from individual E/Ms (Fig. 7.7a-d) or as 
couplings from conventional groupings of E/Ms (Fig. 7.8a-j). In 

this case, the latter approach fully separates the creation of 
magnetic forces and torques but is less general, partly because the 
groups of E/Ms are not mutually exclusive. 

The effects of the spatial variations of the applied fields around 
the origin of axes may be found by recomputing Figs. 7.8 with model 
polarization concentrated into a single dipole at the origin (Figs. 
7.9a-j). It is found that these results are identical to direct 
computation from Eqns. 7.3 etc. with the field capabilities of the 
E/M groupings (calculated) as: 
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Fi.g„ 7.7 C ouplings from + configuration individual E/Ms 

E/Ms as Fig. 7.5 
Model as Fig. 7.6 


Data is shown specifically for only 4 E/Ms. Remaining E/Ms may be 
deduced from this data by symmetry, resulting in the sign table below: 

Signs of force/torque components in: 

Fig.7.7a Fig. 7.7b 



Fig.7.7c Fig.7.7d 
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Fig. 7.7c (ii) 

+ configuration 
E/M 2 (4,6,8) 

Torques 

TORQUE 

(Nm) 
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Fig. 7.8 Couplings from + configuration classical applied fields 


E/Ms as Fig. 7.5 
Model as Fig. 7.6 


Fields are defined by the relationship of E/M currents as shown below: 


Field 

II 

12 

1 3 

14 

15 

16 

1 7 

18 

19 

IlO 

Lift 

+ 

N/A 

- 

N/A 

- 

N/A 

+ 

N/A 

N/A 

N/A 

Drag 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

+ 

+ 

Auxiliary 
drag No.l 

- 

N/A 

" 

N/A 

- ■■■ 

N/A 

“ 

N/A 

N/A 

N/A 

Auxiliary 
drag No.2 

N/A 


N/A 

“ 

N/A 


N/A 


N/A 

N/A 

Sideforce 

N/A 

- 

N/A 

+ 

N/A 

+ 

N/A 


N/A 

N/A 

Pitch 


N/A 

+ 

N/A 

— 

N/A 

+ 

N/A 

N/A 

N/A 

Auxiliary 
pitch No.l 


N/A 


N/A 

+ 

N/A 

+ 

N/A 

N/A 

N/A 

Auxiliary 
pitch No.2 

N/A 

“ 

N/A 

>. 

N/A 

+ 

N/A 

+ 

N/A 

N/A 

Axial 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


+ 

Yaw 

N/A 


N/A 

+ 

N/A 


N/A 

+ 

N/A 

N/A 


- following Fig. 7.5 


This representation of applied fields is derived from the fields 
classically applied with the model at its datum position (lift, pitch 
etc.), augmented so as to encompass all possible combinations of 
currents. 

The separation of Auxiliary Drag into two constituents (Figs.7.8c 
& d) is done since simultaneous application of these two fields with 
one in the reverse sense produces the through-wing field used in SUMSBS 
for SPM roll control. 





Fig. 7.8c 

Auxiliary drag 
No.l field 
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Fig. 7.8d 
Auxiliary drag 
No.2 field 
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Fig. 7.9d 
+ configuration 
Auxiliary drag 
No.2 field 
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+ configuration 
Pitch field 
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Fig. 7.9g 

+ configuration 
Auxiliary pitch 
No.1 field 
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Fig. 7.9h 

+ configuration 
Auxiliary pitch 
No,.2 field 


TORQUE 


(Nm) 
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Fig. 7.9j 

+ configuration 
Yaw field 
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Table 7.2 Field capabilities for + configuration 


Field components (T) / 10 S AT /E/M 


E/M grouping 

B 

XX 

B 

xy o 

B 

xz G 

B 

YYo 

B 

yz 0 

B 

zz 0 

B 

x,y,z 0 

Lift 

O 

0 

0.1969 

o 

o 

O 

0 

Drag 

0.0676 

0 

0 

-0.0338 

0 

-0.0338 

0 

Auxiliary 

0.002 

0 

O 

0.1252 

0 

-0.1272 

0 

drag 1 
Auxiliary 

0.002 

0 

0 

-0.1272 

o 

0.1252 

0 

drag 2 
Sideforce 

0 

0.1969 

0 

0 

0 

0 

0 


B 

X o 

B 

y o 

B 

z o 

All gradients 



Pitch 

0 

0 

-0.1267 

0 




Auxiliary 

-0.1147 

0 

0 

0 




pitch 1 
Auxiliary 

-0.1147 

0 

0 

0 




pitch 2 
Axial 

-0.0808 

0 

0 

0 




Yaw 

0 

0.1267 

0 

0 





The correlation between the couplings for the representative model 
and the zero volume model, the latter being essentially theoretical 
results, is generally fair, the cases where the E/Ms are relatively 
remote from the model (Figs. 7.9b, 7.9h) exhibiting the best agreement. 
The "theoretical" curves may thus be adequate for synthesis of experi¬ 
mental control algorithms, provided generous stability margins are 
available and such algorithms would not then require updating 
following minor model changes. Peak performance of the system is, 
however, unlikely to be approached unless the computed couplings of 
Fig. 7.8 or similar are used directly. 

The variations of the maximum force and torque capability can be 
found from the data in Fig 7.7 using the methods of Section 7.3. 

Again there is difficulty due to the range of possible combinations 
of model forces and torques. However magnetic torque capability, 
apart from roll torque, is seldom a critical factor in MSBSs and 
with model excursions restricted to the pitch plane, sideforce 
capability may be largely neglected. Figs. 7.lO(a-b) thus illustrate 
the maximum force capability for the representative + system with the 
force vector constrained to lie in the plane of pitching and with all 
torques zero. The sense of the force vector is clarified in Fig.7.11. 
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The direction of the magnetic force vector required to oppose the 
net aerodynamic force would normally be expected to lie in the quadrant 
defined by 0$3<90°. Transient forces may be outside this range, as 
may be the model's weight vector. However, for the purpose of 
identifying the absolute maximum force capabilities of a MSBS, and 
since the curves of Fig. 7.10 are repeated for 3*‘3±180° the 0^3<90° 
quadrant is seen to be of principal interest. Some of the content 
of Fig. 7.10 may now be simplified by choosing three components of 
aerodynamic significance, namely pure lift force, pure drag force and 
normal force, all defined in Fig. 7.12, and adding the envelopes 
within which all curves for 0^3<90 0 must lie (Fig. 7.13). The axial 
E/M current limits (E/Ms 9 and 10) have been factored by 0.5984 to 
approximate the measured relative performance of the main and axial 
E/Ms in SUMSBS. Computation for the zero volume model is illustrated 
in Fig. 7.14 and exhibits significant differences in certain areas, 
such as drag forces at low angles of attack. Effects of variations 
in the E/M array geometry are well illustrated with this form of 
presentation. For instance, doubling the current limit applied to 
the axial E/Ms (E/Ms 9 and 10) increases, as would be expected, the 
drag force capability at small angles of attack but also attenuates 
the otherwise rapid fall off in normal force capability with angle of 
attack increasing from small angles (Fig. 7.15). 

The discontinuities in slope of the curves in Figs. 7.13 - 7.15 
are caused by changes in the distribution of E/M current limiting, for 
example at some particular angle of attack E/Ms 1 and 7 (say) may be 
operating below their current limits, with all other E/Ms at their 
respective limits (see Section 7.3). Small increases in angle of 
attack may result in the currents in E/Ms 1 and 7 increasing, other 
currents remaining fixed, and at some angle all E/Ms will be at their 
limiting currents. Further increases in angle of attack will 
generally require E/Ms other than E/Ms 1 and 7 to be reduced in 
current, overall performance then being delineated effectively by 
a different performance curve. This may be illustrated for one case, 
chosen to be the normal force curve from Fig. 7.13, Fig. 7.17 clearly 
showing tne changing E/M current distributions. 

7.4.3 Discussion 

Reserving treatment of the significance of the absolute magnitudes 
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Pig. 7.14 Maximum performance of the + configuration with zero 



ANGLE OF ATTACK CDEG3 


Fig. 7.15 Maximum performance of the + configuration with E/Ms 
9,10 doubled in strength. Otherwise as Fig.7.13 
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of forces and moments in Figs. 7.7 - 7.10 and 7.13 - 7.15 for 
Section 7.7, the computed couplings do nevertheless indicate the 
practicality of the particular + configuration geometry studied 
for high angle of attack testing. 

Figs. 7.10 and 7.13 clearly illustrate that the principal aero¬ 
dynamic forces (appearing in the plane of pitching) can be 
countered relatively effectively at most angles of attack with 
only the normal and drag force minima around 30° angle of attack 
presenting any serious difficultyi Augmentation of the main 
drag field partially alleviates this problem, the severity of which 
is heavily dependent on the precise geometry (see (49)). Without 
indulging in further analysis it is concluded from Figs. 7.8f-j 
that sufficient capability for the purposes of generating pitch and 
yaw torques (roll not being treated here) exists at all angles of 
attack, leaving only one fundamental flaw in the E/M configuration, 
that is the inability to generate sideforce at 90° angle of attack 
(Fig. 7.8e). This is due to the lack of capability in H , which 
was shown in Section 7.2.1 to be required for this purpose. 

The particular geometry studied appears to be practical for 
operation up to at least 60° angle of attack, then being limited 
principally by falling sideforce capability. It is thought that 
this view remains broadly valid for all + configurations of con¬ 
ventional proportions (49), though detail optimisation of perfor¬ 
mance is clearly possible. One example of this might be the in¬ 
corporation of "skew" into the lateral E/Ms (Fig. 7.16), although 
this technique does partially destroy the symmetry of the 
configuration. 

It is important to recall at this stage that the field component 
required for generation of roll torque using the SPM or SIM schemes 
is H y Z > the E/M array were augmented in such a way so as to 

include capability for this component, then operation to 90° angle 
of attack and, in fact, through 360° pitch attitude range Would 
appear possible. 

The somewhat bizarre fluctuations in E/M currents apparently 
required to achieve maximum force capability with rising angles of 
attack (Fig. 7.17) are almost certainly impractical due to the low 
E/M current slewing rates available in realistic MSBSs. However, 
the complete reversals of current occurring with small increases in 
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angle of attack at certain points in Fig. 7.17 generally indicate 
that the E/M in question is relatively ineffective in generating 
force in the required sense at that particular attitude. Approxi¬ 
mation to the absolute maximum performance with appropriate 
restrictions on current fluctuations are thus expected to be possible 
with relatively small performance degradation, though no analysis 
has yet been attempted. 

7.5 "X" E/M configuration study 

7.5.1 Geometry 

The lack of capability in H of the + configuration can be 

yz 

rectified by simply rotating the model's axis system (in the roll 
sense) so as to convert the + configuration to an X configuration. 

A representation of this version of SUMSBS has been computed, using 
identical E/M and model geometry and characteristics as the represen¬ 
tative + system of Section 7.4 (Fig. 7.18). 

7.5.2 Performance 

Following identical analytic procedures to those employed for the 
+ system, Figs. 7.19a-b show the main E/M couplings, those from the 
axial E/Ms (E/Ms 9 and 10) remaining as shown in Fig. 7.7d. Figs. 
7.20a-g show the couplings from classical E/M groupings, with 
couplings to the zero volume model presented in Figs. 7.21a-g. 
Immediately it can be seen that sideforce capability exists at 90° 
angle of attack (Figs. 7.20/7.21d). 

Computing the maximum force capabilities in the plane of pitching 
yields a disappointing result (Fig. 7.22). Although the lift force 
capability at zero angle of attack is augmented, compared to the + 
configuration, as would be expected since all 8 main E/Ms now contri¬ 
bute to lift force, the performance around 45° angle of attack is 
disastrous. The zero volume model exhibits similar results 
(Fig. 7.23), whilst doubling the axial E/M strength effects a 
significant but probably inadequate improvement (Fig. 7.24). 

Eqns. 7.3 etc. may be used to illustrate the reason for the 
poor performance around 45°, using the (calculated) field 
capabilities of the E/M groupings : 
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Fig. 7.19 Couplings from X configuration individual E/Ms 

E/Ms as Fig. 7.18 
Model as Fig. 7.6 

Data is shown specifically for only 2 E/Ms. E/Ms 9 and 10 remain 
as shown in Fig. 7.7d. Remaining E/Ms may be deduced from the data 
shown by symmetry, resulting in the sign table below: 

Signs of force/torque components in Fig. 7.19a 

E/M F F F L M N 

___ _x_ y z 

1 - + + + -+ (shown) 

7 + -* - + 

8 -- + + + + 


Signs of force/torque components in Fig. 7.19b 
E/M F x F y F z L M N 

3 +----- (shown) 

4 + + - + - + 

5 +-- + + + 

6 + + -- + - 
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Fig. 7.19a(ii) 

X configuration 
E/M 1 (2,7*8) 
Torques 
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Fig. 7.19b(ii) 

X configuration 
E/M 3 (4,5,6) 
Torques 



angle; of attack cpegs 
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Fig. 7.21 Couplings from X configuration classical applied fields 
- showing zero volume model 

These Figures are directly comparable with Figs. 7.20, which 
should be consulted for details of model, E/Ms and current relation 
ships. 

SOLID lines represent the zero volume model. 

BROKEN lines show the standard model and are identical to the 
corresponding Figs.7.20. 


Fig. 7.21a 
X configuration 
Lift field 
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Fig. 7.21b 

X configuration 
Auxiliary drag field 
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Fig. 7.2 1c 
X configuration 
Sideforce field 
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X configuration 
Auxiliary sideforce 
field 
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Fig. 7.21e 
X configuration 
Pitch field 
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Fig. 7 - 21£ 

X configuration 
Auxiliary pitch field 


TORQUE 


<Nm) 
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Fig. 7.21g 
X configuration 
Yaw field 
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Fi g. 7.24 Maximum performance of the X configuration. Axial E/Me 
doubled in strength. Otherwise as Fig.7.22 
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Table 7.3 Field capabilities for X configuration 
Field components (T) / 1Q 6 AT / E/M 


E/M grouping 

B 

B 

B 

B 

B 

B 

B 


xx o 

xy o 

xz o 

yy o 

yz o 

zz o 

x,y,z 

Lift 

0 

0 

0.2786 0 

o 

0 

0 

Drag 

0.0676 

o 

O 

-0.0338 

o 

-0.0338 

0 

Auxiliary drag 

0.0039 

0 

0 

-0.0019 

0 

-0.0019 

0 

Sideforce 

O 

0.2786 

0 

0 

0 

0 

0 

Auxiliary 

0 

0 

0 

0 -0 

.2525 

0 

0 

sideforce 









B 

x o 

B 

y Q 

B 

z o 

All gradients 



Pitch 

0 

0 -0.1793 

1 0 




Auxiliary 

-0.2296 

O 

O 

0 




pitch 








Axial 

-0.0808 

O 

O 

o 




Yaw 

0 

0.1793 

0 

0 





Eqn. 7.4 gives, for the x and z components of force in this case: 


and since M^=M Cos© and M z =-MSin9 we have, at 45 angle of attacks 

F = k (H - H ) 

x xx xz 

F = k (H - H ) - < 7 ’ 9 > 


These equations are heavily coupled, directly via the term 

and via the lack of independent capability in or H^, shown in 

the table above. The requirement is generally to develop positive 

F^ and positive F^ simultaneously in some proportion and it is seen 

that since H contributes to F and F in opposite senses and that 
xz x z 

the (coupled) and components are relatively small in the 

chosen configuration the available magnitudes of forces will tend 
to be low. 

7.5.3. Discussion 

Referring to Tables 7.2 and 7.3 it is seen that where H , H 

xx yy 

and components only are developed at the origin, the sum of 
these components is always zero. This corresponds to Maxwell's 
equations for a conservative field with no free poles : 
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V . B B + B + B = 0 - (7.10) 

xx yy 7 , 7 . 

applied in free space. By symmetry it can be seen, using Eqn. 

7.10, that in the case of symmetric axial E/Ms or the complete 
symmetric X array (E/Ms 9 and 10 or X E/Ms 1-8) : 

B ’= B -hB - (7.11) 

yy 7.7 xx 

De< oupling ol Eqn. 7.9 using 7.11 is thus not attractive. In 
the case of the vertical or lateral E/M quadruplets from the t array 
(+E/Ms 1-10) however : 

B ' -B » B - (7.12) 

yy zz xx 

Modifications to the geometry of the X array, such as realignment 
of the E/Ms as shown below, may prove to affect Eqn. 7.11 sufficiently 
powerfully to partially restore the viability of the X configuration 
but no investigations have yet been attempted. 



195 




7.6 E/M 


[figurations for LMSBSs 


The examples of + and X configurations studied herein cannot be 
assumed to be fully representative of the general classes of such 
configurations, since large variations in overall proportions and 
relative E/M capabilities are possible. However, the inability to 
generate H , hence sideforce at 90° angle of attack with the + 
configuration, and the difficulties around 45° angle of attack with 
the symmetric X configuration are thought to be fundamental. 

There is no objection in principle to combining the two configura¬ 
tions in an attempt to combine their advantages and negate their 


disadvantages, indeed this was envisaged (but never carried out) in 
Ref. 49. Numerous detail variations are possible, such as : 



(Fore and aft spacing of E/Ms is not shown 
and may itself be subject to considerable variation) 
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Configurations a) and b) above correspond closely to proposed 
SIM roll arrays (Section 5.8). Indeed, although roll control has 
generally not been treated in this Section it is the case that the 
X array as specified can effectively generate the required through 

wing fields (H ,) for all model angles of attack. Large roll dis- 
yz o 

placements from datum would not be permissible with the X array, 

however, since H , cannot be generated at <J>=45° (Section 4.3.1). 
yz o 

The hybridised arrays above (a-c) do not suffer this limitation. 

The employment of model core geometries other than the conventional 
slender, axially magnetized configuration is possible, particularly 
with permanent magnet core material, but would generally lead to other 
complications, including increases in the complexity and/or severity 
of cross couplings between model degrees of freedom. 

7 .7 Discussion 

The absolute magnitudes of forces (or torques) shown in the various 
Figures in this chapter are of somewhat limited significance. They 
do not represent the maximum attainable forces for the chosen con¬ 
figuration, merely the forces attainable within the constraints of 
given E/M current limits and model characteristics. The absolute 
maximum attainable forces with a permanent magnet core are principally 
set by the magnetic behaviour of the core (demagnetization and usable 
polarization) but these limits could not be identified at the time of 
writing due to the lack of data concerning high coercivity permanent 
magnet materials (see Section 4.3). The corresponding limits for a 
soft iron core will be set by the behaviour of the induced magneti¬ 
zations as the core material enters its saturation region, presenting 
a challenging computational problem (see Section 5). 

The studies made herein are principally intended to indicate the 
general practicality of the + and X configurations, partly by 
identification of any magnetic shortcomings in these configurations, 
and the magnitude of the control difficulties caused by the changing 
couplings between model and E/Ms. Further investigation of the latter 
area appears appropriate but cannot be undertaken here. 

In the design of an extreme attitude capable MSBS it is probably 
more appropriate to consider the E/M array simply as a generator of 
field and field gradient components in the test section, rather than 
an assembly of "lift", "drag", etc., E/Ms, as has been common 
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previously. Usable model attitudes are unlimited if all components 
can be generated effectively and independently. It would be to 
advantage if the applied fields were relatively pure and uniform 
in the region of the model in order to minimise cross coupling 
effects and this requirement may perhaps be met by the use of 
Helmholtz pairs and corresponding optimised quadruplet geometries. 

7.8 Experimental demonstration of suspension with SUMSBS at high 
angles of attack 

7.8.1 Introduction 

Section 7.4 concluded that SUMSBS should be usable, in its + 
configuration, to around 60 angle of attack with conventional models. 
Preliminary attempts at suspending models up to this angle of attack 
have been made and are reported briefly here. 

7.8.2 Special features of SUMSBS for high angle of attack operation 

The position sensing system, fully described in Appendix 1, 

represents only an expedient and very short term solution to the 
difficult problem of accurately monitoring model position and 
attitude over wide ranges of attitude. The system devised achieves 
high angle of attack capability only by mechanical rotation and 
realignment of the four main sensors, axial position being monitored 
by a single fixed sensor (Appendix 1, Fig. 7.26). 

Without modification, SUMSBS was capable of suspending conventional 
axisymmetric models to around 30° angle of attack. 

It was decided to attempt initial high angle of attack suspension 
at around 50°, requiring suitable relocation of position sensors of 
course, and decoupling of the now grossly coupled and quite inoperable 
controller. Decoupling proceeded using the computed couplings of 
Section 7.4, and was performed in model axes, since the position 
sensing system operates essentially in this axis system. The 
precise methods used are described below, but it should be carefully 
noted that they represent a considerable simplification from the 
ideal, made entirely for convenience, it being believed that 
successful stable suspension would be sufficient to verify the 
fundamental theory of this Section, even if achieved on a less than 
wholly representative basis. 
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Fig. 7„8j shows that yaw torque (model axes of course) is only 
weakly affected by pitch rotation, with Fig. 7.8e furthur indicating 
that sideforce is only affected by falling gain, easily restored by 
appropriate increase in controller loop gain. Consulting Fig. 7.7a 
it is seen that if E/Ms 1 and 7 are operated in conjunction, but with 
opposite signs of currents, then only pitch torque is generated and 
relatively effectively. The full pitch torque capability previously 
enjoyed is restored by supplementing the effect of E/Ms 1 and 7 
by the use of E/Ms 2,4,6 and 8 via the Auxiliary Pitch No.2 field, as 
Fig. 7.8h. 

The residual difficulty lies in the generation of forces in the 
required senses and magnitudes in the xz plane, Now Fig. 7.7d 
indicates that the use of E/Ms 9 and 10 (the axial E/Ms) in conjunc¬ 
tion produces a net force approximately perpendicular to the x' axis 
of the model, with the model at 50° angle of attack (Fig. 7.27). 

Forces approximately parallel to the model's x' axis are conveniently 
provided by E/Ms 3 and 5, with these E/Ms operated in conjunction, as 
Fig. 7.7b and Fig. 7.27. It was found that the currents required in 
E/Ms 9 and 10 required to oppose the relevant component of the model's 
weight proved to be very close to the available current limits 
(20 A). A standing current was thus demanded in E/Ms 2,4,6 and 8 
corresponding to the Auxiliary Drag No.2 field (Fig.7.8d), effectively 
reducing the residual model weight required to be supported otherwise. 

Careful adjustment of controller loop gains was necessary before 
stable suspension could be achieved, but the "50°" controller 
subsequently proved capable of suspension over the range of angle of 
attack of 40° to 60°. No calibration or detail examination or 
optimisation of performance was attempted but the Figs. 7.28a-g 
illustrate the full range of angle of attack explored. 

Suspension quality at the higher angles of attack was poor, perhaps 
not surprisingly in view of the highly simplified experimental 
approach, and considerable further study is thought appropriate, some 
aspects of which are discussed in Section 9.2.2. 

However Figs. 7.28 are thought to validate the fundamental and key 
conclusion of this Section, that is that the use of MSBSs for high 
angle of attack or extreme attitude wind tunnel testing is absolutely 
viable. 
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Fig. 7.28a 
-5° angle of attack 


Normal controller and position 
sensors (configuration A) 

See Fig. 7.26 for details of 
position sensor configur- 
Plates ations 

overleaf 


Fig. 7.28b 

0 angle of attack, as 7.28a 

- - 

i 


Fig. 7.28d 

30° angle of attack 

Normal controller 

Position sensors in B configur¬ 
ation. 


Fig. 7.28f 
o 

50 angle of attack 
"50 " controller 
Position sensors in C configur¬ 
ation. 


Fig. 7.28c 

15° angle of attack 

Normal controller 

Position sensors in B configur¬ 
ation. 


Fig. 7.28e 

40° angle of attack 

"50°" controller 

Position sensors in C configur¬ 
ation. 


Fig. 7.28g 
o 

60 angle of attack 
"50°" controller 
Position sensors in C configur¬ 
ation. 





RELIABILITY OF LARGE MAGNETIC SUSPENSION AND BALANCE SYSTEMS 


8.1 Introduction, the need for reliability 

MSBSs are inherently complex systems, incorporating a wide 
variety of sub-systems such as E/Ms, perhaps with associated 
cryogenic plant, power supplies, at least one set of model position 
sensors, analogue or digital closed-loop control systems and, of 
course, a wind tunnel and its ancillaries. Due largely to this 
complexity, all existing MSBSs are rather "unreliable" in the sense 
that control over model position is lost relatively frequently, 
perhaps after every few hours of operation. Model "flyaway" 
inevitably follows, since MSBSs of the types considered to date are 
open-loop unstable in at least one model degree of freedom (Earnshaw's 
Theorem). Flyaway is seldom more than inconvenient with small 
systems since models tend to be small, relatively simple and 
inexpensive and operated with low tunnel dynamic pressures. 

In the case of a LMSBS, however, high dynamic pressures may be 
required to satisfy demands for Reynolds number and Mach number, 
and winged models are likely to be operated in high lift configura¬ 
tions at these dynamic pressures, leading to large net forces. 

Despite the fact that models will be relatively massive, since it is 
essential to occupy a high proportion of model volume with magnetic 
material, the typical maximum aerodynamic force will exceed the 
model's; weight by a large factor (approximately lO in Ref. 5). 

Model accelerations following a worst-case control loss may thus be 
high and it is clear that enormous impact loads would be generated 
if a breakaway model struck a tunnel wall. Whilst it would be 
mandatory to ensure that the tunnel itself could withstand such 
impact, extensive damage to the model would seem inevitable. Con¬ 
ventional large wind tunnel models represent considerable financial 
and temporal investment and it is thus certain that model control 
loss in a LMSBS cannot be permitted to be anything other than a very 
rare event. 

The prime reason for building a LMSBS is to acquire aerodynamic 
data concerning the test model. Precisely the same goal motivates 
the construction of most large wind tunnels of any type. Test 
programs in existing large tunnels tend to be interrupted relatively 
frequently by equipment malfunctions affecting data acquisition. 

This situation is tolerated since aerodynamic test programs can 
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usually accommodate delays or slippages without serious consequence 
and testing is frequently modified and new tests devised on the 
basis of early results. There seems no particular reason to assume 
that testing in a LMSBS should be any different in this respect. 

It would seem appropriate, therefore, to conceive of the "control" 
segments of a LMSBS being highly reliable, with data acquisition and 
ancillary functions remaining relatively unreliable. 

At this juncture the usage of the word "reliable" must be clarified. 
"Failures" Of general systems are often classified into two types, 
these being "unreliability" where hardware mechanically fractures Or 
electrically shorts etc; or "lack of integrity" where no identifiable 
physical failure occurs but nevertheless the system fails to perform 
its design task due to some other factor, perhaps adverse environ¬ 
mental conditions, such as electrical interference. Often the use 
of these two terms overlaps considerably and precise classification 
of particular failures is not always possible. The uncertainty over 
the terms seems especially severe in the case of MSBSs where the 
causes of loss of model control frequently do not involve simple 
mechanical or electrical failure. Herein, therefore, the reader 
should be aware of the use of the terms "reliability" and "integrity" 
somewhat interchangeably, the former usually encompassing the latter, 
though not vice versa, and implying the probabilistic likelihood of 
continuous control of the model being maintained within predetermined 
operational constraints. 

The inherent reliability of key hardware in contemporary MSBSs is 
not especially high and the large number of subsystems required leads 
to low overall reliability unless special measures are taken. 
Development of a LMSBS will involve high technical risk since at 
least the first LMSBS will be unique in overall architecture and 
most of its subsystems can be expected to incorporate significant 
technical innovations and exhibit considerable differences from any 
contemporary hardware in design, duty and environment. It is thus 
expected that considerable design effort is both essential and worth¬ 
while, and that the total system configuration must be especially 
formulated, in order to ensure that a LMSBS be capable of reliable 
operation. 

This Section seeks to explore the impact on overall MSBS archi¬ 
tecture of a requirement for highly reliable operation. Detailed 
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definition of hardware features or performance is not attempted. 

8.2 Reliability engineering in LMSBSs 
8.2.1 Modular design in LMSBSs 

The simplified architecture of a MSBS shown below suggests a 
natural division of the overall system into five major subsystems. 
There may be certain overlap in function between these subsystems, 
for instance, the system controller may include processing of model 
sensor outputs, but the division is valid insofar as the interfaces 
between subsystems are principally unidirectional and either 
fundamentally fixed or entirely optional. 



Module 

Primary interface to 
next subsystem 

Comments 

Model 

Fixed: Dynamics 

- 

Model and 
system sensors 

Optional: probably 
electronic 

May be sensitive 
to environment 

System con¬ 
troller 

Optional: probably 
electronic 

May include sensor 
processing 

Power supplies 

Fixed: Current, Voltage 

Closely coupled to 
•B/Ms 

E/Ms 

Fixed: Magnetic fields 

Mutual coupling 
between E/Ms 


Fig. 8.1 Inherent modularity of MSBSs 
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A highly modular design approach is thought to be particularly 
applicable in LMSBS development where the complete system will be 
novel, highly complex and composed of devices and systems from widely 
differing branches of engineering. Modularisation can be advan¬ 
tageously applied at the subsystem as well as the system level and 
the potential benefits of this approach include:- 

(1) Simplified system development. Devices or subsystems 
from several different working groups may be brought together 
in an agreed and considered overall architecture with minimum 
difficulty. 

(2) Easy upgrading of subsystems. A major advantage in 
LMSBS design since very few will be built and hence each will 
tend to be rather experimental in nature. As operational 
experience is gained and design or construction techniques 
improve, enhancements to overall performance may be made 
most economically by subsystem upgrading. 

(3) Simplified development and maintenance of subsystems. 

LMSBSs will certainly include many similar components within 
many modules or subsystems (such as several power supplies). 
Considerable benefits accrue if groups of similar components 
can be arranged to be fundamentally identical and the number 
of different groups of components reduced. 

(4) Reduced overall cost. If (1) and (3) can be success¬ 
fully implemented the overall manufacturing costs will be 
reduced due to elimination of the duplication of component 
or subsystem development. 

8.2.2. Operational task classification and general reliability 
considerations 

It has already been noted that the tasks undertaken by a MSBS may 
be divided into two major classes, that is "Suspension" and "Balance". 
These tasks are of quite different criticality to the system and 
merit different approaches to their reliability. 

The suspension task is critical and failure to prevent model 
flyaway win be regarded as catastrophic failure of the system. No 
particular constraint need be placed on the quality of model suspen¬ 
sion for fulfilment of this task. 

The balance task is of lesser criticality and may in principle 
be unconditionally aborted at any time and with any frequency 
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without triggering catastrophic failure. The quality of sus¬ 
pension should be high but useful fulfilment of this task merely 
requires that the overall system's mean time between any event 
causing unacceptable degradation of suspension quality be reasonably 
long, with a correspondingly short mean downtime. These criteria 
should be satisfactorily met by a well-designed system without 
recourse to special measures. 

The suspension task requires that the mean time between flyaway 
be very long indeed, perhaps in the order of thousands of hours. 
Taking into account the fact that each LMSBS will be of complex and 
largely novel design there seems no possibility of realising this 
objective without recourse to special design features. Further, if 
component failure rates in a LMSBS follow the classical "bath-tub" 
curve then comparatively large numbers of failures would be 
expected in the early life of that system. Unless massive and 
lengthy test programs, such as those carried out for manned space 
missions, are contemplated, it would be unreasonable to expect that 
all latent design or manufacturing errors or "infant mortality" 
failures could be detected. This assertion leads directly to 
consideration of a design strategy where the system may somehow 
survive at least one random single point failure - a "fault 
tolerant" design strategy. 

Fault tolerant systems may be conveniently subdivided into two 
classes: 

(1) Fail operational 

(2) Fail soft 

The fail operational strategy implies that sufficient spare 
capacity is somehow incorporated so that the system may continue 
operation following a failure, usually at undiminished performance. 
Hardware redundancy is essential if this strategy is to be successful 
in surviving any single point failure. The fail soft strategy 
abandons the requirement for continued operation of the system at 
full performance but seeks to ensure that the consequences of any 


Classical 
bathtub 
curve for 
component 


"Infant 

^mortality" 



failure rates 


Rate 


Constant fail.rate 


"Wear- 

out" 



failure are sufficiently attenuated so as to permit planned and 
graceful degradation of system performance. The requirements 
placed on LMSBSs suggest something of a hybridisation of these 
two philosophies. Critical suspension modules must fail opera¬ 
tional although certain degradation of performance may be per¬ 
missible. Balance hardware need not fail operational but must 
fail soft insofar as disruption of the suspension task must be 
prevented. Both fail operational and fail soft strategies 
generally require that propagation of failures through the system 
("cascade" failures) be prevented. Graceful degradation of 
overall system performance following a failure is quite acceptable, 
perhaps desirable, and could be manifested as an operational 
strategy where the suspension task is rapidly aborted, such as 
by initiation of a rapid shut-down sequence (RSS), on detection 
of any failure in key hardware. It is interesting to note that 
this option is not open to, for instance, the designer of an 
advanced flight control system. In order to achieve good mean time 
between RSSs either the agglomerated mean time between failure of 
all key components and systems must be long or sufficient redundancy 
must be incorporated to permit continued operation at low risk 
following one or more component failures. 

Good serviceability requires that the fraction: 

Mean downtime following RSSs 

Mean time between RSSs 

be small. It is frequently proposed in other contexts that this 
criterion be met in complex systems by holding spare components or 
subsystems in reserve, rapidly to replace failed units. Since it 
has been anticipated that a substantial proportion of the causes of 
RSS initiation might not involve classical hardware failure, 
imposition of requirements for extensive spares holdings is inappro¬ 
priate here. 

It happens that even following loss of model control there will 
exist an identifiable optimum strategy, involving suspension E/Ms 
and pcr^^ps wind tunnel controls, which would "guide" the resultant 
flyaway, so as to cause minimum damage to the model and/or tunnel. 

One type of failure more serious than that described above as 
catastrophic may be contemplated, that is, a sequence of events 
triggering cascade mechanical failures of the suspension E/Ms. 
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This type of failure is technically credible, such as through helium 
supply disruption to an array of certain types of superconducting 
E/Ms, but can and should be rendered virtually inconceivable by 
careful design of relevant systems. Such a gross failure could 
otherwise be justifiably regarded as cataclysmic. 

8.3 Configuration of reliable LMSBSs by incorporation of 
hardware redundancy and some other techniques 

8.3.1 Incorporation of hardware redundancy into the E/Ms 
and some further considerations 

Introduction of some form of redundancy into the E/M array is 
necessary if sufficient force/moment capability to maintain the 
model under control is to remain following an E/M failure. However, 
overheads are incurred on several counts, including : 

(a) Provision of surplus redundant ampere-turns. An E/M 

could incorporate redundancy by being composed of several 
independent windings, each winding creating substantially the 
same useful field, configured such that the E/M is capable of 
supplying its maximum design field with one winding failed. 

The redundant E/M must incorporate more ampere-turns that the 
equivalent non-redundant E/M and the factor : 

Total ampere-turns in redundant E/M 

Ampere-turns in equivalent non-redundant E/M 
represents an E/M cost factor which depends on the number of 
independent windings per E/M. 

(b) Increased complexity in the design of the E/M array. 

This is particularly important since provision of multiple 
separate dewars (superconducting E/Ms) may seriously compromise 
close packing of E/Ms around the working section. Some segmen¬ 
tation of dewars is likely to be required to prevent certain 
serious failures, such as LHe supply loss, cascading through more 
than one winding. Failure to pack E/Ms closely causes larger 
E/Ms to be required since the operational effectiveness of each 
is depleted. 

(c) Lower thermal efficiency of superconducting E/Ms. A 
major path for heat leakage into the LHe environment is along the 
current supply leads to each winding where they enter the inner 
dewar. Further, since the leads cannot be superconducting, 
resistive power is dissipated along their length. Proliferation 
of supply leads will be expected to increase power losses on both 



of the above counts. 

The effect of (a) is to decrease the E/M cost factor with 
increasing number of windings per E/M whereas (b) increases the 
cost factor with increasing number of windings per E/M. It is 
thus seen that there may be an optimum number of windings per 
E/M (Fig. 8.2). (c) demands that the number be low. 

It is thought that the E/M cost factor will be relatively small 
with the optimum number of windings per E/M and that that optimum 
number will be small due principally to three factors : 

(1) Use of Automatic Power Off (APO). Without APO one worst 
case failure is a power supply runaway. In this event effectively 
two windings of a particular E/M are lost after one failure since 
the equivalent of one winding must simply be used to neutralise the 
runaway winding. The reduction in the idealised cost factor with 
APO varies with the number of windings per E/M as shown in Fig. 8.3. 
The system must be designed to cope with the current transient in 
the chosen APO procedure . 

(2) "Passive" Redundancy (PR). For reasons other than 
provision of redundancy it seems certain that a major architectural 
feature of any LMSBS will be the positioning of E/Ms all around the 
tunnel working section with extensive symmetry. Generally, at 
least the mirror image of any E/M suffering failure may take all or 
part of the load previously taken by that E/M. If this fact is 
exploited the complexity of each E/M may be reduced (Fig. 8.4). 

(3) State of the art of E/M dewar design. Continuous progress 
in this area is being made by many teams engaged in the design and 
construction of cryogenic systems for a wide variety of applications. 
It appears that the overall thickness of dewar and thermal insu¬ 
lation required to achieve satisfactorily low heat leakage into the 
LHe environment can be quite small, of the order of 1.5 inches 

(50, 51, 52). Provision of separate dewars on this basis would 
not seriously compromise the packing of E/Ms with the largest con¬ 
templated LMSBSs (perhaps 8 feet test sections). Smaller systems 
would be more seriously affected, the thickness of insulation 
required not being a function of physical scale. However, full 
insulation is not necessarily required between the closest packed 
E/Ms (Fig. 8.5). 

In order to raise the fundamental reliability of each E/M winding 
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If, say, winding A fails, then 
windings B-P inclusive may all 
contribute usefully to creating 
the, say, lift force previously 
created by winding A. 


Fig. 8.4 The concept of passive redundancy 


Vacuum insulation 



Fig. 8.5 Suggested insulation requirements for close-packed E/M 


windings 



it may be decided to incorporate redundancy into the associated 
ancillary systems. These are particularly numerous in the case 
of superconducting E/Ms, typically including cryogenic plant, 
external dump resistors for discharging a winding following failure, 
and instrumentation such as helium level detectors. 

It is relatively straightforward to incorporate redundant plant, 
sensors, etc., into the E/Ms but this would only be beneficial where 
the probability of.an E/M quench being initiated by means external 
to the E/M (perhaps mechanical shock) was low compared to the 
probability due to the ancillaries themselves. 

8.3.2. Incorporation of hardware redundancy into the E/M power 
supplies 

Each winding of each E/M requires some form of power supply. 

The majority, perhaps all, of the E/Ms in a LMSBS would require 
current to be controlled with large amplitude frequency response 
to a few Hz. Terminal voltages of either sign are essential and 
it is likely that various factors would necessitate bipolar 
current capability, forcing choice of a four quadrant supply. 

Precise current and voltage requirements would depend critically 
on a system's detail design and specification but there seems no 
possibility at present of operation of a LMSBS without individual 
power supply capabilities of the order of several hundred amps and 
volts respectively. No supply is currently known to exist which 
simultaneously approaches the above requirements. 

The difficulties associated with switching powerful electrical 
supplies are severe, particularly where large quantities of 
energy (millions of Joules in Ref. 5) can be stored in each E/M. 
Large current and voltage transients might be expected during 
switching processes in this case. Therefore, whereas it is techni¬ 
cally feasible to incorporate redundancy in power supplies by 
straightforward duplication of supplies to each winding, with appro¬ 
priate means to switch to backup units following primary supply 
failure, it may be thought an unattractive possibility on practical 
grounds. The financial penalty associated with duplication of 
power supplies is, of course, a significant further consideration. 

If each E/M winding is operated with a single non-redundant power 
supply, failure of that supply necessitates immediate withdrawal 
of the winding from suspension duties since the winding current is 
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no longer controllable. With superconducting E/Ms, it is not 
necessary to power down a winding following supply failure, since 
the winding's time constant is effectively infinite. However, as 
previously mentioned, the possibility of supply runaway suggests 
immediate application of an APO procedure. The conventional means 
of achieving rapid power down with existing superconducting E/Ms 
is simply to short circuit the E/M terminals with an external dump 
resistor of suitable value. The hardware necessary for this type 
of procedure is easily duplicated, where required, to achieve 
highly reliable operation. 

It may be possible to exploit the inductive coupling that exists 
between windings, particularly between windings of the same E/M, 
which will, following Lens's Law, tend to maintain the net magnetic 
field of any E/M as constant. Controlled shutdown of one E/M 
winding might thus tend to induce a natural increase of current in 
the other windings of that E/M, with little external intervention. 
This, of course, broadly corresponds to the required system response 
following winding failure. 

Interruptions of the power utility are not out of the question 
and some buffering would seem to be essential. conventional means 
of buffering supply interruptions include : 

Emergency generators. 

Motor-generator sets in main power feeds. 

Batteries with static inverters. 

Considering the high transient powers involved, the use of any of 
the above techniques would place a considerable cost overhead on a 
LMSBS. However, since the power dissipation in the load E/Ms will 
be nominally zero, the power supplies may be regarded as operating 
by supplying energy to the E/Ms where current increase is required 
and retrieving it where decrease is required. It then becomes 
highly attractive from the point of view of energy efficiency to 
arrange for a quantity of energy to be somehow stored, perhaps 
within each supply, at least equal to the magnetostatic energy 
stored in the associated E/M windings at peak design currents. 

Energy is then simply transferred between the energy store and E/Ms 
with the utility demand being devoted merely to losses and 
ancillaries. As previously discussed, the supplies used with the 
UVa superconducting E/M MSBS featured energy storage (on internal 
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capacitors) in precisely the fashion described above (7), as do 
the supplies now used with SUMSBS (Appendix 1). It is immediately 
seen that supplies of this type could be made capable of executing 
a controlled RSS following loss of utility power, with only 
ancillaries such as control systems requiring supply buffering. 

8.3.3 Reliability and hardware redundancy in the system 
controller 

The system controller is seen as the key module of any advanced 
LMSBS, since it must encompass a wide variety of functions, which 
together directly influence the configuration and operation of all 
other parts of the total system. These functions include : 

a) Model stabilisation and control. The principal traditional 
feature of MSBS control systems. LMSBSs are certain to be operated 
with test models covering a wide range of variations in aerodynamic 
and/or magnetic configurations and characteristics. The control 
algorithms must, therefore, exhibit considerable versatility and 
robustness if they are to operate successfully, even over short time 
scales, since the aerodynamic characteristics of the test models 
will not be known in detail or to any accuracy before a test program 
commences. 

b) Prevention of model demagnetization. Permanent magnet cores 
and, to a lesser extent, magnetically soft cores, would need to be 
protected from exposure to combinations of applied fields that would 
result in partial (loss of calibration) or extensive (loss of model 
control) demagnetization of these cores. Demagnetization of 
permanent magnet cores will be possible in all realistic LMSBS 
designs and loss of magnetization of magnetically soft cores is 
credible in all MSBSs. 

c) System monitoring. The controller must incorporate 
extensive subsystem health monitoring if it is to respond intelli¬ 
gently to significant failures, such as by initiation of a RSS. 

It is obvious that the monitoring processes themselves should be 
carried out reliably. 

Very highly reliable system control hardware, implementing 
very dependable algorithms is seen as the foundation of a successful 
LMSBS design. 

Since it seems certain that the major part of the control system 
of a LMSBS will be realised with digital computers, there would be 
no particular difficulty in the incorporation of hardware redundancy. 
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Several options are available to achieve this, including : 

(1) Overall redundancy (Fig. 8.6) 

A small number, perhaps three, mainframe or mini¬ 
computers could relatively easily be configured as parallel 
independent channels such that the complete control (suspension) 
task is implemented simultaneously on each channel. Failure 
of one channel could presumably be detected by some form of 
majority voting procedure, with the failed channel being 
switched out by command signals from serviceable channels. 

This approach permits relatively straightforward execution 
in hardware, though at the cost of a considerable surplus of 
computing power (200% in a triplicated system). 

(2) Localised parallel redundancy (Fig. 8.7) 

Series connection of several stages configured as (1) above, 
produces a configuration exhibiting improved tolerance of 
multiple failures but at the expense of increased hardware com¬ 
plexity and no reduction in the required surplus of computing 
power. 

(3) Multiprocessing (Fig. 8.8) 

A relatively large number of versatile processing units can 
be connected by multiple independent data buses as shown, to 
effectively form a single processing system. This system can 
exhibit high tolerance to multiple failures where the task 
previously performed by failed units is rapidly re-allocated, 
either by software or hardware actions, amongst surviving units. 
This architectural approach appears to be favoured currently 
for ultra-reliable flight control systems, and is consequently 
under intensive development for this application (53, 54). 

The hardware and software complexity can be extreme however, 
although typically only a modest surplus of computing power 
(perhaps 50%) is required. 

The MSBS control task is broadly comparable to that of a compre¬ 
hensive flight control system in respect of the nature, complexity 
and frequency of the calculations required. However, the possi¬ 
bility of rapidly aborting the control task, such as by use of some 
form of RSS (Section 8.2.2), following detection of subsystem 
failure, limits the number of independent failures that need be 
tolerated. The development of a specialised multi-processing 
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Fig. 8.6 Overall redundancy in the model control system (schematic) 
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Fig. 8.7 Localised parallel redundancy in the model control system 
(schematic) 
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Multiple parallel independent 
internal data buses 



Multiple versatile processors (B), communicate' ' 
with multiple internal data buses via bus protection 
units (A). Failed processors switched out by command 
signals from surviving units (C). Task reallocated 
automatically amongst serviceable units. 


Fig. 8.8 Multiprocessing in the model control system (schematic) 
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controller specifically for this application may not, therefore, 
be justified. 

Further, whilst certain series segmentation of the control 
task appears possible, for instance : 

(a) Model state estimation, including sensor 

processing. 

(b) Control response calculation. 

(c) Translation to required power supply 

demands. 

the extent to which this segmentation is possible may be too 
limited to favour implementation of (2) above. 

The fact that the control computers are unlikely to represent 
the dominant cost in a LMSBS and since frequent upgrades or 
revisions of software are likely, particularly in the early stages 
of LMSBS development, appears to mitigate towards (1) above. 

Considerable attention will need to be paid to the software used 
in the control computers in order to ensure its own reliability, 
particularly since parallel operation of identical software is 
envisaged (parallel redundancy). This is, however, a very difficult 
area, the science of software testing and reliability being very much 
in its infancy, and cannot be further addressed herein. 

8.3.4 Hardware and other redundancy in the model position sensors 

All model position (also velocity etc.) sensing systems so far 
considered for use with LMSBSs suffer a serious lack of integrity, 
even where the hardware used is reliable. Specifically : 

1) Optical systems. The type of system used to date with 
SUMSBS ( (11) and Appendix 1) is inoperable where rapid variations 
of the transmittivity of the optical paths are encountered. Such 
variations must be expected to occur relatively frequently in a 
cryogenic wind tunnel, for instance with free stream condensation, 
and in all tunnels when flow visualization is undertaken. 

Methods which measure in either an analogue (e.g. T.V. or laser 
scanning) or in a pseudo-digital fashion (solid state arrays) the 
absolute position of either a target affixed to, or of a suitable 
part of, the model, may be designed to tolerate the conditions 
mentioned above (24), but all optical systems must fail in conditions 
of severe degradation of the optical path. 

2) Electromagnetic position sensor (EPS). Any high intensity 
burst of electromagnetic radiation with suitable frequency content 
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will "white-out" an EPS. Such a burst might occur, for instance, 
during electrical storms or catastrophic failure of E/M windings 
or power supplies. At present, considering the problems experienced 
by the known users of EPS systems in filtering E/M power supply noise 
(10, 55) it seems unreasonable to assume that the problem of white- 
out may be easily overcome. 

It is seen that a powerful argument exists, quite separately from 
consideration of the reliability of position sensing hardware, 
favouring the simultaneous use of at least two fundamentally 
different position sensing techniques. As an example, this may 
perhaps be achieved by use of some form of optical system together 
with an EPS, since there is no particular reason to expect simul¬ 
taneous obscuration of the optical paths across the test section and 
electromagnetic disturbance. Without such provision a LMSBS would 
be unprotected against the common-mode loss of all position sensing 
systems arising, from the adverse environmental conditions mentioned 
above. 

Subsystem redundancy may still be attractive in cases where hard¬ 
ware reliability tends to be low (eg. light sources) but the inclusion 
of totally redundant sensing systems of similar type may not be 
fruitful in the context of the previous discussion, a particular case 
in point being the EPS. 

The only method by which sensible overall redundancy could be 
incorporated into the EPS is by use of systems operating at different 
frequencies since, for instance, duplicated sets of sensing coils 
would otherwise pick up the same (possibly erroneous) signals. It 
is not, however, clear whether wide spacing of frequencies is tech¬ 
nically feasible, the available frequency band being limited at low 
frequencies by the inherent electromagnetic noise spectra of the MSBS 
and by EPS coil and other system inductances at high frequencies. 

Now, the EPS hardware physically located within the E/M array of a 
MSBS is limited, merely an assembly of excitation and sensing coils. 
These can be made very reliable indeed. The driving, sensing and 
signal processing electronics will typically be located remotely to 
the MSBS and may relatively easily be duplicated, perhaps with fast 
switching out of failed units. There may, therefore, be no require¬ 
ment for overall redundancy within the EPS. 

It would seem reasonable to consider the application of special 
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analytic techniques to position sensor outputs in order that 
spurious information be rejected. The considerable redundancy 
of information that will exist, even with limited sensor redun¬ 
dancy, recommends the techniques of analytic redundancy, that is 
the detection of false sensor outputs by comparison of all 
available data with predictions derived from a real-time mathe¬ 
matical model. Such techniques are under intensive development 
for other purposes (56, 57). 

8.4 The impact of redundancy requirements on the sizing of the 
E/M array for a representative + configuration 

8.4.1 Theoretical background 

If each E/M is divided into a number of windings, say three, 
with the intention of tolerating the failure of one winding, then 
the E/M does not necessarily require 50% (in this case, assuming 
APO) surplus ampere-turns. The effect of passive redundancy 
(Section 8.3.1) is to reduce the requirement for surplus E/M 
capacity by distributing the load previously taken by the failed 
winding over a number of E/Ms, including, but not exclusively, the 
remaining windings of the partially failed E/M. 

It happens that the methods derived in Section 7.3 may be used 
to examine the effects of E/M failures on the maximum force and 
torque capability of specific MSBSs simply by adjusting (downwards) 
the current limit applied to a specific E/M, in order to simulate 
failure of a winding within that E/M. These methods seek to find 
the absolute optimum distributions of currents in all available (in 
this case, remaining) E/Ms so as to manifest the maximum capability 
for a particular combination of forces and torques. Here, this 
will tend to yield estimates of the minimum conceivable effect of the 
failure in question. Shortcomings in the re-distribution and 
manipulation of E/M currents will result in further performance 
degradation. 

Since there exists an infinite number of combinations of model 
attitude, position and force and torque requirements it is clearly 
impractical to attempt to present a comprehensive study of the 
effects of E/M failure, even within a single specific geometry. 

The approach taken here, therefore, is to examine the effects of E/M 
failure on some of the data from Section 7, specifically the maximum 
force capabilities of the + configuration, with the model simply 
pitching, the resultant force vector constrained to lie in the plane 
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of pitching, and with zero magnetic torques. To recap, this 
system is of essentially similar overall configuration as SUMSBS, 
albeit geometrically and physically simplified to enable FORCE to 
be used to calculate the system’s performance, and at ten times 
the true scale of SUMSBS for convenience. The absolute magnitudes 
of forces emerging from this analysis are not relevant here, but 
the relative effects of various failures, shown later, will be 
closely representative of the effects of the corresponding failures 
in SUMSBS. These studies are of limited significance insofar as 
they may convey somewhat optimistic or pessimistic assessments of 
the overall influences of the relevant failures, sideforce and 
torque generation not being studied, but serve to illustrate the 
powerful effect of passive redundancy. 

Certain simplifications have proved to be necessary in order to 
prevent the computation of each set of performance curves becoming 
excessively laborious. In fact the envelopes within which 
relevant overall performance curves in Section 7 must lie (see 
Figs. 7.13 - 7.15 and accompanying text) are omitted. 

8.4.2 Results 

The configuration here is identical to that used in Section 7, 
and is specified in Figs. 7.5 and 7.6. Fig. 8.9 defines a 
numbering sequence for the E/Ms. With the sense of the three 
force vectors of interest defined in Fig. 8.10, Fig. 8.11 shows 
the system performance with all E/Ms fully operational. This 
Figure is identical to the relevant parts of Fig. 7.13. Figs. 

8.12 to 8.19 contrast this performance to that with one-third or 
two-thirds of each E/M failed. One-third failure corresponds to 
a failure, with APO, of single windings of the main E/Ms of SUMSBS; 
two-thirds to failure, with APO, of two windings in the same E/M, 
or to worst-case failure of single windings without APO. SUMSBS' 
axial E/Ms have only dual windings, but one-third and two-thirds 
failures are computed here for consistency. 

8.5 Discussion 

The content of Figs. 8.12 - 8.19 has great significance to the 
cost of a LMSBS. With one-third of an arbitrary E/M failed, the 
reduction in force capability for the cases computed is nowhere 
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Fig. 8.9 Numbering sequence for E/Ms of + geometry, as Figs. 7.5, 


The magnetic forces 



Fig. 8.10 Sense of force vectors for Figs. 8.11 - 8.19 




Fig. 8.11 Performance of representative + system with all E/Ms 
operational 

Sense of forces as Fig.8.10. Geometry as Figs. 7.5,7.6. 

This Figure shows the MAXIMUM capability for the forces as 
with all other forces and torques zero. 


shown 












greater than one-sixth and is typically around half that value. 
Should this data prove to be representative of all cases, it would 
follow that the surplus E/M capacity required to permit tolerance 
of single arbitrary winding failures with triple E/M windings is 
not greater than 20%. This should be contrasted to the value of 
50% shown in Fig. 8.3, where passive redundancy was not considered. 
Dependent upon which parts of the performance curves, such as Fig. 
8.11, prove to be critical, the requirement for surplus capacity 
may be even less than stated above. 

It is currently thought that the further cost overheads incurred 
due to the extra complexity of, say, triple windings, as compared 
to single windings of similar total ampere-turns, will prove to be 
modest and acceptable. Further, it is seen that dual windings may 
be viable, incurring a basic cost factor increment of not more than 
33%, and that large numbers of windings per E/M are unlikely to be 
worthwhile (see Fig. 8.20). 

Whilst no corroborative computed data is available, it is easily 
seen that increases in the number of individual E/Ms in a MSBS will 
tend to increase the effect of passive redundancy, thus decreasing 
the cost factors for given numbers of windings per E/M, since there 
would then exist greater opportunity for redistribution of the load 
of a failed winding. The number of E/Ms in SUMSBS is 10 and is 
thought to represent a minimum conceivable number for LMSBS designs. 
LMSBS cost factors are thus likely to be lower than those quoted 
above for comparable cases. 

Briefly to summarise this Section, it has been shown that the 
increase in cost, complexity and sophistication of LMSBS hardware 
caused by a requirement for highly reliable operation need not be 
outrageous. E/M cost factors can be kept satisfactorily low with 
apparently modest hardware complexity, appropriate use of existing 
types of position sensor should provide adequate redundancy and 
independence of information, comparatively straightforward control 
system hardware may be used and the overall power consumption of the 
complete system need not be materially affected. 
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Fig. 8.20 The effect of passive redundancy (PR) and APO on E/M 
cost factors, PR factor as Figs. 8.12-8.19 



9„ DISCUSSION AND CONCLUSIONS 


Since the content of this thesis divides into three principal 
subject areas (roll control, extreme attitude testing and reliability) 
the first part of this Section is divided likewise. These sub¬ 
sections will commence with a brief summary of material from the main 
text, followed by specific discussion and conclusions. 

9„1 Roll control 
9 „1.1 Summary 

A new method of roll control has been proposed, based on the 
magnetization, in the spanwise sense and with port-starboard symmetry, 
of the wings of winged aircraft models. The system exhibits two 
principal variations, since the spanwise magnetizations may be 
permanent (SPMs) or induced (SIMs). 

The roll torque generated by SPMs in a suitable uniform applied 
field gradient is shown to be proportional, at low applied field 
levels, to the first moment of core volume about the model's plane 
of symmetry, also, of course, to the applied field strength and the 
intensity of the spanwise magnetizations. The maximum torque capa¬ 
bility of SPMs is of major interest and is expected to depend princi¬ 
pally on the magnetic properties of the core material (almost 
certainly ReCo) but, due to the lack of available data, has proved 
impossible to identify at this stage. However a crude estimate 
indicates that substantial torques may be available. 

The existence of cross couplings between model degrees of freedom 
is recognized but it is shown that with unswept wings there exists 
only one primary coupling and relatively few other couplings of 
significance. 

A pilot model with SPM wings has been suspended in full 6-component 
control (with an all-digital control system, in fact) with no 
difficulty. The demonstrated roll torque capability and the achieved 
levels of stiffness and damping of rolling motion were adequate to 
recommend the system as a first choice for future use in SUMSBS. 

Subject to the provision of adequate E/M power, the SIM system will 
provide greater roll torque capability than SPMs, due partly to the 
higher levels of spanwise magnetization that are attainable. No full 
demonstrations of the system were possible (that is, with 6-component 
model suspension) due to the unsuitability of the available E/M 
configuration. However, comprehensive numerical computations (GFUN) 
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clearly indicate that high torques will be available, restricted 
principally by current E/M technology. Limited experimental 
torque measurements at low applied field levels exhibit excellent 
agreement with corresponding calculations. 

9.1.2 Discussion and conclusions 

It is thought that development of a powerful method of active 
roll control is absolutely essential before construction of a LMSBS 
could be attempted. 

The variant of the Spanwise Magnet scheme that could be fully 
demonstrated (SPMs) may be practical for LMSBS application, subject 
to accurate identification of adequate maximum torque capabilities. 
Some complication of model design is necessary to incorporate 
permanent magnet cores into load bearing wing structure and, princi¬ 
pally for this reason, the system is regarded with less favour than 
the alternative variant of the Spanwise Magnet scheme, that is, SIMs. 

The limited practical experience so far gained with SPMs indicates 
that the system is an excellent choice for SUMSBS and demonstrates 
the general practicality of Spanwise Magnet roll control (SPMs or 
SIMs). Cross couplings between model degrees of freedom due to 
the presence of the spanwise magnetizations are easily handled by 
suitable control algorithms and present no serious difficulty. 
Suspension quality with the SPM winged model was good. The spanwise 
magnetizations and the applied "roll” and other fields do alter 
(mostly adversely) the natural magnetic stiffness in several degrees 
of freedom but modest adjustments to the controller will be sufficient 
to restore the full stiffness and high stability currently achieved 
with unwinged models. Further computations of these effects would 
be appropriate, though the data presented in Section 4.3.1 contains 
the essential elements for a preliminary examination. 

The SIM system is currently viewed as the most promising system 
for LMSBS application. Model design is very straightforward, since 
wings may be machined from the solid with suitably chosen materials 
and the predicted torque capabilities are very high, limited princi¬ 
pally by current E/M technology. 

The level of confidence placed on the available calculations, of 
crucial importance of course, is currently high. However, extensive 
further and alternative calculations and a comprehensive demonstration 
at small scale are seen as essential before a LMSBS is committed to 
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the system. Such verification would necessarily include treatment 
of cross coupling effects, which are expected to be somewhat more 
pronounced than with SPMs. 

Subject to such computation and verification being satisfactory 
it is believed that careful study and optimisation of the roll E/M 
array will prove fruitful. The configurations used herein were 
chosen to be representative of likely LMSBS E/M arrays but were 
dimensioned entirely arbitrarily and are clearly deficient in many 
respects. If the ratio of useful field at, say, the model wing tip 
to the peak field level within or around the individual E/Ms is taken 
as a measure of performance, then the performance of the chosen SIM 
arrays is low. Adjustments to the E/M sizes, locations and winding 
cross-sections are expected to yield significant improvements. 

9.2 Extreme attitude testing 
9■2„1 Summary 

The largely intuitive realisation that the suspended model is 
influenced merely by the local applied magnetic fields and is other¬ 
wise unaffected, apart from the direction of the gravity vector, by 
its location or orientation within the suspension system, is easily 
illustrated by analysis. 

This realisation leads directly to the inherent ability of MSBSs 
to suspend and control models at any attitude and location, provided 
that fields and field gradients in the required senses and magnitudes 
local to the model can be somehow developed. 

It happens that the transformations between applied fields in 
balance axes and those seen in model axes are straightforward, but 
corresponding couplings for field gradients are relatively complex. 
Further difficulties are presented by the effects of the spatial 
variations of the applied fields and field gradients about the model's 
centroid. These are quite pronounced in systems such as SUMSBS 
where the E/M faces are in relatively close proximity to the model, 
but will tend to be of lesser magnitude in LMSBSs. Such effects 
render the adequate separation of model degrees of freedom, essential 
for the operation of contemporary control algorithms, quite complex. 

The prediction that SUMSBS would be usable, in the + configuration, 
from nominally 0° to 60° angle of attack (5-component control) was 
verified by achieving suspension of an axisymmetric model over this 
range. The inability of SUMSBS to operate around 90° angle of attack 
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in the + configuration and the difficulties apparent with the X 
configuration were traced to straightforward deficiencies in the 
E/M array. 

9.2.2 Discussion and conclusions 

It has long been claimed that one major advantage of wind tunnel 
MSBSs was the possibility of rapidly selecting near-arbitrary model 
attitudes. However, it is believed that no analytic or practical 
verification had previously been attempted. 

The problems presented by the gross changes in magnetic coupling 
between E/Ms and model, occurring during large changes of model 
attitude, are difficult but surmountable, the couplings being, of 
course, analytic or near-analytic and repeatable. 

Digital control is essential for sensible operation over wide ranges 
of angle of attack in order to permit adaption of the control loops 
to the changes in magnetic couplings and full digital control was, in 
fact, used for all demonstrations. 

Suspension quality with SUMSBS at the higher angles of attack was 
poor due to three major factors: 

(i) The decoupling of the axial position sensor was 
inadequate, resulting in coupling between pitch, heave 
and axial motions. 

(ii) The already approximate computed magnetic couplings 
(using FORCE), essential for adequate decoupling of model 
degrees of freedom in the controller, were further simplified 
before incorporation into that controller, for convenience. 

(iii) The natural magnetic stiffness of certain model 
degrees of freedom in the quasi-steady applied fields 
supporting the model's deadweight apparently became strongly 
negative at the higher angles of attack. Analysis and com¬ 
putation of this effect should be undertaken before further 
attempts at high angle of attack suspension are made. 

The inability of SUMSBS to operate above about 60° angle of attack 
in the -*• configuration and the poor performance predicted for the X 
configuration indicate that this simple fully symmetric 10 E/M array 
is unlikely to represent a good choice for a LMSBS. However, the 
shortcomings of the configuration, identified theoretically, are 
easily rectified by augmentation of the E/M array, indeed such 
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augmentation would be likely if SIM or SPM roll control was employed, 
and it is strongly believed that unlimited model attitudes, with full 
6-component control, will be easily achievable in realistic LMSBS 
configurations. The most difficult problem in such a system may, 
indeed, prove to be the provision of model position and other sensors 
capable of satisfactorily monitoring wide ranges of model attitudes 
and motion. 

9.3 Reliability 

9.3.1 Summary 

The widespread concern that the capital cost overhead incurred by 
a requirement for highly reliable operation would drive already 
expensive LMSBSs out of reach, has been partly dispelled. The two 
areas of greatest expenditure, that is the E/Ms (with associated 
plant) and the E/M power supplies, need only carry modest levels of 
redundancy and the complication of design, such as segmentation of 
E/Ms into several separate windings, is not severe. Indeed, 
segmentation of E/Ms and corresponding duplication of power supplies 
was featured in the only comprehensive LMSBS design study made to 
date (5), for reasons apparently unconnected with reliability. 

It is shown that the reliability of the remaining MSBS subsystems 
such as model position sensors and control systems may be adequately 
raised, albeit at the expense of some design complication, such as 
by the incorporation of at least two fundamentally different position 
sensing systems. 

9.3.2 Discussion and conclusions 

An unreliable LMSBS would be an expensive burden rather than a 
benefit. Any attempts to argue that a relatively high attrition rate 
of models in a "cheap" LMSBS could be justified is undermined by 
simple application of Murphy's Law, which indicates that the most 
expensive and intricate models would be destroyed before any useful 
aerodynamic data had been obtained. 

Though full technical presentation would be inappropriate here, 
preliminary demonstrations of hardware (E/M and E/M power supply) 
redundancy have already been made with SUMSBS. Specifically, with 
an axisymmetric model in steady suspension (wind-off and low angles 
of attack), ANY power supply, including those contributing towards 
the support of the model's weight, may be switched off with impunity. 
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The transient motion following the switch-off (or switch-on) is 
naturally controlled by the existing control algorithms, the model's 
weight, for instance, simply being redistributed amongst remaining 
operational E/Ms, and motion could be virtually eliminated if the 
controller was able to detect the "failure" of the power supply in 
real time, such as by monitoring the E/M currents. These demon¬ 
strations clearly confirm the existence and the powerful effect of 
Passive Redundancy, in this case not depending in any way on special 
configuration or adaption of a relatively simple MSBS. 

The key to reliable and economical LMSBSs is the EXPLOITATION 
rather than the SUPPRESSION of the many unique features of such 
systems. 

9.4 SUMSBS 

The commissioning of the new features, hardware or otherwise, of 
SUMSBS, has proved remarkably trouble free. 

Incorporation of digital control hardware and digitized control 
algorithms has been straightforward, with no particular difficulty 
in interfacing with the rest of the system. This success owes much 
to the care taken in the hardware design phase to create a system of 
high integrity, such as by the incorporation of complete electrical 
isolation between the control computer and the MSBS rig itself and 
the incorporation of sensible analogue filtering in position sensor 
preamplifiers and power supplies. The general performance of the 
controller is now as good as had ever been achieved previously with 
analogue systems. Changes or adjustments of control parameters may 
be made with precision and repeatability in seconds, rather than 
minutes or hours, and there now seems no reason to consider anything 
other than full digital control for a LMSBS or, indeed, any new MSBS. 

Considerable scope for improvements to the performance of the 
controller certainly remains, including the quality of suspension 
achievable with wind-on, the adaption of the controller to changing 
test conditions and the performance at extreme model attitudes. 

The new transistor switching power supplies, radically different 
from anything known to have been used previously in this application, 
are steadily amassing many hours of operation, to date completely 
faultlessly. The E/M loads (high inductance, low resistance) appear 
to be well suited to this type of supply. The level of radiated 
electrical noise, due to the high frequency switching, is much less 
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than had been feared, and has presented no serious difficulty. 
Accurate monitoring of E/M currents is not straightforward, but 
can be achieved with appropriate techniques. The independent 
and bipolar operation of each E/M has considerably enhanced the 
performance and versatility of the system, with no significant 
disadvantage yet apparent. 

The principal weakness in SUMSBS remains the position sensors. 
Significant changes in model geometry or gross changes of model 
attitude are only possible with mechanical realignment of the 
sensors. This is a serious inconvenience but in no way indicates 
a fundamental difficulty. Superior position sensing systems, 
optical or otherwise, are available and may be gradually incor¬ 
porated into SUMSBS in future years. 
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APPENDIX 1 


THE SOUTHAMPTON UNIVERSITY MAGNETIC SUSPENSION AND BALANCE SYSTEM 
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Al.l Introduction 239 
A1.2 E/M geometry and configuration 239 
A1.3 Power supplies 242 
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A1.5 Models 257 
A1.6 E/M current sensing 257 
A1.7 Control systems 259 


Al.l introduction 

The present status of SUMSBS hardware, represents only an interim 
phase in a more or less continuous program of development. This 
Appendix seeks merely to acquaint the reader with the current major 
features of SUMSBS, particularly those relatively recently intro¬ 
duced or updated. Fuller descriptions of certain hardware may be 
found in some earlier material, such as Refs. 11, 15. 

A1.2 E/M g eometry and configuration , 

The approximate dimensions of the main E/Ms and associated iron 
cores are shown in Figs.Al.l/2. All 8 E/Ms are nominally identical, 

6 in fact being carried over from earlier versions of SUMSBS. 
Previously, the iron cores had been configured rather as yokes (Fig. 
A1.2) with one "upstream" and one "downstream" E/M per assembly. 

The cores have.been modified by removal of the centre portions (shaded 
in Fig. A1.2) in order to reduce the fore-and-aft inductive couplings 
between associated E/Ms. This coupling had not previously been unduly 
troublesome, in fact the presence of the link augmented, for instance, 
the conventional "lift" field by some 30% due to the provision of 
return flux paths, but was considered undesirable for suspension at 
high angles of attack, or with spanwise magnet roll control, where 
more or less independent control of individual E/Ms is mandatory. 

Each E/M is wound with approximately 400 turns of 12 gauge enamelled 
copper wire, uncooled, and split into three sections normally 
connected in series. 
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Cores are approx. 

200 Stalloy laminations 
between 0.125 gauge 
plate facings. 

Overall thickness of 
cores is approx. 

2.9375 

Shaded portions 
now removed 


Fig. A1.2 Dimensions of main E/M iron cores. Dimensions in inches 


View on A 




Fig. Al.3 Dimensions of axial E/Ms of SUMSBS, showing configuration 



The twin axial E/Ms (Fig. A1.3) are somewhat similar in con¬ 
struction but are wound with thinner wire (14 gauge) split into 
two sections of approximately 500 turns each, normally connected 
in parallel. 

The E/M configuration is fully symmetric and is illustrated 
in Figs. A1.4 - A1.8. 

A1.3 Power Supplies 

All 10 E/Ms are currently independently controlled by CSR 
Contraves NC 421 servo controllers (Ref. 58 and Figs. A1.9 - Al.io). 
These controllers, rated at 21A continuous, 45A peak and 100V for 
four quadrant operation, were, in fact, designed for industrial 
D.C. servomotor drives but are well suited to this application. 

The fundamental mode of operation is pulse width modulation, with 
the ability to be driven unconditionally to 0% or 100% modulation. 
Since the switching frequency is nominally 5kHz the forward path 
voltage gain, neglecting the effects of current feedback, can be 
set flat to around 1kHz, representing a considerable advance over 
previous SUMSBS supplies (11). Fixed internal current feedback is 
in fact provided (Fig. Al.ll) and i.s considered to be extremely 
advantageous. Synthesis of model control algorithms is eased, 
countering the lengthening of the E/M time constants caused by 
removal of the series ballast resistors previously used (Section 
2.2.2 and Ref.11), the linearity of response is improved, output 
drift with temperature and time reduced and the effects of voltages 
induced in E/Ms by magnetic coupling from others is significantly 
attenuated. 

The power regulators which feed the controllers (Figs. A1.12, 
A1.13) themselves feature internal energy storage (on capacitors) 
totalling some 280 Joules overall, this being slightly greater than 
the total magnetostatic energy retained in the E/M array with all 
E/Ms at peak continuous current. Despite the fact that only some 
20% of this capacity is normally dynamically available, it is seen 
that since regeneration of load energy is allowed (four quadrant 
operation) the net power consumption of the complete system under 
normal conditions will tend to be low, being mostly due to resistive 
losses in E/Ms and power stages. Some 25 kW (continuous rating) 
of power is available from transformers and has been found to be 
more than adequate for all requirements. 
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A1.4 Position sensors 

The basis of the position sensing system used currently is an 
array of four analogue detectors of traditional form (Ref. 11 and 
Pigs. Al.14 - A1.21). Certain detail improvements have been made, 
most significantly to the light pickup assemblies where the first 
stage of amplification and filtering is now incorporated directly 
behind the sensor (Fig. A1.18). This is done in order to reduce 
the pickup of airborne electrical noise along the cabling between 
the MSBS rig and the main control system, this cabling now carrying 
relatively high level signals. 

High angle of attack capability is effected simply by appropriately 
rotating the aforementioned sensors broadly in the pitch sense (Figs. 
A1.21, A1.22). The range of sensitivity is such that sensing of a 
range of angle of attack from 0 to 60 degrees is practicable (see 
Section 7.8). 

Two sensors differ significantly from those described above, namely 
the axial position and roll attitude sensors. The axial position 
sensor was required to be capable of monitoring relatively large 
apparent axial motions, these being induced during large pitch 
rotations in the chosen configuration (Figs. A1.20, A1.21). A long, 
slender photodiode provides the necessary pickup range but would tend 
to be unduly sensitive to ambient light unless further measures are 
taken. Utilisation of a laser light source (lmW He - Ne) enables a 
narrow band optical filter of suitable wavelength to be installed in 
front of the photodiode, reducing the ambient light pickup by a factor 
of approximately 50 (measured). 

The roll attitude sensor comprises a laser (again lmW He - Ne), 
directed onto a mirror suitably located on the model, the reflected 
beam falling onto a UDT PIN-SC/50 two axis position sensing detector 
(Refs. 60, 61, and Figs. A1.20, A1.21). This device is essentially 
analogue in nature, providing four output signals which, when suitably 
processed, yield information as to the location of the centroid of 
the total incident light on the device's surface. The output can be 
made relatively insensitive to ambient light, partly by use of optical 
filters, and to variations in the intensity of the signal beam. The 
two axis capability of the device may be more fully exploited if 6- 
component high angle of attack suspension is attempted. 
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Fig. A1.14 Twin light source 
PLATES OVERLEAF . assembly 















Fig. A1.17 Light pickup 












Fig. A1.21 Schematic diagram of SUMSBS' position sensors 











A1.5 Models 


Recent recommissioning of SUMSBS has been performed almost 
entirely with simple cylindrical models of arbitrary aerodynamic 
shape (Fig. A1.23) fitted with existing cylindrical Alnico V 
magnetic cores. It has been found that these cores can now be 
demagnetized following model flyaway by certain combinations of 
applied fields, necessitating automatic shut-off of all E/Ms 
whenever position sensor signals indicate that the model is departing 
from their range. Even with this protection, the long term cali¬ 
bration of such a core must be regarded as dubious, recommending 
future changes in the core material. 

Exploratory 6-component suspension with spanwise magnet roll 
control has been achieved with a model basically as described above, 
but crudely adapted to carry arbitrarily dimensioned Alnico V wing 
cores (Fig. A1.24). The initial configuration of the roll angle 
sensor (Figs. A1.20, A1.21) necessitated a somewhat awkward location 
of the model-mounted mirror but this is, of course, not a fundamental 
shortcoming of this general technique. 

A1.6 E/M current sensing 

Since the output voltage from the CSR power supplies appears as a 
pulse width modulated, 5kHz, 1O0 V peak-to-peak square wave, with 
neither output terminal near ground potential (Fig. A1.25), special 
measures are necessary to acquire accurate signals from conventional 
shunts in the E/M supply leads. A shunt signal typically consists, 
under these conditions, of a small (- IV) differential component 
representing load current, together with a common-mode AC signal of 
approximately 100 V peak-to-peak (Fig. A1.25). Fully isolated signal 
amplifiers are necessary to permit safe handling of the peak common¬ 
mode voltages and high rejection of this signal component is essential 
for accurate resolution of load current. 

Configuring dual series shunts as signal sources for a linear 
isolation amplifier (Burr-Brown 3650 HG, Ref. 62), as shown in Fig. 
A1.2G, causes the "common-mode" signal to appear, in fact, as an 
"isolation" signal, that is, wholly between the input common and true 
ground. Adequate rejection capability of this signal is available 
from the chosen class of device, specifically some 140 dB at D.C. 
falling to approximately 90 dB at 5kHz, all figures given for a device 
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Fig. A1.26 Isolated shunt signal acquisition 






voltage gain of 10, from Ref. 62. In fact, the electrically noisy 
environment in which the devices are operated causes some degradation 
of this performance, mostly with regard to output noise, but modest 
filtering of the output provides a satisfactorily clean and effec¬ 
tively real-time signal for data acquisition purposes. 

Al.7 Control systems 

The control loop is implemented wholly on a digital computer 
(PDP - 11/34) with only limited signal conditioning performed 
externally. Local 12 bit A/D and D/A converters communicate with 
the host computer via two fully isolated, high speed, uni-directional 
16 (+ control) bit data buses, the computer in fact being located 
some 90 feet from SUMSBS itself (Figs. A1.27 - A1.28). Treatment of 
the control algorithms used is outside the scope of this report but 
all suspension necessary for the work reported herein has been accom¬ 
plished with relatively simple digital versions (59) of traditional 
analogue algorithms (11), absolute performance of such controllers 
thus being comparable to that achieved by the previous analogue 
systems. 











A/D subsystem 



Fig, A1.28 Sche matic diagrams of A/D and D/A 
converter subsysterns. 
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A2.1 Introduction 

The forces and torques acting on isolated bodies whose 
magnetization distribution is known, are amenable to approximate 
numerical calculation where the applied field can be somehow founds 
The program FORCE is an example of one simple method, in fact 
carrying out elementary numerical integration of the relevant 
vector products (Eqns. 3.1) over the volume of the core. 

The core is converted to an array of dipoles, each dipole 
representing the magnetization of suitably chosen model elements. 

Such representation of the core is only straightforward if the core 
magnetization is essentially uniform, whence each dipole can be 
placed at the centroid of the corresponding element. Uniform 
magnetization is frequently assumed in the case of ReCo materials, 
where the material's high coercive force precludes self-demagnet¬ 
ization effects. However, computations may be approximately 
applied to other cases, such as slender axially magnetized iron cores. 

Computation of applied E/M fields is generally only straight¬ 
forward where the E/Ms are air-cored. In the case of MSBSs where 
the E/Ms are relatively remote from the model it is usually adequate 
to represent the E/Ms as assemblies of straight line wire elements, 
whence the external field is found by simple summation of the contri¬ 
bution of each element. The field of individual elements is 
easily found, such as by use of the Biot-Savart law (63 and various). 

A program segment that had been shown to reliably compute the field 
from straight line wire elements was extracted from the program 
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TABLE, courtesy Massachusetts Institute of Technology (64, 65), 
and is used in FORCE essentially unmodified. 

A2.2 Program structure 

The fundamental program structure is shown in Fig. A2.1. The 
program is heavily modularised to permit implementation on a 
minicomputer (PDP 11-34) with relatively limited available memory. 

Main model and E/M data storage is in virtual arrays. Memory 
requirements are comfortably below the available 28K words (main 
code) and 32K words (virtual arrays) with 1000 model and 1000 E/M 
elements (see Fig. A2.2). 

The main (root) segment of the program simply performs the 
function of calling the operational subroutines in sequence, this 
approach allowing these subroutines to be overlaid in memory, as 
Fig. A2.2. 

A2.3 Subroutine descriptions 

MODIN controls the filling of the model core dipole array via calls 
to user named datafiles or the subroutines FUSIN or WING. Model 
core data can be output to the console device or to user named data¬ 
files. The model can be located arbitrarily in space and looping 
can be called whereby one locating parameter, say pitch attitude, is 
stepped by chosen increments over a chosen range, with force and torque 
calculations repeated at each step. 

The fundamental equations controlling the creation of the dipole 
array are as follows: 


M = / M 6V 


where the integration is over 
volume of a single element. 


leading to modified forms of Eqns. 
F ' = Z (M* . V ) H. 

— l—i —l 


(A2.1) 


FUSIN permits representation of simple axisymmetric fuselage geometries 
as shown in Fig A2.3 Where FORCE is used simply to calculate fields 
due to particular E/M arrays, FUSIN fills the CORE array with dummy 
model elements of zero magnetization (Fig. A2.4). 
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Pig. A2.1 Fundamental block structure of program FORCE 



Approx. 
memory 
location 


Fig. A2.2 Memory overlay structure of program FORCE 



















Each point represents 
the location of a 
dipole, in turn 
representing the 
magnetization of a 
single element 


Fig. A2.4 Field grid. Fuselage input option = 2. Zero magnetization 


all element; 



WING permits representation of simple slab wings, magnetized 
symmetrically (SPMs, see Section 4), as shown in Fig. A2.5 

COIL creates an array of straight line wire elements representing 
the assembly of E/Ms. Symmetry is used extensively to reduce the 
quantity of user input per program run and the program’s main 
memory requirement. Data is input for a single "master" E/M, 
which can later be reflected and/or rotated in certain ways to 
generate other E/Ms in that particular "set" (Fig. A2.6). Up to 
five E/M sets may be input sequentially, permitting analysis of 
relatively complex MSBS configurations. Four types of E/M may be 
specified, single wire loops (Fig. A2.7), E/Ms effectively composed 
of straight blocks of conductor of quadrilateral cross section 
(Fig. A2.8) and "lateral" or "axial" pseudo-circular E/Ms (Figs. 

A2.9, A2.10). Pseudo-circular E/Ms are reduced to assemblies of 
circumferential straight line elements, positioned at suitable 
radii so as to create a field at the centre of each loop identical 
to that created by the equivalent circular loop carrying the same 
current. E/M data may be output to user named datafiles or to the 
console device. Data for a complete array (several E/M sets where 
required) may be input from a user named datafile. 

E/M ampere-turn products are input individually at the console 
device in a standard sequence (Fig. A2.ll). 

FACE generates the co-ordinates of straight line wire elements 
distributed over the cross section of a block of conductor, per¬ 
mitting E/Ms to be input as a series of co-ordinates representing 
the interfaces where blocks of conductor join as Fig. A2.8. 

FACE is derived from the TABLE subroutine COORS (MIT, (65) ), 
but differs in two important respects. Firstly, since COORS only 
generated wire element distributions that correctly simulated 
uniform current density over the winding cross section when the 
interfaces were quadrilateral, FACE generates the fourth corner of 
each interface internally. Secondly, COORS used the input inter¬ 
face co-ordinates to define directly the location of the outermost 
wire elements. FACE uses the interface co-ordinates to describe 
an outer envelope, inside which all elements are distributed. 

Thus, identical input data to COORS and FACE will result in different 
wire element distributions. 
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DLONG 



Fig. A2.5 Slab wing geometry. Wing input option = 1 














''ANGLE I 



NDIVY or 
NDIVX 

* • \ 

NDIVR 


Section on AA showing 
wire loop distribution 

As shown: 

NDIVR=1, NSEG=12 
NDIVY(lateral) =2 
NDIVX(axial) =2 


AXIAL and LATERAL 


Fig. A2.9 and A2.10 Lateral and Axial Pseudo-circular E/Ms. E/M input 






CALC controls the calculation of E/M fields, via calls to the 
subroutine FIELD, summation of model forces and torques and the 
routing of calculated data. The position and orientation of 
model elements is updated using previously input data (Fig. A2.12). 
FIELD then provides the field and field gradient components at 
each model element in turn, the relevant vector products of field 
and magnetization being calculated and summed. Model forces and 
torques are expressed about the model's own co-ordinate origin, 
in translated but unrotated axes (Fig. A2.13). The calculated 
field at each model element may be output to datafile or the 
console device. Force and torque data may be routed likewise. 
Calculations are repeated if stepping of the model position or 
orientation had been previously selected. 

FIELD calculates the field and field gradient components at the 
single points required, due to the specified E/M array. Individual 
E/M set calculations are made first for the master E/M, then for 
E/Ms generated by the chosen reflection/rotation symmetries. Set 
calculations are repeated for all specified sets with the final 
field data being the sum over all E/Ms of all sets. 


273 - 



MSBS fixed axes : 
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A2.4 


GENERAL INSTRUCTIONS FOR USERS OF PROGRAM FORCE *******JK**)*:*sJo{o}c>k**3KJf:* 
THE MAXIMUM NUMBER OF MODEL ELEMENTS CURRENTLY PERMITTED IS 1000 

the number of elements required for a particular model may be calculated 
as follows: 

N=<NOX+i>*<NOY+l>*<NOZ+i> - FIELD GRID 

=NL*NA*NR — CYL. FUSE. 

+NA*NR*NA1#2 -ELLIP. ENDCAPS 

+NS*NC*NT*2 ■-■■■ ■-SLAB WINGS 


THE MAXIMUM NUMBER OF E/M WIRE ELEMENTS CURRENTLY PERMITTED IS 1000 
THE NUMBER OF ELEMENTS IN A PARTICULAR E/M ARRAY MAY BE CALCULATED 


as follows: 

M=SUM 


OVER ALL SETS OF: 
NINT 

NINT*MM*NN 

NDIVR#NDIVY*NSEG 

NDIVR*NDIVX#NSEG 


- SINGLE LOOP 

- INTERFACE ARRAY 

— LATERAL CIRCULAR 
-— AXIAL CIRCULAR 


ONLY THE MASTER E/M IN EACH SET AFFECTS THE VALUE OF M. E/MS GENERATED 
BY SYMMETRY OPTIONS DO NOT REQUIRE FURTHUR STORAGE ALLOCATION. 

THE NUMBER OF MODEL AND E/M ELEMENTS SHOULD ALWAYS BE CHOSEN TO BE AS 
SMALL AS WILL YIELD ACCEPTABLE RESULTS. CHOICE OF UNNECESSARILY LARGE 
NUMBERS OF ELEMENTS WILL RESULT IN INCONVENIENTLY LONG EXECUTION TIMES 
AND PERHAPS IN THE ACCUMULATION OF ROUNDING ERRORS. 

COMPUTAIQN TIME* TAKEN FROM THE COMPLETION OF THE LAST USER 
INPUT TO FINAL PROGRAM COMPLETION* NEGLECTING HARDCOPY PRINT TIMES* 

HAS BEEN FOUND TO BE APPROXIMATELY: 

T=N*MM*<STEPS+1>*0.01 

WHERE MM-TGTAL NUMBER OF E/M ELEMENTS SUMMED OVER ALL 
E/MS OF ALL SETS. 

PDF- 11-34 WITH HARDWARE FLOATING POINT OPTION* 
RT~11 OPERATING SYSTEM* GIVEN IN SECONDS. 

the: STANDARD FORMAT FOR THE TITLES OF USER NAMED DATAFILES UNDER RT-11 

is as follows: 

******.DAT 

WHERE ****** IS AN ARBITRARY 6 CHARACTER ALPHANUMERIC NAME 


ALL INPUT/OUTPUT IS IN SI UNITS* SPECIFICALLY: 

DIMENSIONS IN METRES 

FIELDS IN TESLA 

FIELD GRADIENTS IN TESLA/METRE 

ANGLES IN DEGREES (CONVERTED INTERNALLY TO RADIANS) 


275 - 



A2.5 


USER INPUT VARIABLES FOR PROGRAM FORCE ***»C5»:**5lc>|c*>k****>lo|c>lo|c******>lc*>lc*** 

REFER TO FIGS.A2.1 TO A2.13 FOR DEFINITION OF VARIABLES 
NOT CITED BELOW 

MODEL GEOMETRY OPTION (MOPT) 

ab where: 

A-WING INPUT OPTION 

0 NO WING 

1 SLAB WING 

B=FUSELAGE INPUT OPTION 

0 FROM DATAFILE (WHOLE MODEL) 

1 CYLINDRICAL FUSELAGE 

2 FIELD GRID 

3 AS 14-ELLIPSOID ENDCAPS 

EG. MOPT”11 SPECIFIES CYLINDRICAL FUSE. WITH SLAB WINGS 
MOPT=l SPECIFIES CYLINDRICAL FUSE. ONLY 

MODEL OUTPUT OPTION (IOUTOP) 

0 NONE 

1 TO CONSOLE 

2 TO DATAFILE 

MODEL DISPLACEMENT OPTION (MDISPO) 

0 NONE 

1 SINGLE SPECIFIED DISPLACEMENT 

2 LOOPING OF CHOSEN PARAMETER FROM SPECIFIED INITIAL 

DISPLACEMENT 

-2 AS 2 BUT FROM NULL INITIAL DISPLACEMENT 

WHERE MDISP0=2»-2 

STEPS NO. OF INCREMENTS IN CHOSEN PARAMETER NOT INCLUDING 
START POSITION 

PARM 1-6 SELECTING X>YrZ»YAW,RITCH*ROLL DISPL. RESPECTIVELY 
SIZE SIZE OF STEFS (METRES OR DEGREES) 

NO. OF E/M SETS (NSET) 

NO. OF INDEPENDANTLY SPECIFIED SETS OF E/MS 
0 INPUT COMPLETE ARRAY FROM DATAFILE 

(E/M) INPUT OPTION (INPOPT) 

1 SINGLE LOOP 

2 INTERFACE ARRAY 

3 LATERAL PSEUDO-CIRCULAR 

4 AXIAL PSEUDO-CIRCULAR 

(E/M) SYMMETRY (ISYMM(I*J>> 

1 NONE (MASTER E/M ONLY) 

2 REFLECTION IN Y»Z AND XrZ PLANES 

3 1 ROTATION ABOUT X AXIS - BY 180 DEGREES 


4 3 ROTATIONS 90 

5 5 60 

6 7 45 

7 REFLECT IN Y»Z PLANE 

8 REFLECT IN XrY X,Z AND Y>2 PLANES 

? AS 3 WITH REFLECTION IN Y»Z PLANE 

10 AS 4 YrZ 

11 AS 5 Y»Z 

12 AS 6 YfZ 


E/M OUTPUT (IOUTOP) 

0 NONE 
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INTERFACE ARRAY COORDINATES ONLY 
AS 1+WIRE ELEMENT COORDS. 

WIRE ELEMENTS TO DATAFILE 


E/M AMPERE TURNS 

TOTAL AMPERE-TURN PRODUCT FOR EACH E/M. INPUT IN STANDARD SEQ. 

(WIRE ELEMENT CURRENTS DERIVED APPROPRIATELY) 

THE SIGN FOR CURRENT IN SYMMETRY-DERIVED E/MS IS SUCH THAT 
THE COMPONENT OF FIELD RADIALLY IN TOWARDS THE ORIGIN IS 
MAINTAINED CONSTANT (IN MAGNITUDE. NOT DIRECTION) WITH 
CONSTANT E/M CURRENTS 

(DATA) OUTPUT OPTION (.TOUT) 

=AB where: 

A=FIELD LISTING OPTION 
0 NONE 

1 TO CONSOLE 

2 TO DATAFILE (2FIELD.DAT) 

B=FORCE/TORUUE OUTPUT OPTION 

0 TO CONSOLE 

1 TO DATAFILE (2F0RCE.DAT) 

EG.21 ROUTES FIELD AND FORCE/TORUUE DATA TO FILES 

0 ROUTES FORCE/TORUUE DATA TO CONSOLE t FIELD DATA LOST 


Refer to subroutine principal variable lists for description of 
detail input variables. 
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FORTRAN IV 

C ROOT SEGMENT OF PROGRAM FORCE >{c * >{ofc>Jc Jfc * * >K>K « «^ >5? 5{< 5{< 5i< * ^ Jit * >{« >K « 

C NO VARIABLES FUNCTIONAL IN THIS SEGMENT 
C 

0001 VIRTUAL CX<1000r2>>CY'<1000>2)?CZ(1000 r 2) ? CORE<1000r6) 

0002 DIMENSION DISP (7) ? FORCE ( 9) ? ISYMM(5 r 3 > t CURR<16*5) » J.DISP < 3 > 

C 

C MODEL DATA INPUT 
C 

0003 CALL MODIN(CORE t DISP»IDISP > 

C 

C COIL DATA INPUT 
C 

0004 CALL COIL (CX y CY r CZ y ISYMM y CURR) 

C 

C FORCE ? TORQUE AND FIELD CALCULATION 
C 

0005 CALL CALCCCXyCYrCZyCOREs'DISPi’ISYMMyCURRi’IDISP) 

0006 STOP ' HAVE A NICE DAY ' 

0007 END 


A2a 6 Subroutine listings 



FORTRAN IV 


V02, 
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0001 SUBROUTINE MODIN(CORE,DISP,IDISP) 

0002 VIRTUAL C0RE(1000»6> 

0003 DIMENSION DISPI7)>IDISP(3) 

0004 PI=3.141592654 

C 

C ZEROISE CORE AND DISPL. ARRAYS 
C 

0005 DO 10 J=1 »3 

0006 10 IDISP!J)=0 

0007 DO 20 1=1,6 

0008 DO 20 J=l,1000 

0009 20 CORE!J»I)=0. 

0010 NE=0 

C 

C MODEL INPUT SECTION 
C 

0011 30 TYPE 40 

0012 40 FORMAT!' MODEL GEOMETRY OPTION?') 

0013 READ!5»*»ERR=30> MOPT 

0014 IF(MQPT.EG.O) 60 TO 60 

0016 MGPT2=M0PT/10 

0017 M0PT1=M0PT-M0PT2*10 

0018 CALL.. FUSIN< MOPT 1 >NE> CORE) 

00.19 CALL WING(M0PT2»NE*C0RE> 

0020 IDISP(3)=NE 

0021 GO TO 50 

C 

C MODEL INPUT FROM DATAFILE 
C 

0022 60 TYPE 70 

0023 70 FORMAT<' FILENAME?'/) 

0024 CALI... ASSIGN< 1 , *-l , 'OLD' 

0025 DEFINE FILE 1 (1001,12»UfNREC) 

0026 DO 00 1 = 1 * 1000 

0027 00 READ! .1' I) (CORE (I»J)»J=1»6> 

0028 READ(l'lOOl) IDISP(3) 

0029 CL0SE<UNIT=1) 

C 

C MODEL OUTPUT OPTION 
C 

0030 50 TYPE 90 

0031 90 FORMAT!' MODEL OUTPUT?') 

0032 READ ! 5»# * ERR=50) IOUTOP 

C 

C MODEL OUTPUT 
C 

0033 GO T0(100fll0fl20)»(IOUTOP+1) 

0034 110 TYPE 130 

0035 TYPE 140 

0036 .130 FORMAT!' MODEL CORE COORDS. AND MAG.') 

0037 140 FORMAT(4X r1HX»7X»1HY * 7X * 1HZ * 8X»6H(J*V)X*8X>6H(J#V)Y 

lv8X*6H(.J*V)Z) 

0038 DO 150 1=1 r IDISP<3) 

0039 150 TYPE 160,(CORE!I,J),J=1 ,6) 
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0040 160 FORMAT(3F8.4.3E14.4) 

0041 GO TO 100 

0 

C MODEL OUTPUT TO DATAFILE 

C 

0042 120 TYPE 70 

0043 CALL ASSIGN(2. .-1.'NEW'» » ) 

0044 DEFINE FILE 2 (1001.12.U.NREC) 

0045 DO 170 1=1r1000 

0046 170 WRITEI2'1)(CORE <I.J)»J=1.6) 

0047 WRITE(2'1001)IDISP<3) 

0048 CLOSE(UNIT=2) 

C 

C GET MODEL POSN. AND ATTITUDE 
C 

0049 100 TYPE 180 

0050 180 FORMAT<' MODEL DISP. OPTION?') 

005.1 READ< 5»He.ERR= 100) MDISPO 

0052 IF (MDISF’Q.GT .0) GO TO 190 

C 

C MDISP0=0.-2 
C 

0054 DO 200 J=1.7 

0055 200 DISP(J > =0 * 

0056 IF(MDISPO.EQ.-2) GO TO 220 

0058 GO TO 210 

C 

C MDISP0=1 
C 

0059 190 TYPE 230 

0060 230 FORMAK' MODEL DISP. <X»Y» Z >» ORI (YAW. PITCH»ROLL)?' ) 

0061 READ(5»#»ERR=.t90) (DISP(J).J=1>6) 

0062 DISP<7)=0. 

0063 DO 240 J=4»6 

0064 240 DISP<J)=DISP<J)*PI/1B0. 

0065 IF < MDISPO. Etl.l) GO TO 210 

C 

C MDISF’0=2. -2 
C 

0067 220 TYPE 250 

0068 250 FORMAK' STEPS. PARAMETER. SIZE?' ) 

0069 READ<5.#»ERR=220) IDISP<1).IDISP < 2 >.DISP < 7) 

0070 IF(IDISP(2).LT.4) GO TO 210 

0072 DISP (7) "-DISP (7) *P.I/180. 

0073 210 RETURN 

0074 END 
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0001 

0002 

0003 

0004 


OOOS 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 


0016 

0017 

0018 

0019 

0021 


0022 

0023 

0024 

0025 

0026 


SUBROUTINE FUSIN (MOF'T 1 > NE » CORE ) 

VIRTUAL CORE!1000»6) 

PI-3*141592654 

GO TO <30»20» 10»20) » (MO.PTl+1) 

C 

C FIELD GRID SECTION 
C 

10 TYPE 40 

40 FORMAT <' XOrYOrZO, XINC»YINC»ZINC> NOXuNOYuNOZ?") 

READ(5r«) X0»YO»ZO>XINC»YINC)>ZINC>NOXnNOYi>NQ2! 

DO SO IX=1»NOX 
DO 50 I Y~i »NOY 
DO 50 IZ-lrNOZ 
NE-NE+1 

CORE (NE»1 >--X0+< IX-1>#XINC 
C0RE(NE»2)=Y0+( IY--1 )#YINC 
50 CORE (NE , 3)*Z0+ (IZ-1)*ZINC 

GO TO 30 
C 

C CYLINDRICAL SECTION 
C 

20 TYPE 60 

60 FORMAT <' CYLIND.FUSE.tFLEN»FRAD»NA»NR?') 

READ < 5 f %) FLEN > FRAD,NL » NA t NR 
IF< NL * EQ♦0) GO TO 80 
ANGL=PI*2./NA 
C 

C LOOP FOR RADIUS 
C 

DO 70 IRAD=0»NR-1 
R2~<IRAD+1>*FRAD/NR 
Rl=IRAD*FRAD/NR 
XAREA=ANGL*< R2*R2-R1*R1)/2. 

D=4.*<R2*#3--Rl**3)*SIN<ANGL/2. )/(3.-#<R2*R2-Rl*Rl >"*ANGL) 
C 

C LOOP FOR ANGLE 
C 


0027 

0028 


0029 

0030 

0031 

0032 

0033 

0034 

0035 


0036 


DO 70 IANG=0»NA“1 
ANG"1 ANG#ANGL 
C 

C LOOP FOR AXIAL POS. 

C 

DO 70 ILEN=0*NL-1 
NE=NE+1 

CORE < NE»1) "-FLEN/2. + < IL.EN+0.5) #FLEN/Nl. 
CORE(NE,2)=D*COS(ANG) 

CORE(NE»3>«D*SIN(ANG) 

CORE (NE>4) ^XAREAXcFLEN/NL 
70 CONTINUE 

C 

C END OF CYLINDRICAL SECTION 
C ELLIPSOID ENDCAPS SECTION 
C 

80 l'F<M0PT1 * NE >3) GO TO 30 
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0038 

0039 

0040 

0041 


0042 

0043 

0044 

0045 

0046 


0047 

0048 


0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


0059 

0060 


0061 

0062 


TYPE 90 

90 FORMAT(' ELLIP.ENDCAPSJARfNAI?') 

READ<5»*> AR»NA1 
THETA=PI/<NA1*2> 

C 

C LOOP FOR ANGLE FROM AXIS 
C 

DO 100 ITHETA=OfNA1-1 
THETA1=THETA#1THETA 
THETA2-THETA*<ITHETA+1) 

D-CQS<THETA1)-COS < THETA2) 

T=C0SCTHETA1*2.>~C0StTHETA2*2*) 

C 

C LOOP FOR ANGLE AROUND AXIS 
C 

DO 100 IANG=0»NA-1 
ANG ; ~I ANG&ANGL 
C 

C LOOP FOR FRONT/BACK AND PSEUDO RADIUS 
C 

DO 100 IFRBAK=-1»1*2 

DO 100 IRAD==0»NR-1 

R1 ~ IR ADStFRAD/NR 

R2=(IRAD+1)*FRAD/NR 

R-- < R2*#4-R1**4 ) / (R2**3-R1**3 ) 

NE"=NE+1 

C=3.*SIN < ANGL/2.> *R* <THETA+T)/(4.JKANGL^D) 

CORE(NE f 1) = IFRBAK*(-FLEN/2.-AR*R#T*3./(D*16.)) 
C0RE<NEr2)=C*C0S<ANG> 

C0RE<NEr3)=C*SIN(ANG> 

CORE < NE f 4) - < R2*#3-R1**3) #AN6L*AR:l(Ii/3. 

100 CONTINUE 
C 

C END OF ELLIP. ENDCAPS 
C 

30 RETURN 
END 


-282 




FORTRAN IV 


002.5 


PAGE 001 


0001 SUBROUTINE WING(M0PT2*NE,C0RE> 

0002 VIRTUAL C0RE<1000r6> 

0005 PI=3.141592654 

0004 IF<M0PT2.EQ.0) GO TO 60 

C 

C SLAB WING INPUT 
C 

0006 10 TYPE 20 

0007 20 FORMAT(' SLAB WING *1/2SPAN > 1/2GAP t CHORD t THICK t TAPER t SWEEP?') 

0008 READ(5 * *»ERR=10) SPAN FRAD t CHORD * THICKrTAPER»SWEEP 

0009 30 TYPE 40 

0010 40 FORMAT(' DLONG»DVERT»NS»NC»NT»IWIDIR?') 

001.1 READ (5 1 * t ERR=30> DLONG»OVERT r NS , NC * NT r IWIDIR 

0012 - SWEEP-SWEEPfcP1/180. 

0013 SINSWE~SIN<SWEEP) 

0014 COSWE=COS<SWEEP) 

0015 SPAN2-SPAN-FRAD 

C 

C LOOP FOR PORT AND STB. WINGS 
C 

0016 DO 50 IWING=-lfl»2 

C 

C LOOP FOR SPAN AND THICKNESS 
C 

0017 DO 50 ISPAN-0,NS-1 

0018 DO 50 ITHICK=0»NT-1 

C 

C CALCULATE CORE DIPOLES 
C 

0019 DO 50 ICH0RD=0»NC-1 

0020 NE--NE+1 

0021 T= (< .1. -TAPER) * < ISPAN+1) ~NS) / ( < 1. -TAPER) *ISPAN-NS) 

0022 TYPO-FLOAT(ISPAN)/NS+ <1.-2.* <1-T)/3.)/< <1.+T)*NS) 

0023 CHORDL=CHORD*<1.-<1.-TAPER)#TYPO) 

0024 V0L=CH0RD#THICK*SPAN2* <1.+T)/< 2. #NC#NT*NS ) 

0025 CORE(NE.1)-DLONG-CHORDL* < 0.5~ < 0.5+ICHORD)/NC)-SPAN2*TYP0 

1.YSINSWE/C0SWE 

0026 CORE(NE r 2)a(SPAN2*TYPO+FRAD)*IWING 

0027 CORE<NE.3)=DVERT-THICK*< 0.5-< 0.5+1THICK)/NT) 

0028 CORE(NE * 4)~VOL#SINSWE#IWIDIR 

0029 C0RE<NEf5)=VQL*C0SWE*IWIDIR*IWING 

0030 50 CONTINUE 

0031 60 RETURN 

0032 END 
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0001 SUBROUTINE COIL(CX »CY,CZ,ISYMM,CURR) 

0002 VIRTUAL CX!1000,2),CY(1000,2),CZ(1000,2) 

0003 DIMENSION XX!4,10),YY(4,10)»ZZ!4,10),ISYMM(5,3),CURR<16,5) 

0004 DIMENSION IN0C0(12> 

0005 DATA INOCO/1,4,2,4,6,0,2,8,4,8,12,16/ 

0006 NSET=0 

0007 IUIRE=1 

000(3 I0UTQP--0 

0009 PI=3.141592654 

0010 DO 10 J=1,5 

0011 ISYMM <J,1)~0 

0012 DO 10 1=1,16 

0013 10 CURR(I,J)=1. 

0014 20 TYPE 30 

0015 30 FORMAT<' NO. E/M SETS?') 

0016 READ!5,*,ERR=20) NSET 

0017 IF(NSET.GT.O) GO TO 40 

C 

C COMPLETE SYSTEM DATAFILE INPUT 
C 

0019 TYPE 50 

0020 50 FORMAT(' FILENAME?'/) 

0021 CALL ASSIGN(1»»-1»'0LD'»,> 

0022 DEFINE FILE 1 !1003,12,U,M) 

0023 DO 60 1=1,3 

0024 60 READ!1'I)(ISYMM!J,I),J=1,5) 

0025 DO 70 K=1,1000 

0026 70 READ(1'K+3XCX(K,J),CY(K,J>,CZ!K»J),J=1,2> 

0027 CL0SE(UNIT=1> 

0028 DO 80 1=1,5 

0029 IF(ISYMM!1,1).EQ.O) GO TO 80 

0031 NSET=NSET+1 

0032 80 CONTINUE 

0033 GO TO 90 

C 

C E/M SET LOOP START 
C 

0034 40 DO 100 ISETI=1,NSET 

C 

C GET MASTER E/M GEOMETRY FOR THIS SET 
C 

0035 TYPE 110,ISETI 

0036 110 FORMAT!' E/M SET',15) 

0037 120 TYPE 130 

0038 130 FORMAT!' INPUT OPTION,SYMMETRY OPTION?') 

0039 READ(5,>R»ERR=120) INPOPT, ISYMM! ISETI , 1) 

0040 GO TO (140,150,160,170) INPOPT 

C 

C SINGLE LOOP INPUT SECTION 
C 

0041 140 TYPE 180 

0042 1(30 FORMAT!' NO. ELEMENTS?') 

0043 READ(5»#,ERR=140) NINT 

0044 190 TYPE 200 
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0045 200 FORMAT (' ELEMENT END POINTS?') 

0046 READ<5,*rERR=190> (CX(I*l>»CY<I»t) »CZ< I»l)»I~li>NINT> 

0047 NEND-NINT+IWIRE-1 

0048 DO 210 J==IWIRE»(NEND-1) 

0049 CXIJr2I=CX(J+l»l) 

0050 CY(J»2)=CY(J+l»1) 

0051 210 CZ(J»2>=CZ(JH»1> 

0052 CX (NEND *2>-CX< IWIREn 1) 

0053 CY < NEND»2)-CY(IWIRE*1) 

0054 CZ( NEND »2)==CZ< IWIRE* 1) 

0055 ISYMMJISETI>3)~1 

0056 GO TO 220 

C 

C INTERFACE ARRAY INPUT SECTION 
C 

0057 150 TYPE 230 

0058 230 FORMAT<' NUMBER OF LAYERS» LOOPS/LAYERi> INTERFACES?') 

0059 READ(5y#»ERR--150) MM»NN»NINT 

0060 240 TYPE 250 

0061 250 FORMAT<' INTERFACE COORDINATES?') 

0062 DO 260 J«1»NINT 

0063 260 READ(5r*rERR=240) <XX(I»J)rYY<I»J)fZZ<I»J)»1=1»3) 

0064 NEND=MM*NN*NINT+IWIRE-1 

0065 CALL FACE ININTfMM fNN»XX»CX»IWIRE) 

0066 CALL FACE(NINT t MM,NNrYY,CY>IWIRE) 

0067 CALL FACE(NINT * MM»NN > ZZ»CZ»IWIRE > 

0068 ISYMMIISETIf3)=MM*NN 

0069 GO TO 220 

C 

C END OF INTERFACE ARRAY 
C LATERAL CIRCULAR E/M SECTION 
C 

0070 160 TYPE 280 

0071 280 FORMAT(' DX> DY»RADI»RAD2>DYi»DY2 >NDIUR >NDIVY rNSEG t ANGLE?') 

0072 READ< 5 * #>ERR-160) DX r DY r RADI t RAD2 1 DY1>DY2»NDIVR* NDIVY * NSEG t ANGLE 

0073 CQSANG=COS< ANGLE#PI/180«> 

0074 SINANG-SIN <ANGLEKPI/180.) 

C 

C LOOP FOR AXIAL DEPTH 
C 

0075 DO 290 I=1»NDIVY 

0076 Y-DY1+(1-0.5)#(DY2-DY1)/NDIVY 


C LOOP FOR RADIAL THICKNESS 
C 

0077 DO 290 M~1»ND1VR 

0078 R-RAD1+(M-0.5)*(RAD2-RAD1)/NDIVR 

0079 R»R*NSEG*SIN < PI/NSEG >/< PDKCQS ( PI/NSEG)') 

C 

C LOOP AROUND CIRCUMFERENCE 
C 

0080 DO 290 J~1 i NSEG 

0081 K — 1 < 1-1 >*NDIVR+M-1 >*NSEGhJ+IWIRE-l 

0082 ANG-PDK2./NSEG 
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0083 ZT1~SIN<ANG*<J~1))*R 

0084 ZT2=SIN(ANG*J)*R 

0085 CX(K * 1)=DX+COS< ANG*(J-l))*R 

008A CX<K»2)=DX+C0S<ANG*J)*R 

0087 CY(K»1)=Y*C0SANG~ZT1*SINANG+DY 

0088 CY < K r 2) ==Y*C0SANG-ZT2*SINANG+DY 

0089 CZIKr 1>=Y*SINANG+ZT1*C0SANG 

0090 CZ<K»2)=Y*SINANG+ZT2*C0SANG 

0091 290 CONTINUE 

0092 NEND==NSEG#NDIVR*NDIVY+IUIRE-1 

0093 ISYMM<ISETIr3)=NDIVR*NDIVY 

0094 GO TO 220 

C 

C END OF LATERAL CIRCULAR E/M 
C AXIAL CIRCULAR E/M SECTION 
C 

0095 .170 TYPE 310 

0096 310 FORMAT( ' XI,X2»RAD1iRAD2,NDIVRrNDlVX*NSEG?'> 

0097 READ(5r*rERR=170) XI rX2»RAD1»RAD2rNDIVR r NDIVX»NSEG 

C 

C LOOP FOR AXIAL DEPTH 
C 

0098 DO 320 1=1»NDIVX 

0099 X=Xl+<1-0.5)* <X2-X1)/NDIVX 

C 

C LOOP FOR RADIAL THICKNESS 
C 

0100 DO 320 M=lfNDIVR 

0101 R=RAD1 + (M- 0.5)* <RAD2-RAD1)/NDIVR 

0102 R=R*NSEG*8IN<PI/NSEG)/<PI*C0S<PI/NSEG)) 

C 

C LOOP AROUND CIRCUMFERENCE 
C 

0103 DO 320 J=1»NSEG 

0104 K=<(I~1>*NDIVR+M-1)*NSEG+J+IWIRE~1 

0105 ANG=PI*2./NSEG 

0106 CX < l<» 1) =X 

0107 CX<K»2)=X 

0108 CY < K»1)-COS(ANG* < J-1))*R 

0109 CY(K , 2)=CQS(ANG*J)*R 

0.110 CZ<K»1)-SIN(ANG*(J~1) )*R 

01.1.1 CZ < K r 2) -SIN ( ANG* J ) *R 

0112 320 CONTINUE 

0113 NEND=NSEG*NDIVX*NDIVR+IWIRE-1 

0114 ISYMM( ISETI »3>=ND.IVX*NDIVR 

C 

C END OF AXIAL CIRCULAR E/M 
C STORE E/M ELEMENT NO. 

C 

0115 220 ISYMM( ISETI r2)“NEND-IW'IRE+l 

C 

C OUTPUT SECTION 
C 

011A IF(I0UT0P.EQ.3) GO TO 435 
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0118 340 TYPE 350 

0119 350 FORMAT( ' E/M OUTPUT?') 

0120 READ(5 * # *ERR=340) IOUTOP 

0121 GO TO (435*360*370*435). (IOUTOP+1) 

0122 370 TYPE 380 

0123 TYPE 381 

0124 380 FORMAT(' WIRE COORDINATES') 

0125 381 FORMAT(5X * 2HX1 * 8X * 2HY1* 8X,2HZ1* 8X » 2HX2 * 8X * 2HY2 * 8X * 2HZ2 > 

0126 DO 390 I=IWIRE*NEND 

0127 390 TYPE 400*<CX<I *J)*CY(I *J)*CZ<I *J)*J=1*2) 

0128 400 FORMAT(6F10.5) 

0129 360 IF<INPOPT.NE.2) GO TO 435 

0131 TYPE 410 

0132 TYPE 411 

0133 410 FORMAT(' INTERFACE NO. AND COORDS') 

0134 411 FORMAT(3X,4HINT.*9X*1HX*14X*1HY*14X*1HZ) 

0135 DO 420 J-=l *NINT 

0136 420 TYPE 430*(J*XX(I *J)*YY<I *J)*ZZ<I *J)* 1=1*4> 

0137 430 FORMAT<I5*3F15.5) 

C 

C INCREMENT WIRE ELEMENT COUNTER 
C 

0138 435 IWIRE=NEND+1 

0139 100 CONTINUE 

C 

C END OF E/M SET LOOP 

C COMPLETE SYSTEM DATAFILE OUTPUT 

C 

0140 440 IF(IOUTOP.NE.3) GO TO 90 

0142 TYPE 50 

0143 CALL ASSIGN<2**-l*'NEW'* *) 

0144 DEFINE FILE 2 (1003*12*U*M) 

0145 DO 450 1=1*3 

0146 450 WRITE(2'I)(ISYMM(J*I)*J=1*5) 

0147 DO 460 K=1* 1000 

0148 460 WRITE(2'K+3) (CX(K*vl) *CY<K> J) *CZ(K* J) * J=1 *2) 

0149 CLOSE<UNIT=2> 

C 

C READ AHEAD E/M CURRENTS 
C 

0150 90 DO 470 ISETI=1*NSET 

0151 NOC=INOCO<ISYMM<ISETI* 1)> 

0152 TYPE 110*ISETI 

0153 490 TYPE 500,NOC 

0154 500 FORMAT (' E/M AMP. TIJRNS/1000? <',I2,'>') 

0155 READ(5 * * *ERR=490)(CURR(I *ISETI)*1 = 1 * NOC) 

0156 DO 470 1=1*NOC 

0157 CURR<I *ISETI)=CURR<I *ISETI)#1000* 

0158 470 CONTINUE 

0159 RETURN 

0160 END 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 


0009 

0010 

0011 

0012 

0013 

0014 

0016 

0017 

0018 

0020 

0021 


0022 

0024 

0025 


SUBROUTINE CALC<CXyCYyCZ,CORErDISP*ISYMM,CURR,IDISP) 
VIRTUAL CX < 1000 »2> r CY < 1000 f 2> f CZ<1000 f 2)f CORE(1000,6) 
DIMENSION POS(6>,BB<3),GB(9),COSR(3)fSINR(3),T(3),DPOS(3> 
DI MENS ! ON DISP ( 7 ) , FORCE < 6 > r I SYMM < 5 ,3 >, CURR <16 » 5)» I B I SP < 3 > 
PI-3*141592654 
XMU-0♦00000125664 
LUNOUT-5 
LUNFLD-5 


C OUTPUT ROUTING 

C OPEN DATAFILES IF REQUIRED 

C 

10 TYPE 20 

20 FORMAT(' OUTPUT OPTION?') 

READ(5,*,ERR-10) IOUT 
IFLIST-IOUT/10 
I OF ILE-1 OUT-1FL. I ST* 10 
IF(IOFILE.NE.l) GO TO 30 
LUNOUT-1 

OPEN(UNIT-1,NAME='2F0RCE.DAT',TYPE-'NEW' 
1,RECORDSIZE-78) 

30 IF<IELIST*NE.2) GO TO 40 

LUNFLD-2 

OPEN(UNIT-2,NAME-'2FIELD.DAT',TYPE-'NEW' 
1,RECORDSIZE-120) 


FORM-' FORMATTED' 

■ FORM-'FORMATTED' 


40 

60 


IF <IFL1ST.EQ•0.OR * LUNFLD.EQ.2) GO TO 50 
WRITE(LUNFLD,60) 

FORMAT(I HO,5X,lHX,9X,1HYr 9X,1HZ,8X,2HBX,8X,2HBY,8X,2HBZ,8X 
1,3HBXX, 7X, 31-IBXY, 7X ,3HBXZ, 7X, 3HDYY , 7X , 3HBYZ, 7X, 3HBZZ ) 


C 

C POSITION/ORIENTATION LOOP START 


C 

0026 50 


DO 70 K-0,IDISP(1) 


0028 80 FORCE(N)-0. 

C 

C CONVERT MODEL ORIENTATION 
C 

0029 DO 90 J-l,3 

0030 COSR(J)-COS(DISP(Jf3 >) 

003.1 90 SINR(J)=SIN < DISP(J+3)) 

C 

C MODEL. ELEMENT . LOOP START 
C 

0032 100 DO 110 N-l,IDISP<3) 

C 

C ROTATE ELEMENT 
C 

0033 DO 120 JJ-1,2 

0034 

0035 POS(J)-CORE(N,J)*COSR(1)*COSR(2)-CORE(N,J+1> *SINR(1) 

1*COSR(2)fCORE <NfJ+2) *SINR< 2) 

0036 POS(J+1)-CORE < N f J)* < SINR(1)*COSR(3)+COSR(1)#SINR< 2)* 

IS I NR (3 ) ) (CORE < N, J f 1) * (COSR (.1) *COSR < 3) -SI NR (1) *S (NR (2) #S I NR < 3 ) ) 
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2+CORE (N»J+2) *C0SR (2 > *SINR(3 > 

0037 120 POS( J+2)-COFit' (NrJ)*(SINR(1)*SINR(3)-COSR(1>*SINR(2>*COSR( 3> ) 

l+CORE(NrJ+l)*<C0SFi(l)*SINR(3)+SINR(l>*SINR(2>*C0SR(3>> 

2+CORE< N r J+2)*COSR(2)*COSR(3) 

C 

C GET MODEL ELEMENT ABSOLUTE POSITION 

C 

0038 DO 130 J-lr 3 

0039 130 DPOS < J)=POS(J)+DISP < J) 

C GET FIELD FOR THIS ELEMENT 

0040 CALL FIELD(BBrGBrCXrCYrCZrlSYMMrCURRrDPOS) 

C 

C COMPUTE ELEMENT FORCES 
C 

0041 DO .140 J=1 r3 

0042 JJ~< J-1>*3+1 

0043 T (J) - < POS (4)*GB(JJ > +POS < 5 > *GB <JJ+1) +POS (6) #GB (J J+2)) 

0044 140 FORCE(J)"FORCE(J)+T < J)/XMU 

C 

C AND ELEMENT TORQUES AND MOMENTS OF FORCES 
C 

004S FORCE(4>-FORCE (4 > + (( POS(2)*T(3)"POS(3)*T(2)) + 

1(P0S(5)*BB(3)~P0S(6)*BB<2>))/XMU 
0046 FORCE < 5)-FORCE(5 > + < < POS(3)*T(1)-POS <1)*T(3)) + 

1< POS< 6)#BB(1)-POS(4)*BB(3)))/XMU 
0047 FORCE(6)=FORCE(6 > + (< POS(1> *T< 2)-POS(2)*T <1) ) + 

1< POS(4)*BB(2>-POS < 5)*BB(1)>)/XMU 
C 

C FIELD OUTPUTS? 

C 

0048 IF(IFLIST.EQ.O) GO TO 110 

0050 WRITE(LUNFLDr150) POS(1)r POS < 2)rPOS < 3 > r BB <1>»BB < 2 > r BB < 3) 

lrGBCI.) rGB(2> r6B(3) rGB(S) rGB(6) rGB<9) 

0051 150 FORMAT <12F10.4) 

0052 110 CONTINUE 

C 

C END OF ELEMENT LOOP 
C FINAL OUTPUTS 
C 

0053 IF<K•GT.0.OR.IOFILE«EQ.1) GO TO 160 

0055 WRITE<LUNOUT > 170) 

0056 170 FORMAT(7Xr6HF0RCEXr 7X r 6HF0RCEYr 7X r 6HF0RCEZr 7X r 6HT0RQUX r7X 

1 1 6HT0RQUY r 7X >6HT0RQUZ) 

0057 .160 WRITE (LUNOUT r ISO) (FORCE (I)rl = lr6) 

0058 180 FORMAT(6F13.2) 

C 

C INCREMENT DISPLACEMENT 
C 

0059 DISP(IDISP<2))=DISP(IDISP(2)>+DISP(7> 

0060 70 CONTINUE 

C 

C END OF DISPLACEMENT LOOP 
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0061 

0063 

0064 190 

0066 

0067 200 

0068 


IF(LUNOUT.NE.l) GO TO 190 
CLOSE (UNITED 
IF(LUNFLD.NE•2) GO TO 200 
CLOSE<UNIT=2> 

RETURN 

END 


FORTRAN IV V02.5 PAGE 001 

0001 SUBROUTINE FACE(NINT»M»N*WW,W,IWIRE> 

C THIS SUBROUTINE DERIVED FROM MIT PROGRAM "TABLE"» BUT WITH TWO 
C SIGNIFICANT ALTERATIONS IN FUNCTION. REFER TO PROGRAM DESCRIPTION 
C 

0002 VIRTUAL W(1000»2) 

0003 DIMENSION WW<4*10> 

0004 IND--IWIRE-1 

C 

C GENERATE FOURTH CORNERS OF INTERFACES 


0005 

0006 


0007 

0008 

0009 

0010 

0011 


0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


C 

DO 10 I**1,NINT 

10 WW< 4»I>~WW<3*I>-WW<2»I)+WW<1 tl) 

C 

C GENERATE INTERMEDIATE COORDINATES FOR EACH INTERFACE 
C 

DO 20 1=1 !>M 
DO 20 J”1 r N 
DO 20 K"-l fNINT 
IND-IND+l 

W( INDr 1 ) = (WW( 1 rK>#(2)|((N~J)+1 )#(2#(M-I > + l) +WWC4 »K (2&J--1) 
l*<2#<M-I) + l)+WW(2»K>*<2#I-l>*<2*<N--J> + l)+WW<3rlO)|c<2*I--l) 
2*<2*J-1> )/<4*M*N) 

20 CONTINUE 

C 

C GENERATE W2 FROM W1 SO THAT ENDPOINTS OF WIRES ARE TOGETHER 
C 

DO 30 I-l f M*N 

J==<I--1)*NINT 

Wa*NINTr2)=W<JH»l> 

DO 30 K~1>NINT-1 
W(J+Kf2)~W(J+K+l»1) 

30 CONTINUE 

RETURN 
END 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 


0011 

0012 

0013 

0015 

0016 

0018 


SUBROUT INE FIELD(BB> GB»CX,CY r CZ >ISYMM »CURR»DPOS> 
VIRTUAL CX<1000>2)»CY<1000»2>>CZ<1000»2> 

BIMENSION S(3) t T(3)»DC(3) t DD(3)>DG<3)>DP< 3)> DS(3) 
DIMENSION BB<3)»GB<6)rCURR<16»5) > ISYMMC5»3) r DP0S<3> 
IWIRE™! 

PI“3.141592654 
DO 10 N=1 >9 
10 GB(N)=0. 

DO 20 N-l>3 
20 BB(N)=0. 

C 

C COIL. SET LOOP START 
C 

DO 30 ISET-1»5 
INCOIL-O 

IF(ISYMM(ISETc1)«E0.0) GO TO 30 
IREFL=1 

IF< ISYMM( ISET»1) »l.T*7) GO TO 40 
IREFL---1 
C 

C CALCULATION VECTORING AND CONTROL SECTION 


0019 40 DO 50 IX-"1»I REEL*-2 


0020 

0021 

0022 


0025 

0026 

0027 

0029 


0031 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 


C MASTER COIL 
C 

DY=1. 

DZ~1. 

DYZ-0* 

INCOIL"INCOIL+1 
ICURR-IX 

ASSIGN 60 TO IPOINT 
GO TO 1000 

60 IF(ISYMM(ISET»1)*E0»1♦OR*ISYMM(ISET»1)»EQ»7) GO TO 50 

IF<ISYMM(ISET»1).EQ.2.OR♦ISYMM<ISET »1).EO♦8) GO TO 70 
C 

C ROTATIONAL SYMMETRIES 
C 

IF<ISYMM<ISET,1>.GT.A) GO TO 80 
I ROT-- (I SYMM <ISET»1> ~2 > *2 
GO TO 90 

80 IROT”(ISYMM(ISET »l>-8)#2 

90 DO 100 N“.1. * IRGT-1 
ANGL--2. #PI *N/1 ROT 
DY=COS(ANGL> 

DZ--DY 

DYZ”SIN(ANGL) 

INCOIL=INCOIL+1 
ICURR”IX 

ASSIGN 100 TO IPOINT 
GO TO 1000 
100 CONTINUE 
GO TO 50 
C 
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C REFLECTION SYMMETRIES 
C 

0047 70 bz=~i. 

0048 DY=1. 

004V DYZ=0. 

0050 INC0IL=INC0IL+1 

0051 ICURR--IX 

0052 ASSIGN 110 TO IPOINT 

0053 GO TO 1000 

0054 .110 DY---1. 

0055 DZ=1. 

0056 INC0IL"=INC0IL+1 

0057 ICURR--IX 

0058 ASSIGN 120 TO IPOINT 

0059 GO TO 1000 

0060 120 DY--1• 

0061 DZ=~1. 

0062 INC0II..-INC0IL+1 

0063 ICURR-IX 

0064 ASSIGN 105 TO IPOINT 

0065 GO TO 1000 

0066 105 CONTINUE 

0067 50 CONTINUE 

C 

C INCREMENT WIRE COUNTER 
C 

0068 IWIRE=IWIRETISYMM <ISET,2 > 

0069 30 CONTINUE 

C 

C END OF COIL SET LOOP 
C COMPLETE DUMMY VARIABLES 
C 

0070 GB(4-)=GB<2) 

0071 GB ( 7) “-GB (3 > 

0072 GB<8)=GB<6) 

0073 RETURN 

C 

C FIELD AT POINT SECTION 
C 

C CORRECT WIRE CURRENTS 
C 

0074 1000 CURREN-1CURRflCURR(INCOIL 1 1SET)/1SYMM <ISET 1 3) 

C 

C Jtc******^*****#****************#************^***^***^****^^*** 
C THE REMAINDER OF THIS SUBROUTINE IS TAKEN WITH NO SIGNIFICANT 
C ALTERATION FROM THE PROGRAM "TABLE* t COURTESY MASSACHUSETTS 
C INSTITUTE OF TECHNOLOGY. 

C 

C WIRE ELEMENT LOOP START 
C 

0075 DO 2000 L.==IWIRE , IWIRE+ISYMM< ISETi-2)-1 

C 
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C CALCULATE ArBrCrDrErF 
C 

0076 A=<IX*CX<L»1)'-DP0S<1>> 

0077 B=<IX*CX<L»2)-DP0S(1)> 

0078 C~(DY*CY(L,1)-DYZ*CZ< L»1)-DPOS< 2)) 

0079 D= < DY*CY (L> 2 > -DYZ#CZ (L r 2) -DPOS < 2)) 

0080 E = (DZ*CZ <Lr1) +DYZ*CY< L. r 1 >-DPOS< 3 > > 

0081 F=(DZ*CZ(L ,2 )+DYZ*CY< L r2)-DPOS(3)) 

C 

C SUBSCRIPT ArBrCrDrErF FOR LATER USE 
C 

0082 S (.1) - A 

0083 S(2)-C 

0084 S(3)™E 

0085 T(1)=B 

0086 T < 2)=D 

0087 T(3)=F 

C 

C CALCULATE UrVrU 
C 

0088 U=C#F-D*E 

0089 V«E*B~F*A 

0090 W-A#D™B*C 

C 

C CALCULATE RHOlr RH02 AND THEIR DOTr CROSS PRODUCTS 
C 

0091 R1-(A*A+C*C+E*E)**0.5 

0092 R2-(B*B+D*D+F*F)**0.5 

0093 R8-R1+R2 

0094 RM=R1*R2 

0095 RDR=A*B+C*D+E*F 

0096 RXR=U+V+U 

C 

C CALCULATE THE DERI OAT IVES OF THE StJMr ETC. OF RH01 AND RII02 
0097 DO 1010 M“.L»3 

0098 DP<N)=-(S < M)*R2/R1IT(M)*R1/R2) 

0099 DS (M)---(S(M) /R1+T < M > /R2 ) 

0100 DD(M)~~ < S(M)+T< M)) 

0101 1010 CONTINUE 

0102 DC <1)"F-E+C-D 

0103 DC < 2)-E-F+8-A 

0104 DC(3)=D-C+A~B 

C CALCULATE AND TEST H TO DETERMINE EON. FOR G TO BE USED 
C 

0105 H-(RM+RDR)/RM 

0106 IF(H.LT.O.Ol) GO TO 1020 

0108 G-RS/(RM*(RM+RDR)) 

C 

C CALCULATE G AND ITS XrYrZ DERIVATIVES 
C 

0109 DO 1030 M-l»3 

0110 DGA-RM*(RM+RDR)*DS < M) 
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011 1 DGB--RS*(RM*(DP<M) +DD(M>> +DP< M ) * <RM+RDR) > 

01 1 2 1030 DG < M) - (DGA-DGB)/(RM*C RM+RDR))**2 

0113 GO TO 1040 

0114 1020 G=<(RS)*<RM-RDR))/<RM*RXR*RXR) 

0115 DO 1050 M=1>3 

011 6 DGA=(RS*(DP(M)-DD< M))+DS(M)*<RM-RDR)>*RM*RXR**2 

011 7 DGB=RS* ( RM-RDR) * (RM#2. *RXR*DC (M ) +DP < M ) *RXR**2> 

0118 1050 DG<M)=<DGA-DGB)/<RM*RXR**2>**2 

C CALCULATE THE FIELD AND GRADIENT CONTRIBUTIONS OF EACH 
C CURRENT ELEMENT 

C 

0119 1040 DGX-DG(l) 

0120 DGY-DG<2> 

0121 DGZ=DG<3> 

0.122 CURP=CURREN/10000000* 

0123 CURM-CURP*G 

C 

C SUM THE INDIVIDUAL CONTRIBUTIONS TO THE FIELD AND GRADIENT TO 
C GET THE TOTAL. FIELD AND GRADIENT 


0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 


2000 


0.134 

0.135 


BB(1)«BB<1)+Cl)RM*U 

BB < 2)-BB(2 >+CURM*V 

BB < 3)-BB < 3)+CURM*W 

GB(1)-GB<1)+CURP*U#DGX 

GB (2 ) -GB < 2 ) +CL.JRP* < G* <E-F ) +U*D6Y ) 

GB(3)-GB(3)+CURP*(G*(D-C)+U*DGZ) 

GB(5)-GB(5)+CURP*V*DGY 

GB(6)-GB( 6 )+CURP*(G*< A-B)+V*DGZ) 

GB(9>-GB(9)+CURP#DGZ#U 

CONTINUE 

OF WIRE ELEMENT LOOP 

GO TO IPOINT 
END 
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VARIABLE LIST FOR SUBROUTINE MODIN ********************* 
I» INTO MOBIN 

D OUT OF MOBIN 

SPfIDISP LOADED WITH MODEL CHARACTERISTICS 


MODEL ELEMENT POSITION AND MAGNETIZATIONS 
X POSITION 


X WISE HAGNETXZATION*ELEMENT VOLUME 


REQUIRED MODEL DISPLACEMENT FROM ZERO POSITION 
X DISPLACEMENT 


YAW DISPLACEMENT 


STEP SIZE FOR DISPLACEMENT LOOP WHERE REQUIRED 
LOOP CONTROL PARAMETERS 

NO. OF REQUIRED STEPS IN POSITION/ORIENTATION LOOP 
PARAMETER (IN ARRAY DISP> TO BE INCREMENTED 
TOTAL NO. OF MODEL ELEMENTS 
MODEL INPUT OPTION 
INTO 

FUSELAGE INPUT OPTION 

Nl'MBFR OF MODEL ELEMENTS (PASSED TO IDISP<3>> 




DISPLACEMENT 



‘rincipal variable lisl 





HALF WING CENTRE GAP <FUSELAGE RADIUS?) 

WING ROOT CHORD 

WINS thickness 

TIP CHORD/ROOT CHORD 
WING surer‘BACK 

X WISE DISPLACEMENT OF WING FROM MODEL DATUM 
2 WISE DISPLACEMENT OF WING FROM MODEL DATUM 

(DLONG.BVERT ONLY USEFUL WHERE FUSELAGE PRESENT) 
NO* OF ELEMENTS ALONG SINGLE WING PANEL 
NO. OK ELEMENTS ALONG WING CHORD 
NO. OF ELEMENTS THROUGH WING THICKNESS 



Ft EMFNT VOLUME 

EQUIVALENT VALUE OK TAPER FOR SINGLE ELEMENT 





PRINCIPAL VARIABLE LIST FOR SUBROUTINE COIL *W******>mr-m>m)m**## 


VARIABLES PASSED INTO COIL 
NONE 

VARIABLES PASSED OUT OF COIL 

CX,CY,CZ,ISYMM»CURR e/m array data 


CX,CY,CZ X,Y,Z COORDINATES OF WIRE ELEMENTS 

(N,l) WIRE ELEMENT START 

<N,2) WIRE ELEMENT END 

ISYMM E/M SYMMETRY CODES AND SET DATA 

(N, 1 > SYMMETRY OPTION FOR THIS SET 
(N,2> NO, OF WIRE ELEMENTS PER E/M IN THIS SET 

(N,3) NO. OF LOOPS OF ELEMENTS IN THIS SET 

CURR E/M TOTAL AMPERE-TURNS 

(1 ,N) IN FIRST E/M OF SET 

(2»N) SECOND 

ETC. ETC. 

XX,YY,ZZ COORDINATES OF E/M INTERFACES (WHERE USED) 

(i ,N) FIRST CORNER 

(2,N> SECOND 

(3, N > THIRD 

<4,N> FOURTH (PROGRAM GENERATED) 

INOCO NO. OF E/MS PER SET 

<N> NO. IN SET WITH SYMMETRY OPTION N 
NSET NO. OF E/M SETS 

ISETI POINTER TO E/M SET DATA 

IWIRE POINTER TO FIRST WIRE ELEMENT OF E/M 

NEND POINTER TO LAST WIRE ELEMENT OF E/M 

IOUTOP E/M OUTPUT OPTION 

INPOPT E/M INPUT OPTION 

NOC NO. OF E/MS IN THIS SET 

SOME VARIABLES ACTIVE IN THE SINGLE LOOP INPUT SECTION 
NINT NO. OF ELEMENTS 

SOME VARIABLES ACTIVE IN THE INTERFACE ARRAY INPUT SECTION 

MM,NN NO. OF LAYERS »LOOPS/LAYER OF WIRE ELEMENTS 

NINT NO. OF INTERFACES 

SOME VARIABLES ACTIVE IN THE LATERAL CIRCULAR E/M INPUT SECTION 

DX»DY X,Y WISE DISPLACEMENT OF ORIGIN OF E/M AXIS 

RADI,RAD2 INNER,OUTER RADII OF E/M 

DY1,DY2 NORMAL DIST. FROM X AXIS TO E/M FRONT,REAR FACE 

NDIVR,NDIVY NO. OF RADIAL,AXIAL LOOPS 

NSEG NO. OF WIRE ELEMENTS/LOOP (CIRCUMFERENTIALLY) 

ANGLE ANGLE FROM Y AXIS (IN YZ PLANE) TO E/M AXIS 

K COUNTER OF WIRE ELEMENTS IN THIS E/M 

R RADIUS FROM E/M AXIS TO END OF WIRE ELEMENT 

SOME VARIABLES ACTIVE IN THE AXIAL CIRCULAR E/M INPUT SECTION 
XI,X2 AXIAL DISTANCE TO E/M INNER,OUTER FACE 

NDIVR,NDIVX NO. OF RADIAL,AXIAL LOOPS 

NSEG,RADI, 

PAD2»R,K AS LATERAL CIRCULAR E/M SECTION 
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PRINCIPAL VARIABLE LIST FOR SUBROUTINE CALC :m**Jm***)m*:mm*#**#**# 


VARIABLES PASSEB INTO CALC 

CX »CY »CZ»ISYMM»CURR FROM COIL 

DISPfIDISP»CORE FROM MODIN 

VARIABLES PASSED OUT OF CALC 
NONE 


XMU PERMEABILITY OF FREE SPACE(SI) 

LUNOUT FORTRAN LOGICAL UNIT NUMBER FOR MAIN OUTPUTS 

LUNFLD FORTRAN LOGICAL UNIT NUMBER FOR FIELD LISTINGS 

IOUT OUTPUT OPTION 

DECODED INTO 

IFLIST =2 FIELD LISTING TO DATAFILE 

=1 FIELD LISTING ON CONSOLE DEVICE 

=0 NONE 

IOFILE =1 MAIN OUTPUTS TO DATAFILE 

=0 MAIN OUTPUTS TO CONSOLE DEVICE 

K DUMMY VARIABLE CONTROLLING POSITION/ORIENTATION LOOP 

FORCE MAIN RESULTS ARRAY 

(1) FORCE IN X DIRECTION 

(2) Y 

(3) Z 

(A) TORQUE ABOUT X AXIS 

<5) Y 

<£) Z 

COSR COSINES OF MODEL ROTATIONS 

<1) YAW 

<2> PITCH 

<3> ROLL 

SINR SINES OF MODEL ROTATIONS 

(1) YAW 

< 2) PITCH 

(3) ROLL 

N MODEL ELEMENT NUMBER 

POS POSITION AND MAGNETIZATION COMPONENTS OF MODEL ELEMENT 

(1) X POSITION 

(2) Y 

<3) Z 

<4) X WISE MAGNETIZATION 

<5) Y 

( 6 ) 2 

DPOS ABSOLUTE POSITION OF TRANSLATED MODEL ELEMENT 

(1 > X 

<2> Y 

<3) Z 
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PRINCIPAL VARIABLE LIST FOR SUBROUTINE FIELD 


VARIABLES PASSED INTO FIELD 
CXnCY tCZ> ISYMM»CURR 


DPDS 


FROM COIL 
FROM CALC 


VARIABLES PASSED OUT OF FIELD 
BB.GB 


FIELD AT REQUIRED POINT 


BB 

< .1) 
(2) 
( 3 ) 
GB 

( 1 ) 
<2> 
< 3 ) 
< 4 > 
<S) 
< 6 ) 
( 7 ) 
< 8) 
( 9 ) 


I SET 
INCOIL 
I WIRE 
I REEL 


BXrDYjDYZ.IX 

I ROT 

ICURR 

CURREN 

IPOINT 


TOTAL FIELD AT POINT (TESLA) 

X COMPONENT 
Y 
2 

FIELD GRADIENTS AT POINT (TESLA/METRE) 

GRADIENT OF X FIELD IN X DIRECTION 
X Y 

X Z 

Y X 

Y Y 

Y Z 

Z X 

Z Y 

Z Z 

BY FREE SPACE FIELD EQUATIONS* 

( 4 ) “ < 2 ) > < 7 )-( 3 ) , ( 8 )=-< 6 ) 
VECTOR FOR E/M SET DATA IN ISYMM* ETC. 

VECTOR FOR E/M DATA IN CURR 

VECTOR FOR WIRE ELEMENT DATA IN CX» ETC. 

AUXILIARY SYMMETRY CODE 

NO REFLECTION IN YZ PLANE 

REFLECTION IN YZ PLANE 

CONTROL PARAMETERS FOR REFLECTIONS AND ROTATIONS 
NO. OF E/M ROTATIONS ABOUT X AXIS 
CURRENT DIRECTION CORRECTION 
CURRENT IN WIRE ELEMENT 

RETURN VECTOR FROM MIT DERIVED SUBPROGRAM 


VARIABLES SPECIFIC TO MIT DERIVED SUBPROGRAM 
NOT CLASSIFIED 




.RUN 2F0RCE 

MODEL GEOMETRY OPTION? 

CYLIND.FUSE•S FLEN* FRAP *NA * NR? 

I.0*0.1.3*4.1 

SLAB UING:1/2SPAN * 1/2GAP * CHORD* THICN .TAPER.SWEEP? 



.NT.IUIOIR? 



O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO. 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 



O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 



ISH)1 Of-1 ION*SYNhiIKY option? 

3.10 

I»X • OY. RAD1*RAD2 .PY1 . PY 2 .NOIVRrNDIUY * NSEC* ANGLE? 
2.0*0..1.0.1 .5*1.5*2.0*2 .2*8. 43. 

E/M OUTPUT? 



























B etc. 

xy 

BXY BXZ BYY BYZ BZZ 

-0.0021 -0.1844 0.0072 0.0018 -0.0122 

-0.0000 -0.1873 -0.0000 0.0000 -0.0000 

0.0021 -0.1844 0.0072 -0.0018 -0.0122 

-0.0000 -0.1808 0.0069 0.0000 -0.0122 

-0.0000 -0.1870 -0.0002 -0.0000 0.0004 

-0.0000 -0.1369 0.0072 -0.0000 -0.0115 

0.0021 -0.1844 0.0072 -0.0018 -0.0122 

-0.0000 -0.1873 0.0000 -0.0000 -0.0000 

-0.0021 -0.1844 0.0072 0.0013 -0.0122 

-0.0000 -0.1869 0.0072 0.0000 -0.0115 

-0.0000 -0.1870 -0.0002 0.0000 0.0004 

0.0000 -0.1808 0.0069 -0.0000 -0.0122 

0.0004 -0.1891 0.0000 -0.0004 -0.0001 

-0.0004 -0.1891 0.0000 0.0004 -0.0001 

0.0006 -0.1939 0.0000 -0.0006 -0.0001 

-0.0006 -0.1939 0.0000 0.0006 -0.0001 

0.0007 -0.2003 -0.0000 -0.0007 -0.0001 

-0.0007 -0.2003 -0.0000 0.0007 -0.0001 

-0.0004 -0.1891 0.0000 0.0004 -0.0001 

0.0004 -0.1891 0.0000 -0.0004 -0.0001 

-0.0006 -0.1939 0.0000 0.0006 -0.0001 

0.0006 -0.1939 0.0000 -0.0006 -0.0001 

-0.0007 -0.2003 -0.0000 0.0007 -0.0001 

0.0007 -0.2003 -0.0000 -0.0007 -0.0001 

TORQUZ 

0.00 















A2.9 Diagnostics 


RT-11 FORTRAN generates WARNING diagnostics for the program 
segments MODIN, COIL and FIELD. These may be ignored 


Diagnostics for program unit MODIN 

Warning in line 0026 Possible modification of index "I" 

Warning in line 0045 Possible modification of index "I" 

Diagnostics for program unit COIL 

Warning in line 0023 Possible modification of index "I" 

Warning in line 0144 Possible modification of index "I" 

Diagnostics for program unit FIELD 

Warning in line 0012 DO loop entry at label "60" 

Warning in line 0019 DO loop entry at label "60" 

Warning in line 0036 DO loop entry at label "100" 



APPENDIX 3 


Scaling of results to other physical sizes 

The majority of the theoretical analysis presented herein is 
unspecific as to the choice of physical scale. However, numerical 
computations have generally been made at particular fixed scales 
for convenience. Scaling of computed performance is thus of interest 
and is relatively straightforward provided some care is taken. 

Maximum force and torque capabilities 

Exact magnetostatic equivalence of the model (any core material) 
will only be obtained if the model's magnetic characteristics and the 
applied field remain equal at all corresponding points during scaling. 
The former condition may be taken for granted, the latter generally 
requires adjustment of E/M characteristics. With consistent E/M 
geometry and characteristics (notably current density) the field 
strength at geometrically similar points varies as: 

H a Scale 

Thus to preserve field equivalence, the E/M current density, for 
instance, could be varied as: 

J otScale 

With such a condition it is clear that for fixed magnetic 
configurations: 

Model forces a Scale 2 (Vol.* Field gradient) 

Model torques a Scale 3 (Vol.* Field) 

This conflicts with the commonly quoted "Square-Cube" law (66 
and numerous others), which implied that magnetic force capability 
varied as core volume (Scale 3 ) whereas aerodynamic force varied as wing 
2 

area (Scale ). With regard to the maximum force and torque capability 
of MSBSs, the relations presented above (the "Square-Square" law?) 
are undoubtedly valid where permanent magnet cores are used, since here 
the performance of the core is limited by demagnetization, which in 
turn depends only on the total local field strengths and is not other¬ 
wise scale dependent. The relevance for soft iron cores is less 
obvious, since there is no particular restriction to operation of the 
core beyond saturation. However, when the core is far into saturation 
there is some evidence to support the validity of the Square-Square law 



as a measure of maximum capability, for instance Section 5.5.1, 
where torque capability was found to vary as Scale 3 , implying of 

2 

course that force capability may vary as Scale . 

The absolute maximum capability for iron cores may prove to be 
unroalisable with existing technology E/Ms, as appears to be the 
cast' for the SIM computations in Section 5, thence the practical 
scaling law may lie between the Square-Cube and the Square-Square 
laws, inclusive. 

General 

Where maximum capabilities, core demagnetization etc. are not 
critical factors, the following conditions are easily identified 
(15, 49). 

For consistent balance geometry, model magnetizations, E/M 
current densities, etc.: 

Model forces a Scale 3 

4 

Model torques a Scale 

'Technological limitations on E/M performance are functions of 
scale, in particular the usable current density falls with increasing 
scale, but the trends are not continuous, apparently involving 
enforced abandonment of particular conductor technologies at specific 
limiting scales. Further treatment of these effects cannot be 
attempted here. 
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APPENDIX 4 


M agnetic units and definitions 

All equations are given in SI units, whence: 

13 = y JH in free space (y = 4 x 10 7 Hm 1 : B in Tesla; 

° 0 H in A/m) 

However, two alternative subsystems exist, the Kennelley and 
Sommerfeld systems. The Kennelley system is used in this report, 

in the belief that it is somewhat more convenient where permanent 
magnetic material is present. The key definitions of this system are: 

B = y Q + J (flux through permanent magnet material) 

5f = J . Vh<5v (force on a dipole) 

6t = J_ x H_ 6v (torque on a dipole) 

Both formulations lead acceptably to the classical: 

B = Vr H J = y oXm H y r = 1 + X m 

In this report J is replaced by M where the equations are relatively 
unspecific regarding choice of material (permanent magnet or iron core). 
This is a non-standard symbology since the Sommerfeld system leads to 
(introducing primes to distinguish M' (traditional, in A/m) from M 

(Tesla) ): 

B = y (H + M' ) 

It is clear that M in fact equals y^M/. The use of M partially 
avoids the unfortunate clash between the standard symbology for 
polarization (J) and current density (J). There should be no con¬ 
textual difficulty in distinguishing the two otherwise. 

The first subscript to any field property describes the component 
under consideration, the second (where present) describes the gradient 
direction, e.g: 

= field strength in direction a 

H ab s gradient of in the b direction 

In free space H , = H, 
ab ba 

Occasionally a gradient of a gradient is required, this being 
given thus: 

(H , ) s the rate of change (gradient) of the H , 
ab c ab 

component taken in the c direction. 
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APPENDIX 5 


Characteristic symmetry of classical applied fields 

It happens that classical applied fields in MSBSs with high levels 
of E/M symmetry may be identified and classified by the character¬ 
istic symmetry of those fields taken about the origin of axes and, 
strictly, in balance axes. This method of representing fields has 
been used extensively in Section 4 to derive relatively generalised 
expressions for the magnetic couplings occuring with simple SPM wing 
cores. 

When a model is displaced or rotated away from the usual datum 
position and orientation, the strict conditions of applied field 
symmetry break down. However, the adjustments to E/M currents that 
would be required to maintain the principal components of applied 
field, hence the generated force and torque vectors, as fixed in 
model axes, will tend to approximately preserve the symmetry 
conditions cited herein, also in model axes. Analysis based on 
applied field symmetry thus at least holds good for small displacements 
and rotations of the model. 

The E/M configuration employed in this Appendix is identical to 
that used in Section 4 for the SPM roll coupling study, and readers 
should refer to this Section (Figs. 4.8 - 4.10) for precise 
definition of this geometry. However, the fundamental results remain 
valid for any configuration of vaguely conventional form and all 
geometries thereof, provided the E/M array incorporates the required 
E/M symmetry. Fig. A5.1 thus illustrates the general E/M configura¬ 
tion with Table A5.1 showing the definition of the applied fields for 
this configuration. 

Using FORCE, the field and field gradient components around 
the origin of axes (at 4^ 0.1m in x,y and z directions in fact) are 
easily found (Fig. A5.2). The respective signs of the field and 
field gradient components are of principal interest, the magnitudes 
being included here for information only. 

The applied field symmetry matrices may be deduced directly from 
Fig. A5.2 and represented as matrices of Odd (O) or Even (E) 
functions of x,y and z, following Section 4, and Equations 4.8 and 
4.9. The full derived matrices are presented in Fig. A5.3. 
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Fig. A5.1 Schematic diagram showing E/M configuration for field 

symmetry analysis. Refer to Figs. 4.8 - 4.10 for precise 
array geometry. 


Table A5.1 Definition of E/M current directi ons for classical 
applied fields. After Fig. A5.1 


Classical 

field 

Primary 

component 

h T 2 h *4 *5 X 6 *7 *8 *9 *10 

Lift 

Drag 

Sideforce ;1 
Pitch 

Yaw 

Magnetising 

(axial) 

B 

xz B 

B 

Xy B 

B 

y B 

X 

+ + 

I l 

I l 

+ + 

+ l 

l + 

l + 

+ i 

+ + 

I + 



I 11 r i2 I 13 I 14 

1 

Roll 

B 

YZ 

+ - + 


Current directions are chosen to make the classical primary field 
component for each applied field positive. 



LIFT 

FIELD 


Position (m) 

-0.1000 -0.1000 


z 

1000 

Field 

BX 

-0.0242 

component 

BY 

0.0002 

-0.1000 

-0.1000 

0 

1000 

0.0242 

-0.0002 

-0.1000 

0.1000 

-0 

1000 

-0.0242 

-0.0002 

-0.1000 

0.1000 

0 

1000 

0.0242 

0.0002 

0.1000 

-0.1000 

-0 

1000 

-0.0242 

-0.0002 

0.1000 

-0.1000 

0 

1000 

0.0242 

0.0002 

0.1000 

0.1000 


1000 

-0.0242 

0.0002 

0.1000 

0.1000 

0 

1000 

0.0242 

-0.0002 


(T) 

Field gradient component 

(T/m) 

BZ 

BXX 

BXY 

BXZ 

BYY 

-0.0244 

-0.0022 

-0.0016 

0.2434 

-0.0016 

-0.0244 

0.0022 

0.0016 

0.2434 

0.0016 

-0.0244 

-0.0022 

0.0016 

0.2434 

-0.0016 

-0.0244 

0.0022 

-0.0016 

0.2434 

0.0016 

0.0244 

0.0022 

-0.0016 

0.2434 

0.0016 

0.0244 

-0.0022 

0.0016 

0.2434 

-0.0016 

0.0244 

0.0022 

0.0016 

0.2434 

0.0016 

0.0244 

-0.0022 

-0.0016 

0.2434 

-0.0016 


PYZ 

-0.0016 

-0.0016 

0.0016 

0.0016 

0.0016 

0.0016 

-0.0016 


BZZ 

0.0037 

-0.0037 

0.0037 

-0.C037 

0.0037 

-0.0037 

0.0037 


1 DRAG 


FIELD 
(B ) 


-0.1000 -0.1000 

-0.1000 -0.1000 

- 0.1000 0.1000 

- 0.1000 0.1000 

0.1000 -0.1000 

0.1000 -0.1000 

0.1000 0.1000 

0.1000 0.1000 


Z BX 

-0.1000 -0.0226 

0.1000 -0.0226 

-0.1000 -0.0226 

0.1000 -0.0226 

-0.1000 0.0226 

0.1000 0.0226 

-0.1000 0.0226 

0.1000 0.0226 


BY 

0.0113 

0.0113 

-0.0113 

-0.0113 

0.0113 

0.0113 

-0.0113 

-0.0113 


BZ 

0.0113 

-0.0113 

0.0113 

-0.0113 

0.0113 

-0.0113 

0.0113 

-0.0113 


BXX BXY 

0.2265 -0.0007 

0.2265 -0.0007 

0.2265 0.0007 

0.2265 0.0007 

0.2265 0.0007. 

0.2265 0.0007 

0.2265 -0.0007 

0.2265 -0.0007 


BXZ BYY 

-0.0007 -0.1133 

0.0007 -0.1133 

-0.0007 -0.1133 

0.0007 -0.1133 

0.0007 -0.1133 

-0.0007 -0.1133 

0.0007 -0.1133 

-0.0007 -0.1133 



0.0002 

- 0.0002 

-0.0002 

0.0002 


BZ7 

-0.1133 
-0.1133 
-0.1133 
-0.1133 
-0.1133 
-0.1133 
-0.1133 
-0.1133 


SIDEFORCE 

FIELD 

(B ) 

xy ; 


-0.1000 
-0.1000 
-0.1000 
-0.1000 
0.1000 
0.1000 
0.1000 
0.1000 


-0.1000 

-0.1000 

0.1000 

0.1000 

-0.1000 

-0.1000 

0.1000 

0.1000 


-0.1000 
0.1000 
-0.1000 
0.1000 
-0.1000 
0.1000 
-0.1000 
0.1000 


BX 

-0.0242 

-0.0242 

0.0242 

0.0242 

-0.0242 

-0.0242 

0.0242 

0.0242 


-0.0244 

-0.0244 

-0.0244 

-0.0244 

0.0244 

0.0244 

0.0244 

0.0244 


0.0002 

- 0.0002 

- 0.0002 

0.0002 

- 0.0002 

0.0002 

0.0002 

- 0.0002 


BXX 

- 0.0022 

- 0.0022 

0.0022 

0.0022 

0.0022 

0.0022 

- 0.0022 

- 0.0022 


BXY BXZ 

0.2434 -0.0016 

0.2434 0.0016 

0.2434 0.0016 

0.2434 -0.0016 

0.2434 -0.0016 

0.2434 0.0016 

0.2434 0.0016 

0.2434 -0.0016 


BYY 

0.0037 

0.0037 

-0.0037 

-0.0037 

-0.0037 

-0.0037 

0.0037 

0.0037 


PYZ 

-0.0016 
0.0016 
-0.0016 
0.0016 
0.0016 
-0.0016 
0.0016 
-0.0016 


BZZ 

-0.00.16 

0.0016 

0.0016 

0.0016 

0.0016 

-0.0016 

-0.0016 


Fig. A 5 .2 Fiel d components around the origin of balance axes for 
classical appl i ed fields, following Fig. A5.1 









































PITCH 

FIELD 

( V 






ROLL 

FIELD 



-O.IOOO -0.1000 -0.1000 
-0.1000 -0.1000 0.1000 
-0.1000 0.1000 -0.1000 
-0.1000 0.1000 0.1000 
0.1000 -0.1000 -0.1000 
0.1000 -0.1000 0.1000 
0.1000 0.1000 -0.1000 
0.1000 0.1000 0.1000 


BX BY BZ 

0.0001 -0.0069 -0.0069 

-0.0001 0.0069 -0.0069 

-0.0001 -0.0069 0.0069 

0.0001 0.0069 0.0069 

-0.0001 -0.0069 -0.0069 

0.0001 0.0069 -0.0069 

0.0001 -0.0069 0.0069 

-0.0001 0.0069 0.0069 


BXX BXY BXZ 

-0.0008 -0.0008 -0.0000 

0.0008 0.0008 -0.0008 

0.0008 -0.0008 0.0008 

-0.0008 0.0008 0.0008 

-0.0008 0.0008 0.0008 

0.0008 -0.0008 0.0008 

0.0008 0.0008 -0.0008 

-0.0008 -0.0008 -0.0008 


byy byz bzz 

0.0004 0.0691 0.0004 
-0.0004 0.0691 -0.0004 
-0.0004 0.0691 -0.0004 
0.0004 0.0691 0.0004 
0.0004 0.0691 0.0004 
-0.0004 0.0691 -0.0004 
-0.0004 0.0691 -0.0004 
0.0004 0.0691 0.0004 


Fig. A5, 











































Fig. A5, 















Fig. A5.3 Derived symmetry matrices for classical applied fields 


Following Figs. A5.1, A5.2 

The field symmetry is presented as in Eqns. 4.8, 4.9 
LIFT . 



For the remaining fields the matrix alone is presented. 


DRAG SIDEFORCE 
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PITCH 


YAW 
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